/[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.1.1.1 by cafagna, Thu Jul 11 16:02:14 2002 UTC revision 3.3 by bottai, Tue May 2 12:01:08 2006 UTC
# Line 1  Line 1 
1  *  *
2  * $Id$  * $Id: gpuspe.F,v 3.2 2006/02/05 22:13:44 cafagna Exp $
3    *
4    * $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
9    * First GPAMELA release on CVS
10  *  *
 * $Log$  
11  *  *
12  *CMZ :  2.01/00 06/03/2000  13.07.03  by  Francesco Cafagna  *CMZ :  2.01/00 06/03/2000  13.07.03  by  Francesco Cafagna
13  *CMZ :  2.00/00 03/03/2000  15.39.06  by  Francesco Cafagna  *CMZ :  2.00/00 03/03/2000  15.39.06  by  Francesco Cafagna
# Line 9  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 37  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),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 65  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)
89             IF(NUMVOL(1).EQ.0) NVOLSPE(1)=6
90             NVOLSPE(2) = NUMVOL(2)
91             NVOLSPE(3) = NUMVOL(3)
92        ENDIF        ENDIF
93  *  *
94  * Track is exiting current volume  * Track is exiting current volume
# Line 107  Line 128 
128  * Store the hit  * Store the hit
129  *  *
130           IF(IDHIT.EQ.-1) THEN           IF(IDHIT.EQ.-1) THEN
131              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)
134              IF (IDHIT.EQ.0) THEN              IF (IDHIT.EQ.0) THEN
135                 WRITE(CHMAIL,10000) IS,ID,IT                 WRITE(CHMAIL,10000) IS,ID,IT
136                 CALL GMAIL(1,0)                 CALL GMAIL(1,0)
# Line 159  Line 182 
182  * Store the hit  * Store the hit
183  *  *
184              IF(IDHIT.EQ.-1) THEN              IF(IDHIT.EQ.-1) THEN
185                 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)
188                 IF (IDHIT.EQ.0) THEN                 IF (IDHIT.EQ.0) THEN
189                    WRITE(CHMAIL,10200) IS,ID,IT                    WRITE(CHMAIL,10200) IS,ID,IT
190                    CALL GMAIL(1,0)                    CALL GMAIL(1,0)
Legend:
Removed from v.3.1.1.1  
changed lines
  Added in v.3.3
  ViewVC Help
Powered by ViewVC 1.1.23