/[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.2 by cafagna, Mon Jul 25 11:53:21 2005 UTC revision 3.5 by bottai, Tue Jan 24 13:52:01 2006 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gpdspe.F,v 3.1.1.1 2002/07/11 16:02:13 cafagna Exp $  * $Id: gpdspe.F,v 3.4 2006/01/18 09:41:33 bottai Exp $
3  *  *
4  * $Log: gpdspe.F,v $  * $Log: gpdspe.F,v $
5    * Revision 3.4  2006/01/18 09:41:33  bottai
6    *  Some correction to the spe plane ID due to the new geometry tree
7    *
8    * Revision 3.3  2005/12/05 12:15:21  pam-ba
9    * new spectrometer geometry and internal magnetic field
10    *
11    * Revision 3.2  2005/07/25 11:53:21  cafagna
12    * Several updates. See history for details
13    *
14  * Revision 3.1.1.1  2002/07/11 16:02:13  cafagna  * Revision 3.1.1.1  2002/07/11 16:02:13  cafagna
15  * First GPAMELA release on CVS  * First GPAMELA release on CVS
16  *  *
# Line 23  Line 32 
32  ************************************************************************  ************************************************************************
33  #include "gcunit.inc"  #include "gcunit.inc"
34  #include "gpphit.inc"  #include "gpphit.inc"
35          PARAMETER (RADDEG=57.295779)
36        PARAMETER (NH=NHSPE,NHM=NHMSPE,NVOL=3)        PARAMETER (NH=NHSPE,NHM=NHMSPE,NVOL=3)
37  #include "gppdigi.inc"  #include "gppdigi.inc"
38  c +SEQ,GPDSPE.  c +SEQ,GPDSPE.
# Line 39  c +SEQ,GPDSPE. Line 49  c +SEQ,GPDSPE.
49        INTEGER NTRA,IDIG,LTRA,NPLANE,LKDN,NTDIGI,JJ        INTEGER NTRA,IDIG,LTRA,NPLANE,LKDN,NTDIGI,JJ
50        INTEGER KDIGI(4),NUMBVD(NVOL)        INTEGER KDIGI(4),NUMBVD(NVOL)
51        REAL XI,YI,XO,YO,PATH,DE,GPGAUS        REAL XI,YI,XO,YO,PATH,DE,GPGAUS
52          REAL VPOS(3),BMAGNET(3),ANGLEX,ANGLEY,XMULT,YMULT,DXPOS,DYPOS
53  *  *
54        DATA NTRA,LTRA/1,1/        DATA NTRA,LTRA/1,1/
55  *  *
# Line 138  c +SEQ,GPDSPE. Line 149  c +SEQ,GPDSPE.
149           PATHSPE(I) = VHITS(8,I)           PATHSPE(I) = VHITS(8,I)
150           IPARSPE(I) = VHITS(9,I)           IPARSPE(I) = VHITS(9,I)
151           P0SPE (I) = VHITS(10,I)           P0SPE (I) = VHITS(10,I)
152    C
153    C ADD RESOLUTION SIMULATION TO XAVSPE
154    C        
155             ANGLEX=0.
156             ANGLEY=0.
157             IF(ZOUTSPE(I).NE.ZINSPE(I)) THEN
158                ANGLEX=RADDEG*ATAN( ABS(XOUTSPE(I)-XINSPE(I) )/
159         +              ABS(ZOUTSPE(I)-ZINSPE(I) ) )
160                ANGLEY=RADDEG*ATAN( ABS(YOUTSPE(I)-YINSPE(I) )/
161         +              ABS(ZOUTSPE(I)-ZINSPE(I) ) )          
162             ENDIF
163             VPOS(1)=XAVSPE(I)
164             VPOS(2)=YAVSPE(I)
165             VPOS(3)=ZAVSPE(I)      
166             CALL GUFLD(VPOS,BMAGNET)
167    
168             CALL GPRSPE(ANGLEX,ANGLEY,BMAGNET,XMULT,YMULT,DXPOS,DYPOS)
169    
170             NXMULT(I)=INT(XMULT)
171             NYMULT(I)=INT(YMULT)
172    
173             XAVSPE(I) = XAVSPE(I) + DXPOS
174             YAVSPE(I) = YAVSPE(I) + DYPOS
175    
176    C
177    C ADD THE SHIFT INDUCED BY THE MAG. FIELD
178    C
179    
180             IF(ITRPB(I).EQ.1) THEN
181                XAVSPE(I) = XAVSPE(I) - 2.64e-4/4.5*BMAGNET(2)
182             ELSE IF(ITRPB(I).EQ.2) THEN
183                XAVSPE(I) = XAVSPE(I) + 2.64e-4/4.5*BMAGNET(2)
184             ENDIF
185    
186        ENDDO        ENDDO
187  *  *
188  * Loop on the hits retrived to fill the CWN common and the DIGI bank  * Loop on the hits retrived to fill the CWN common and the DIGI b
189  *  *
190  C fc       IF (HF) NTHSPE = NHITS  C fc       IF (HF) NTHSPE = NHITS
191  C fc       DO I=1,NHITS  C fc       DO I=1,NHITS
Legend:
Removed from v.3.2  
changed lines
  Added in v.3.5
  ViewVC Help
Powered by ViewVC 1.1.23