/[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.1.1.1 by cafagna, Thu Jul 11 16:02:13 2002 UTC revision 3.7 by bottai, Tue Feb 7 11:38:52 2006 UTC
# Line 1  Line 1 
1  *  *
2  * $Id$  * $Id: gpdspe.F,v 3.6 2006/02/06 14:17:21 bottai Exp $
3    *
4    * $Log: gpdspe.F,v $
5    * Revision 3.6  2006/02/06 14:17:21  bottai
6    *  adapt again the resolution simulation to the ITRPB volume identification
7    *
8    * Revision 3.5  2006/01/24 13:52:01  bottai
9    * prepare to store the x and y multiplicity of the hit
10    *
11    * Revision 3.4  2006/01/18 09:41:33  bottai
12    *  Some correction to the spe plane ID due to the new geometry tree
13    *
14    * Revision 3.3  2005/12/05 12:15:21  pam-ba
15    * new spectrometer geometry and internal magnetic field
16    *
17    * Revision 3.2  2005/07/25 11:53:21  cafagna
18    * Several updates. See history for details
19    *
20    * Revision 3.1.1.1  2002/07/11 16:02:13  cafagna
21    * First GPAMELA release on CVS
22  *  *
 * $Log$  
23  *  *
24  *CMZ :  2.01/00 05/04/2000  14.35.18  by  Marialuigia Ambriola  *CMZ :  2.01/00 05/04/2000  14.35.18  by  Marialuigia Ambriola
25  *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 10  Line 28 
28  *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
29  *-- Author :    Francesco Cafagna   23/02/96  *-- Author :    Francesco Cafagna   23/02/96
30        SUBROUTINE GPDSPE        SUBROUTINE GPDSPE
31  ************************************************************************  *************************************************************************
32  *                                                                      *  *                                                                       *
33  * To perform SPE digitization                                          *  * To perform SPE digitization                                           *
34  *                                                                      *  *                                                                       *
35  * Called by: GUDIGI                                                    *  * Called by: GUDIGI                                                     *
36  * Author: Francesco Cafagna, 23/02/96 12.50.38                         *  * Author: Francesco Cafagna, 23/02/96 12.50.38                          *
37  *                                                                      *  *                                                                       *
38  ************************************************************************  *************************************************************************
39  #include "gcunit.inc"  #include "gcunit.inc"
40  #include "gpphit.inc"  #include "gpphit.inc"
41        PARAMETER (NH=NHSPE,NHM=NHMSPE,NVOL=4)        PARAMETER (RADDEG=57.295779)
42          PARAMETER (NH=NHSPE,NHM=NHMSPE,NVOL=3)
43  #include "gppdigi.inc"  #include "gppdigi.inc"
44  c +SEQ,GPDSPE.  c +SEQ,GPDSPE.
45  #include "gprhit.inc"  #include "gprhit.inc"
# Line 36  c +SEQ,GPDSPE. Line 55  c +SEQ,GPDSPE.
55        INTEGER NTRA,IDIG,LTRA,NPLANE,LKDN,NTDIGI,JJ        INTEGER NTRA,IDIG,LTRA,NPLANE,LKDN,NTDIGI,JJ
56        INTEGER KDIGI(4),NUMBVD(NVOL)        INTEGER KDIGI(4),NUMBVD(NVOL)
57        REAL XI,YI,XO,YO,PATH,DE,GPGAUS        REAL XI,YI,XO,YO,PATH,DE,GPGAUS
58          REAL VPOS(3),BMAGNET(3),ANGLEX,ANGLEY,XMULT,YMULT,DXPOS,DYPOS
59  *  *
60        DATA NTRA,LTRA/1,1/        DATA NTRA,LTRA/1,1/
61  *  *
# Line 107  c +SEQ,GPDSPE. Line 127  c +SEQ,GPDSPE.
127        NTHSPE = NHITS        NTHSPE = NHITS
128        DO I=1,NHITS        DO I=1,NHITS
129           ITRPB(I)   = NUMBV(1,I)           ITRPB(I)   = NUMBV(1,I)
130           ITRSL(I)   = NUMBV(3,I)           ITRSL(I)   = NUMBV(2,I)
131           ITSPA(I)   = NUMBV(4,I)           ITSPA(I)   = NUMBV(3,I)
132  #if defined(GPAMELA_INTER)  #if defined(GPAMELA_INTER)
133           XAVSPE(I) = VHITS(1,I)           XAVSPE(I) = VHITS(1,I)
134           YAVSPE(I) = VHITS(2,I)           YAVSPE(I) = VHITS(2,I)
# Line 135  c +SEQ,GPDSPE. Line 155  c +SEQ,GPDSPE.
155           PATHSPE(I) = VHITS(8,I)           PATHSPE(I) = VHITS(8,I)
156           IPARSPE(I) = VHITS(9,I)           IPARSPE(I) = VHITS(9,I)
157           P0SPE (I) = VHITS(10,I)           P0SPE (I) = VHITS(10,I)
158    C
159    C ADD RESOLUTION SIMULATION TO XAVSPE
160    C        
161             ANGLEX=0.
162             ANGLEY=0.
163             IF(ZOUTSPE(I).NE.ZINSPE(I)) THEN
164                ANGLEX=RADDEG*ATAN( ABS(XOUTSPE(I)-XINSPE(I) )/
165         +              ABS(ZOUTSPE(I)-ZINSPE(I) ) )
166                ANGLEY=RADDEG*ATAN( ABS(YOUTSPE(I)-YINSPE(I) )/
167         +              ABS(ZOUTSPE(I)-ZINSPE(I) ) )          
168             ENDIF
169             VPOS(1)=XAVSPE(I)
170             VPOS(2)=YAVSPE(I)
171             VPOS(3)=ZAVSPE(I)      
172             CALL GUFLD(VPOS,BMAGNET)
173    
174             CALL GPRSPE(ANGLEX,ANGLEY,BMAGNET,XMULT,YMULT,DXPOS,DYPOS)
175    
176             NXMULT(I)=INT(XMULT)
177             NYMULT(I)=INT(YMULT)
178    
179             XAVSPE(I) = XAVSPE(I) + DXPOS
180             YAVSPE(I) = YAVSPE(I) + DYPOS
181    
182    C
183    C ADD THE SHIFT INDUCED BY THE MAG. FIELD
184    C
185    
186             IF(ITRPB(I).LT.6) THEN
187                XAVSPE(I) = XAVSPE(I) - 2.64e-4/4.5*BMAGNET(2)
188             ELSE IF(ITRPB(I).EQ.6) THEN
189                XAVSPE(I) = XAVSPE(I) + 2.64e-4/4.5*BMAGNET(2)
190             ENDIF
191    
192        ENDDO        ENDDO
193  *  *
194  * 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
195  *  *
196  C fc       IF (HF) NTHSPE = NHITS  C fc       IF (HF) NTHSPE = NHITS
197  C fc       DO I=1,NHITS  C fc       DO I=1,NHITS
Legend:
Removed from v.3.1.1.1  
changed lines
  Added in v.3.7
  ViewVC Help
Powered by ViewVC 1.1.23