1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
|
*DECK MCGPT1
SUBROUTINE MCGPT1(N2SOU,N2REG,NZP,NR2D,INDREG,Z,NOM2D,T2D,I,K,ZI,
1 TK,TIN,CPOI,SPOI,TPO,TPOI,N3D,NOM3D,H3D)
*
*-----------------------------------------------------------------------
*
*Purpose:
* Reconstruct 3D track for a 3D prismatic geometry from a 2D track.
* polar angle in [0, pi/2] case.
*
*Copyright:
* Copyright (C) 2002 Ecole Polytechnique de Montreal
* This library is free software; you can redistribute it and/or
* modify it under the terms of the GNU Lesser General Public
* License as published by the Free Software Foundation; either
* version 2.1 of the License, or (at your option) any later version
*
*Author(s): R. Le Tellier
*
*Parameters: input
* N2SOU number of external surfaces in the 2D tracking.
* N2REG number of regions in the 2D tracking.
* NZP number of z-planes.
* NR2D number of segments corresponding to regions for this 2D track.
* INDREG region/surface index to go from the 2D to the 3D geometry.
* Z z-plan coordinates.
* NOM2D vector containing the region number of the different segments
* of this 2D track.
* T2D vector containing the local coordinates of the segments
* boundaries for this 2D track.
* CPOI inverse of the polar cosine.
* SPOI inverse of the polar sine.
* TPO polar tangent.
* TPOI polar cotangent.
*
*Parameters: input/output
* I starting/ending z plan.
* K starting/ending x-y tracking segment.
* ZI strating/ending z coordinate.
* TK starting/ending x-y tracking coordinate.
* TIN orientation of the starting/surface.
*
*Parameters: output
* N3D number of segments for this 3D track.
* NOM3D vector containing the region number of the different segments
* of this 3D track.
* H3D vector containing the length of the different segments of this
* 3D track.
*
*-----------------------------------------------------------------------
*
IMPLICIT NONE
*----
* SUBROUTINE ARGUMENTS
*----
INTEGER N2SOU,N2REG,NZP,NR2D,INDREG(-N2SOU:N2REG,0:NZP+1),
1 NOM2D(NR2D+2),I,K,N3D,NOM3D(*),TIN
REAL Z(0:NZP)
DOUBLE PRECISION ZI,TK,CPOI,SPOI,TPO,TPOI,H3D(N3D),T2D(0:NR2D)
*---
* LOCAL VARIABLES
*---
!! INTEGER II,N3DP
DOUBLE PRECISION DELZ,DELT,X,L
!! REAL Tstart,Zstart
*
!!!! Tstart=TK
!!!! Zstart=ZI
!!!! N3DP=N3D
DO WHILE ((I.LE.NZP).AND.(K.LE.NR2D))
N3D=N3D+1
NOM3D(N3D)=INDREG(NOM2D(K+1),I)
IF (TIN.EQ.0) THEN
* track enters the region through bottom boundary
DELZ=Z(I)-Z(I-1)
DELT=T2D(K)-TK
ELSE
* track enters the region through left boundary
DELZ=Z(I)-ZI
DELT=T2D(K)-T2D(K-1)
ENDIF
X=DELT/DELZ
IF (TPO.LT.X) THEN
* track leaves the region through the top boundary
ZI=Z(I)
I=I+1
TK=TK+DELZ*TPO
L=DELZ*CPOI
TIN=0
ELSE
* track leaves the region through the right boundary
TK=T2D(K)
K=K+1
ZI=ZI+DELT*TPOI
L=DELT*SPOI
TIN=1
ENDIF
IF (L.GT.0.0) THEN
H3D(N3D)=DBLE(L)
ELSE
N3D=N3D-1
ENDIF
ENDDO
N3D=N3D+1
H3D(N3D)=0.5
NOM3D(N3D)=INDREG(NOM2D(K+1),I)
!!!! call MCGPTP(8,N3D-N3DP+1,NOM3D(N3DP),Tstart,ZStart,CPO,SPO,
!!!! 1 H3D(N3DP),1)
*
RETURN
*
END
!!!!*DECK MCGPTP
!!!! SUBROUTINE MCGPTP(IPOUT,N3D,NOM3D,TS,ZS,CPO,SPO,H3D,DIR)
!!!!*
!!!! IMPLICIT NONE
!!!!*
!!!! INTEGER IPOUT,N3D,NOM3D(N3D),DIR
!!!! REAL TS,ZS,H3D(N3D)
!!!! DOUBLE PRECISION CPO,SPO
!!!!*
!!!! INTEGER II
!!!! REAL ZE,TE
!!!!*
!!!! WRITE(IPOUT,201) TS,ZS,NOM3D(1)
!!!! DO II=2,N3D-1
!!!! TE=TS+H3D(II)*SPO
!!!! IF (DIR.EQ.1) ZE=ZS+H3D(II)*CPO
!!!! IF (DIR.EQ.2) ZE=ZS-H3D(II)*CPO
!!!! WRITE(IPOUT,100) TS,TE,ZS,ZE
!!!! WRITE(IPOUT,200) 0.5*(TS+TE),0.5*(ZS+ZE),NOM3D(II)
!!!! TS=TE
!!!! ZS=ZE
!!!! ENDDO
!!!! WRITE(IPOUT,202) TS,ZS,NOM3D(N3D)
!!!!*
!!!! RETURN
!!!!*
!!!! 100 FORMAT(
!!!! 1 7H line([,E16.8,1H,,E16.8,3H],[,E16.8,1H,,E16.8,2H],,
!!!! 2 14H 'Marker','+')/)
!!!! 200 FORMAT(
!!!! 1 6H text(,E16.8,1H,,E16.8,1H,,1H',I4,2H',,
!!!! 2 32H 'HorizontalAlignment','center')/)
!!!! 201 FORMAT(
!!!! 1 6H text(,E16.8,1H,,E16.8,1H,,1H',I4,2H',,
!!!! 2 27H 'VerticalAlignment','top',,
!!!! 3 31H 'HorizontalAlignment','right')/)
!!!! 202 FORMAT(
!!!! 1 6H text(,E16.8,1H,,E16.8,1H,,1H',I4,2H',,
!!!! 2 30H 'VerticalAlignment','bottom',,
!!!! 3 30H 'HorizontalAlignment','left')/)
!!!! END
|
