/[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.9 by bottai, Tue May 2 12:00:24 2006 UTC revision 3.11 by pam-fi, Fri Apr 23 10:03:06 2010 UTC
# Line 1  Line 1 
1  *  *
2  * $Id: gpdspe.F,v 3.8 2006/04/09 23:29:19 cafagna Exp $  * $Id: gpdspe.F,v 3.10 2007/05/08 15:05:18 bottai Exp $
3  *  *
4  * $Log: gpdspe.F,v $  * $Log: gpdspe.F,v $
5    * Revision 3.10  2007/05/08 15:05:18  bottai
6    *  avoid simulating Lorentz angle drift in SPE spread positions
7    *  since new reconstruction correct it at strip level
8    *
9    * Revision 3.9  2006/05/02 12:00:24  bottai
10    *  call gpucollection for strip simulation, condition on FFIELD repositioned
11    *
12  * Revision 3.8  2006/04/09 23:29:19  cafagna  * Revision 3.8  2006/04/09 23:29:19  cafagna
13  * Several new things, among this: ND and CARD  * Several new things, among this: ND and CARD
14  *  *
# Line 184  c +SEQ,GPDSPE. Line 191  c +SEQ,GPDSPE.
191  *  *
192        NTHSPE = NHITS        NTHSPE = NHITS
193        DO I=1,NHITS        DO I=1,NHITS
194           ITRPB(I)   = NUMBV(1,I)           NSPEPLANE=INT((NUMBV(2,I)-1)/6)+1
195           ITRSL(I)   = NUMBV(2,I)           NUMSENS=NUMBV(2,I)-(NSPEPLANE-1)*6  
196             ITRPB(I)   = NSPEPLANE
197             ITRSL(I)   = NUMSENS
198           ITSPA(I)   = NUMBV(3,I)           ITSPA(I)   = NUMBV(3,I)
199  #if defined(GPAMELA_INTER)  #if defined(GPAMELA_INTER)
200           XAVSPE(I) = VHITS(1,I)           XAVSPE(I) = VHITS(1,I)
# Line 242  C       Line 251  C      
251  C  C
252  C ADD THE SHIFT INDUCED BY THE MAG. FIELD  C ADD THE SHIFT INDUCED BY THE MAG. FIELD
253  C  C
254           IF(FFIELD.NE.0) THEN  cv
255            cv 8/05/2007 not anymore needed since the effect is corrected in reco
256              IF(ITRPB(I).LT.6) THEN  cv before calculating the positions for fit
257                 XAVSPE(I) = XAVSPE(I) - 2.64e-4/4.5*BMAGNET(2)  cv
258              ELSE IF(ITRPB(I).EQ.6) THEN  cv         IF(FFIELD.NE.0) THEN
259                 XAVSPE(I) = XAVSPE(I) + 2.64e-4/4.5*BMAGNET(2)  cv      
260              ENDIF  cv            IF(ITRPB(I).LT.6) THEN
261           ENDIF    cv               XAVSPE(I) = XAVSPE(I) + 2.64e-4/4.5*BMAGNET(2)
262    cv            ELSE IF(ITRPB(I).EQ.6) THEN
263    cv               XAVSPE(I) = XAVSPE(I) - 2.64e-4/4.5*BMAGNET(2)
264    cv            ENDIF
265    cv         ENDIF  
266        ENDDO        ENDDO
267    
268  *  *
Legend:
Removed from v.3.9  
changed lines
  Added in v.3.11
  ViewVC Help
Powered by ViewVC 1.1.23