/[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.7 by bottai, Tue Feb 7 11:38:52 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.6 2006/02/06 14:17:21 bottai 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  * Revision 3.6  2006/02/06 14:17:21  bottai
9  *  adapt again the resolution simulation to the ITRPB volume identification  *  adapt again the resolution simulation to the ITRPB volume identification
10  *  *
# Line 45  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 155  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           CALL GPRSPE(ANGLEX,ANGLEY,BMAGNET,XMULT,YMULT,DXPOS,DYPOS)              
180                CALL GPRSPE(ANGLEX,ANGLEY,BMAGNET,XMULT,YMULT,DXPOS,DYPOS)
          NXMULT(I)=INT(XMULT)  
          NYMULT(I)=INT(YMULT)  
181    
182           XAVSPE(I) = XAVSPE(I) + DXPOS              NXMULT(I)=INT(XMULT)
183           YAVSPE(I) = YAVSPE(I) + DYPOS              NYMULT(I)=INT(YMULT)
184                
185                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).LT.6) 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.6) 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.7  
changed lines
  Added in v.3.8
  ViewVC Help
Powered by ViewVC 1.1.23