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

Diff of /gpamela/gpcal/gpucal.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:04 2005 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gpucal.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $  * $Id: gpucal.F,v 3.2 2003/12/17 11:32:48 pamela Exp $
3  *  *
4  * $Log: gpucal.F,v $  * $Log: gpucal.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 52  Line 55 
55        PARAMETER (CALIB=0.0001085200)        PARAMETER (CALIB=0.0001085200)
56  *END EM.  *END EM.
57        INTEGER IACT,IS,ID,IPA,IT,NUMVOL(20),IMEC,IDHIT        INTEGER IACT,IS,ID,IPA,IT,NUMVOL(20),IMEC,IDHIT
58  CML17/04/03      REAL DELOSS,PATH,TRAPAR(7),VHIT(NHCAL),XYZM(3),XYZD(3)        REAL DELOSS,PATH,TRAPAR(7),VHIT(NHCAL),XYZM(3),XYZD(3)
59        REAL DELOSS,PATH,TRAPAR(7),XYZM(3),XYZD(3)  *      REAL DELOSS,PATH,TRAPAR(7),XYZM(3),XYZD(3)
60        REAL DELTA2,SIGMA2,SIGMA,GRNDM0,GPGAUS,PATHINI        REAL DELTA2,SIGMA2,SIGMA,GRNDM0,GPGAUS,PATHINI
61        LOGICAL SAVE        LOGICAL SAVE
62  *  *
# Line 80  c         print*,'GPUCAL:last event=',ie Line 83  c         print*,'GPUCAL:last event=',ie
83  cccc            VHIT(4) = (DELOSS/CALIB) + VHIT(4)  cccc            VHIT(4) = (DELOSS/CALIB) + VHIT(4)
84  * END EM  * END EM
85  C ML  C ML
86  CML17/04/03            VHIT(4) = DELOSS + VHIT(4)              VHIT(4) = DELOSS + VHIT(4)
87              ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))=              ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))=
88       +           ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))+DELOSS       +           ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))+DELOSS
89  *END ML  *END ML
# Line 91  CML17/04/03            VHIT(4) = DELOSS Line 94  CML17/04/03            VHIT(4) = DELOSS
94  *  *
95        IF(IACT.EQ.1) THEN        IF(IACT.EQ.1) THEN
96           SAVE = .FALSE.           SAVE = .FALSE.
97  CML17/04/03         CALL VZERO(VHIT,NHCAL)           CALL VZERO(VHIT,NHCAL)
98  CML17/04/03         DO I=1,3           DO I=1,3
99  CML17/04/03            VHIT(I) = TRAPAR(I)              VHIT(I) = TRAPAR(I)
100  CML17/04/03         ENDDO           ENDDO
101           IF(IMEC.GT.0) THEN           IF(IMEC.GT.0) THEN
102              SAVE = .TRUE.              SAVE = .TRUE.
103  *EM:  *EM:
104  c            VHIT(4) = DELOSS/CALIB  c            VHIT(4) = DELOSS/CALIB
105  *END EM.  *END EM.
106  CML17/04/03            VHIT(4) = DELOSS              VHIT(4) = DELOSS
107  C ML:  C ML:
108              ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))=              ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))=
109       +           ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))+DELOSS       +           ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))+DELOSS
# Line 123  c         ENDIF Line 126  c         ENDIF
126  *EM:  *EM:
127  c            VHIT(4) = (DELOSS/CALIB) + VHIT(4)  c            VHIT(4) = (DELOSS/CALIB) + VHIT(4)
128  *END EM  *END EM
129  CML17/04/03            VHIT(4) = DELOSS + VHIT(4)              VHIT(4) = DELOSS + VHIT(4)
130  *ML:  *ML:
131              ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))=              ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))=
132       +           ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))+DELOSS       +           ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))+DELOSS
# Line 156  CML17/04/03            IF(VHIT(4).LT.0.) Line 159  CML17/04/03            IF(VHIT(4).LT.0.)
159           IF(IDHIT.EQ.-1) THEN           IF(IDHIT.EQ.-1) THEN
160  *ML:  *ML:
161  *            IF(SAVE) CALL GSCHIT (IS,ID,1,NUMVOL,VHIT,4,IDHIT)  *            IF(SAVE) CALL GSCHIT (IS,ID,1,NUMVOL,VHIT,4,IDHIT)
162  CL            IF(SAVE) CALL GSCHIT (IS,ID,1,NUMVOL,VHIT,1,IDHIT)              IF(SAVE) THEN
163                   CALL GSCHIT (IS,ID,1,NUMVOL,VHIT,1,IDHIT)
164    c               CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT)
165                ENDIF
166  *END ML.  *END ML.
167              IF (IDHIT.EQ.0) THEN              IF (IDHIT.EQ.0) THEN
168                 WRITE(CHMAIL,10000) IS,ID,IT                 WRITE(CHMAIL,10000) IS,ID,IT
# Line 180  CL            IF(SAVE) CALL GSCHIT (IS,I Line 186  CL            IF(SAVE) CALL GSCHIT (IS,I
186  *EM:  *EM:
187  c               VHIT(4) = (DELOSS/CALIB) + VHIT(4)  c               VHIT(4) = (DELOSS/CALIB) + VHIT(4)
188  * END EM.  * END EM.
189  CML17/04/03               VHIT(4) = DELOSS + VHIT(4)                 VHIT(4) = DELOSS + VHIT(4)
190  *ML:  *ML:
191                 ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))=                 ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))=
192       +              ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))+DELOSS       +              ENEMAT(NUMVOL(1),NUMVOL(2),NUMVOL(3))+DELOSS
# Line 212  CML17/04/03               IF(VHIT(4).LT. Line 218  CML17/04/03               IF(VHIT(4).LT.
218  *  *
219              IF(IDHIT.EQ.-1) THEN              IF(IDHIT.EQ.-1) THEN
220  *ML  *ML
221  *               IF (SAVE) CALL GSCHIT (IS,ID,1,NUMVOL,VHIT,4,IDHIT)  C               IF (SAVE) CALL GSCHIT (IS,ID,1,NUMVOL,VHIT,4,IDHIT)
222                   IF (SAVE) THEN
223                      CALL GSCHIT (IS,ID,1,NUMVOL,VHIT,1,IDHIT)
224                   ENDIF
225  CL               IF (SAVE) CALL GSCHIT (IS,ID,1,NUMVOL,VHIT,1,IDHIT)  CL               IF (SAVE) CALL GSCHIT (IS,ID,1,NUMVOL,VHIT,1,IDHIT)
226  *END ML.  *END ML.
227                 IF (IDHIT.EQ.0) THEN                 IF (IDHIT.EQ.0) THEN
Legend:
Removed from v.3.2  
changed lines
  Added in v.3.3
  ViewVC Help
Powered by ViewVC 1.1.23