]>
Commit | Line | Data |
---|---|---|
fe4da5cc | 1 | * |
2 | * $Id$ | |
3 | * | |
4 | * $Log$ | |
5 | * Revision 1.1.1.1 1995/10/24 10:19:42 cernlib | |
6 | * Geant | |
7 | * | |
8 | * | |
9 | #include "geant321/pilot.h" | |
10 | *CMZ : 3.21/02 29/03/94 15.41.31 by S.Giani | |
11 | *-- Author : | |
12 | SUBROUTINE CGBRIK(DX,DY,DZ,LCG,CG) | |
13 | ************************************************************************ | |
14 | * * | |
15 | * Name: CGBRIK * | |
16 | * Author: E. Chernyaev Date: 24.03.88 * | |
17 | * Revised: * | |
18 | * * | |
19 | * Function: Create CG-object for BRICK * | |
20 | * * | |
21 | * References: CGSIZE * | |
22 | * * | |
23 | * Input: DX - X-size * | |
24 | * DY - Y-size * | |
25 | * DZ - Z-size * | |
26 | * LCG - max-size of CG-object * | |
27 | * * | |
28 | * Output: CG - CG-object * | |
29 | * CG(1) - length of CG-object * | |
30 | * = 0 if error in parameters * | |
31 | * < 0 if no space * | |
32 | * * | |
33 | * Errors: none * | |
34 | * * | |
35 | ************************************************************************ | |
36 | #include "geant321/cggpar.inc" | |
37 | REAL CG(*),XYZ(3,8) | |
38 | INTEGER IORDER(4,6) | |
39 | SAVE XYZ,IORDER | |
40 | DATA XYZ/24*0./ | |
41 | DATA IORDER/1,2,3,4, | |
42 | + 1,5,6,2, 2,6,7,3, 3,7,8,4, 4,8,5,1, | |
43 | + 8,7,6,5/ | |
44 | *- | |
45 | CG(KCGSIZ) = 0. | |
46 | * T E S T P A R A M E T E R S C O R R E C T N E S S | |
47 | IF (DX.LE.0. .OR. DY.LE.0. .OR. DZ.LE.0.) GOTO 999 | |
48 | * C O M P U T E S I Z E O F C G - O B J E C T | |
49 | NFATYP = 1 | |
50 | NFAC = 6 | |
51 | NEDG = 4 | |
52 | CG(KCGSIZ) = CGSIZE(LCG,NFATYP,NFAC,NEDG) | |
53 | IF (CG(KCGSIZ) .LE. 0.) GOTO 999 | |
54 | * P R E P A R E N O D E S C O O R D I N A T E S | |
55 | XYZ(2,2) = DY | |
56 | XYZ(1,3) = DX | |
57 | XYZ(2,3) = DY | |
58 | XYZ(1,4) = DX | |
59 | DO 100 I=1,4 | |
60 | XYZ(1,I+4) = XYZ(1,I) | |
61 | XYZ(2,I+4) = XYZ(2,I) | |
62 | XYZ(3,I+4) = DZ | |
63 | 100 CONTINUE | |
64 | * C R E A T E C G - O B J E C T F O R B R I K | |
65 | CG(KCGATT) = 0. | |
66 | CG(KCGNF) = 6. | |
67 | JCG = LCGHEA | |
68 | DO 300 NF=1,6 | |
69 | CG(JCG+KCGAF) = 0. | |
70 | X1 = XYZ(1,IORDER(2,NF))-XYZ(1,IORDER(1,NF)) | |
71 | Y1 = XYZ(2,IORDER(2,NF))-XYZ(2,IORDER(1,NF)) | |
72 | Z1 = XYZ(3,IORDER(2,NF))-XYZ(3,IORDER(1,NF)) | |
73 | X2 = XYZ(1,IORDER(3,NF))-XYZ(1,IORDER(2,NF)) | |
74 | Y2 = XYZ(2,IORDER(3,NF))-XYZ(2,IORDER(2,NF)) | |
75 | Z2 = XYZ(3,IORDER(3,NF))-XYZ(3,IORDER(2,NF)) | |
76 | A = Y1*Z2 - Y2*Z1 | |
77 | IF (A .GT. 0.) A = 1. | |
78 | IF (A .LT. 0.) A =-1. | |
79 | B = X2*Z1 - X1*Z2 | |
80 | IF (B .GT. 0.) B = 1. | |
81 | IF (B .LT. 0.) B =-1. | |
82 | C = X1*Y2 - X2*Y1 | |
83 | IF (C .GT. 0.) C = 1. | |
84 | IF (C .LT. 0.) C =-1. | |
85 | D =-(A*XYZ(1,IORDER(1,NF)) + | |
86 | + B*XYZ(2,IORDER(1,NF)) + | |
87 | + C*XYZ(3,IORDER(1,NF)) ) | |
88 | CG(JCG+KCGAA) = A | |
89 | CG(JCG+KCGBB) = B | |
90 | CG(JCG+KCGCC) = C | |
91 | CG(JCG+KCGDD) = D | |
92 | CG(JCG+KCGNE) = 4. | |
93 | JCG = JCG + LCGFAC | |
94 | DO 200 NE=1,4 | |
95 | NE2 = NE + 1 | |
96 | IF (NE .EQ. 4) NE2 = 1 | |
97 | CG(JCG+KCGAE) = 0. | |
98 | CG(JCG+KCGX1) = XYZ(1,IORDER(NE,NF)) | |
99 | CG(JCG+KCGY1) = XYZ(2,IORDER(NE,NF)) | |
100 | CG(JCG+KCGZ1) = XYZ(3,IORDER(NE,NF)) | |
101 | CG(JCG+KCGX2) = XYZ(1,IORDER(NE2,NF)) | |
102 | CG(JCG+KCGY2) = XYZ(2,IORDER(NE2,NF)) | |
103 | CG(JCG+KCGZ2) = XYZ(3,IORDER(NE2,NF)) | |
104 | JCG = JCG + LCGEDG | |
105 | 200 CONTINUE | |
106 | 300 CONTINUE | |
107 | * | |
108 | 999 RETURN | |
109 | END |