/[PAMELA software]/gpamela/gpspe/gpuspe.F
ViewVC logotype

Diff of /gpamela/gpspe/gpuspe.F

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

revision 3.2 by cafagna, Sun Feb 5 22:13:44 2006 UTC revision 3.3 by bottai, Tue May 2 12:01:08 2006 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gpuspe.F,v 3.1.1.1 2002/07/11 16:02:14 cafagna Exp $  * $Id: gpuspe.F,v 3.2 2006/02/05 22:13:44 cafagna Exp $
3  *  *
4  * $Log: gpuspe.F,v $  * $Log: gpuspe.F,v $
5    * Revision 3.2  2006/02/05 22:13:44  cafagna
6    * A bug in the SPEctrometr hit definition fixed
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 12  Line 15 
15  *CMZ :  1.01/00 23/05/96  15.01.53  by  Francesco Cafagna  *CMZ :  1.01/00 23/05/96  15.01.53  by  Francesco Cafagna
16  *CMZ :  1.00/02 04/04/96  16.47.16  by  Francesco Cafagna  *CMZ :  1.00/02 04/04/96  16.47.16  by  Francesco Cafagna
17  *-- Author :    Francesco Cafagna   21/02/96  *-- Author :    Francesco Cafagna   21/02/96
18        SUBROUTINE GPUSPE(IACT,IS,ID,IPA,IT,DELOSS,PATH,        SUBROUTINE GPUSPE(IACT,IS,ID,IPAA,IT,DELOSS,PATH,
19       +                  TRAPAR,NUMVOL,IMEC)       +                  TRAPAR,NUMVOL,IMEC)
20  ************************************************************************  ************************************************************************
21  *                                                                      *  *                                                                      *
# Line 40  Line 43 
43  ************************************************************************  ************************************************************************
44  #include "gphit.inc"  #include "gphit.inc"
45  #include "gcunit.inc"  #include "gcunit.inc"
46    #include "gpgene.inc"
47    #include "gctrak.inc"
48  *  *
49        INTEGER IACT,IS,ID,IPA,IT,NUMVOL(20),NVOLSPE(3),IMEC,IDHIT        INTEGER IACT,IS,ID,IPAA,IT,NUMVOL(20),NVOLSPE(3),IMEC,IDHIT
50        REAL DELOSS,PATH,TRAPAR(7),VHIT(NHSPE)        REAL DELOSS,PATH,TRAPAR(7),VHIT(NHSPE)
51        REAL DELTA2,SIGMA2,SIGMA,GRNDM0,GPGAUS        REAL DELTA2,SIGMA2,SIGMA,GRNDM0,GPGAUS
52        LOGICAL SAVE        LOGICAL SAVE
53    
54    *ET , SB 27/01/06
55          IF(IMEC.GT.0.AND.DELOSS.GT.0.) THEN        
56             CALL GPUDIFFUSION(IACT,TRAPAR,NUMVOL,DELOSS,STEP,IPAA)
57          ENDIF
58    *end ET, SB  27/01/06
59    
60    
61    
62  *  *
63  * Track is inside a volume  * Track is inside a volume
64  *  *
# Line 68  Line 82 
82              VHIT(7) = DELOSS              VHIT(7) = DELOSS
83           ENDIF           ENDIF
84           VHIT(8) = PATH           VHIT(8) = PATH
85           VHIT(9) = IPA           VHIT(9) = IPAA
86           VHIT(10)= TRAPAR(7)           VHIT(10)= TRAPAR(7)
87           IDHIT = -1           IDHIT = -1
88           NVOLSPE(1) = NUMVOL(1)           NVOLSPE(1) = NUMVOL(1)
# Line 115  Line 129 
129  *  *
130           IF(IDHIT.EQ.-1) THEN           IF(IDHIT.EQ.-1) THEN
131  c$$$            IF(SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT)  c$$$            IF(SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT)
132    cv            IF(SAVE) WRITE(6,*) 'CALL GSAHIT'
133              IF(SAVE) CALL GSAHIT (IS,ID,IT,NVOLSPE,VHIT,IDHIT)              IF(SAVE) CALL GSAHIT (IS,ID,IT,NVOLSPE,VHIT,IDHIT)
134              IF (IDHIT.EQ.0) THEN              IF (IDHIT.EQ.0) THEN
135                 WRITE(CHMAIL,10000) IS,ID,IT                 WRITE(CHMAIL,10000) IS,ID,IT
# Line 168  c$$$            IF(SAVE) CALL GSAHIT (IS Line 183  c$$$            IF(SAVE) CALL GSAHIT (IS
183  *  *
184              IF(IDHIT.EQ.-1) THEN              IF(IDHIT.EQ.-1) THEN
185  c$$$               IF(SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT)  c$$$               IF(SAVE) CALL GSAHIT (IS,ID,IT,NUMVOL,VHIT,IDHIT)
186    cv               IF(SAVE) WRITE(6,*) 'CALL GSAHIT'
187                 IF(SAVE) CALL GSAHIT (IS,ID,IT,NVOLSPE,VHIT,IDHIT)                 IF(SAVE) CALL GSAHIT (IS,ID,IT,NVOLSPE,VHIT,IDHIT)
188                 IF (IDHIT.EQ.0) THEN                 IF (IDHIT.EQ.0) THEN
189                    WRITE(CHMAIL,10200) IS,ID,IT                    WRITE(CHMAIL,10200) IS,ID,IT
Legend:
Removed from v.3.2  
changed lines
  Added in v.3.3
  ViewVC Help
Powered by ViewVC 1.1.23