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

Diff of /gpamela/gpspe/gpdspe.F

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

revision 3.4 by bottai, Wed Jan 18 09:41:33 2006 UTC revision 3.8 by cafagna, Sun Apr 9 23:29:19 2006 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gpdspe.F,v 3.3 2005/12/05 12:15:21 pam-ba Exp $  * $Id: gpdspe.F,v 3.7 2006/02/07 11:38:52 bottai Exp $
3  *  *
4  * $Log: gpdspe.F,v $  * $Log: gpdspe.F,v $
5    * Revision 3.7  2006/02/07 11:38:52  bottai
6    *  test of mail delivery
7    *
8    * Revision 3.6  2006/02/06 14:17:21  bottai
9    *  adapt again the resolution simulation to the ITRPB volume identification
10    *
11    * Revision 3.5  2006/01/24 13:52:01  bottai
12    * prepare to store the x and y multiplicity of the hit
13    *
14    * Revision 3.4  2006/01/18 09:41:33  bottai
15    *  Some correction to the spe plane ID due to the new geometry tree
16    *
17  * Revision 3.3  2005/12/05 12:15:21  pam-ba  * Revision 3.3  2005/12/05 12:15:21  pam-ba
18  * new spectrometer geometry and internal magnetic field  * new spectrometer geometry and internal magnetic field
19  *  *
# Line 19  Line 31 
31  *CMZ :  1.00/02 27/03/96  18.35.42  by  Francesco Cafagna  *CMZ :  1.00/02 27/03/96  18.35.42  by  Francesco Cafagna
32  *-- Author :    Francesco Cafagna   23/02/96  *-- Author :    Francesco Cafagna   23/02/96
33        SUBROUTINE GPDSPE        SUBROUTINE GPDSPE
34  ************************************************************************  *************************************************************************
35  *                                                                      *  *                                                                       *
36  * To perform SPE digitization                                          *  * To perform SPE digitization                                           *
37  *                                                                      *  *                                                                       *
38  * Called by: GUDIGI                                                    *  * Called by: GUDIGI                                                     *
39  * Author: Francesco Cafagna, 23/02/96 12.50.38                         *  * Author: Francesco Cafagna, 23/02/96 12.50.38                          *
40  *                                                                      *  *                                                                       *
41  ************************************************************************  *************************************************************************
42  #include "gcunit.inc"  #include "gcunit.inc"
43  #include "gpphit.inc"  #include "gpphit.inc"
44        PARAMETER (RADDEG=57.295779)        PARAMETER (RADDEG=57.295779)
# Line 36  c +SEQ,GPDSPE. Line 48  c +SEQ,GPDSPE.
48  #include "gprhit.inc"  #include "gprhit.inc"
49  #include "gpspe.inc"  #include "gpspe.inc"
50  #include "gpsed.inc"  #include "gpsed.inc"
51    #include "gpkey.inc"
52  *  *
53        INTEGER IS,ID,NV,NVM,GPNHIT,NTHIT,NHFRST        INTEGER IS,ID,NV,NVM,GPNHIT,NTHIT,NHFRST
54        CHARACTER*4 CSET        CHARACTER*4 CSET
# Line 146  c +SEQ,GPDSPE. Line 159  c +SEQ,GPDSPE.
159           PATHSPE(I) = VHITS(8,I)           PATHSPE(I) = VHITS(8,I)
160           IPARSPE(I) = VHITS(9,I)           IPARSPE(I) = VHITS(9,I)
161           P0SPE (I) = VHITS(10,I)           P0SPE (I) = VHITS(10,I)
162             IF(FFIELD.NE.0) THEN
163  C  C
164  C ADD RESOLUTION SIMULATION TO XAVSPE  C ADD RESOLUTION SIMULATION TO XAVSPE
165  C          C        
166           ANGLEX=0.              ANGLEX=0.
167           ANGLEY=0.              ANGLEY=0.
168           IF(ZOUTSPE(I).NE.ZINSPE(I)) THEN              IF(ZOUTSPE(I).NE.ZINSPE(I)) THEN
169              ANGLEX=RADDEG*ATAN( ABS(XOUTSPE(I)-XINSPE(I) )/                 ANGLEX=RADDEG*ATAN( ABS(XOUTSPE(I)-XINSPE(I) )/
170       +              ABS(ZOUTSPE(I)-ZINSPE(I) ) )       +              ABS(ZOUTSPE(I)-ZINSPE(I) ) )
171              ANGLEY=RADDEG*ATAN( ABS(YOUTSPE(I)-YINSPE(I) )/                 ANGLEY=RADDEG*ATAN( ABS(YOUTSPE(I)-YINSPE(I) )/
172       +              ABS(ZOUTSPE(I)-ZINSPE(I) ) )                 +              ABS(ZOUTSPE(I)-ZINSPE(I) ) )          
173           ENDIF              ENDIF
174           VPOS(1)=XAVSPE(I)              
175           VPOS(2)=YAVSPE(I)              VPOS(1)=XAVSPE(I)
176           VPOS(3)=ZAVSPE(I)                    VPOS(2)=YAVSPE(I)
177           CALL GUFLD(VPOS,BMAGNET)              VPOS(3)=ZAVSPE(I)    
178                CALL GUFLD(VPOS,BMAGNET)
179                
180                CALL GPRSPE(ANGLEX,ANGLEY,BMAGNET,XMULT,YMULT,DXPOS,DYPOS)
181    
182           CALL GPRSPE(ANGLEX,ANGLEY,BMAGNET,XMULT,YMULT,DXPOS,DYPOS)              NXMULT(I)=INT(XMULT)
183                NYMULT(I)=INT(YMULT)
184           XAVSPE(I) = XAVSPE(I) + DXPOS              
185           YAVSPE(I) = YAVSPE(I) + DYPOS              XAVSPE(I) = XAVSPE(I) + DXPOS
186                YAVSPE(I) = YAVSPE(I) + DYPOS
187    
188  C  C
189  C ADD THE SHIFT INDUCED BY THE MAG. FIELD  C ADD THE SHIFT INDUCED BY THE MAG. FIELD
190  C  C
191    
192           IF(ITRPB(I).EQ.1) THEN              IF(ITRPB(I).LT.6) THEN
193              XAVSPE(I) = XAVSPE(I) - 2.64e-4/4.5*BMAGNET(2)                 XAVSPE(I) = XAVSPE(I) - 2.64e-4/4.5*BMAGNET(2)
194           ELSE IF(ITRPB(I).EQ.2) THEN              ELSE IF(ITRPB(I).EQ.6) THEN
195              XAVSPE(I) = XAVSPE(I) + 2.64e-4/4.5*BMAGNET(2)                 XAVSPE(I) = XAVSPE(I) + 2.64e-4/4.5*BMAGNET(2)
196           ENDIF              ENDIF
197             ENDIF  
198        ENDDO        ENDDO
199    
200  *  *
201  * Loop on the hits retrived to fill the CWN common and the DIGI b  * Loop on the hits retrived to fill the CWN common and the DIGI b
202  *  *
Legend:
Removed from v.3.4  
changed lines
  Added in v.3.8
  ViewVC Help
Powered by ViewVC 1.1.23