This commit was generated by cvs2svn to compensate for changes in r2,

[u/mrichter/AliRoot.git] / GEANT321 / ggeom / gnotr1.F

diff --git a/GEANT321/ggeom/gnotr1.F b/GEANT321/ggeom/gnotr1.F
new file mode 100644 (file)
index 0000000..462b1bd
--- /dev/null
+++ b/GEANT321/ggeom/gnotr1.F
@@ -0,0 +1,78 @@
+*
+* $Id$
+*
+* $Log$
+* Revision 1.1.1.1  1995/10/24 10:20:53  cernlib
+* Geant
+*
+*
+#include "geant321/pilot.h"
+*CMZ :  3.21/02 29/03/94  15.41.29  by  S.Giani
+*-- Author :
+      SUBROUTINE GNOTR1(P)
+C.
+C.    ******************************************************************
+C.    *                                                                *
+C.    *    Routine to extend  the  user supplied parameters P  to      *
+C.    *    internal parameters   for the  trapezoidal shape TRAP       *
+C.    *                                                                *
+C.    *         Called by : GSPOSP, GSVOLU , GGDPAR                    *
+C.    *         Author  R.Nierhaus  *********                          *
+C.    *                                                                *
+C.    ******************************************************************
+C.
+*     GNOTR1  computes the coordinates for the 8 vertices
+*     X(8), Y(8) and Z(8).
+*     Then GNOTR2 is called for each of the six boundary surfaces.
+*
+      DIMENSION P(35),X(8),Y(8),Z(8)
+*
+      HZ  =P(1)
+      SX  =P(2)
+      SY  =P(3)
+      HY1 =P(4)
+      HX11=P(5)
+      HX12=P(6)
+      SX1 =P(7)
+      HY2 =P(8)
+      HX21=P(9)
+      HX22=P(10)
+      SX2 =P(11)
+      IF (HY1.LT.1.E-6) HY1=1.E-6
+      IF (HX11.LT.1.E-6) HX11=1.E-6
+      IF (HX12.LT.1.E-6) HX12=1.E-6
+      IF (HY2.LT.1.E-6) HY2=1.E-6
+      IF (HX21.LT.1.E-6) HX21=1.E-6
+      IF (HX22.LT.1.E-6) HX22=1.E-6
+      X(1)=-HZ*SX-HY1*SX1-HX11
+      X(2)=-HZ*SX-HY1*SX1+HX11
+      X(3)=-HZ*SX+HY1*SX1-HX12
+      X(4)=-HZ*SX+HY1*SX1+HX12
+      X(5)=+HZ*SX-HY2*SX2-HX21
+      X(6)=+HZ*SX-HY2*SX2+HX21
+      X(7)=+HZ*SX+HY2*SX2-HX22
+      X(8)=+HZ*SX+HY2*SX2+HX22
+      Y(1)=-HZ*SY-HY1
+      Y(2)=-HZ*SY-HY1
+      Y(3)=-HZ*SY+HY1
+      Y(4)=-HZ*SY+HY1
+      Y(5)=+HZ*SY-HY2
+      Y(6)=+HZ*SY-HY2
+      Y(7)=+HZ*SY+HY2
+      Y(8)=+HZ*SY+HY2
+      Z(1)=-HZ
+      Z(2)=-HZ
+      Z(3)=-HZ
+      Z(4)=-HZ
+      Z(5)=+HZ
+      Z(6)=+HZ
+      Z(7)=+HZ
+      Z(8)=+HZ
+      CALL GNOTR2(X,Y,Z,1,2,4,3,P(12))
+      CALL GNOTR2(X,Y,Z,5,7,8,6,P(16))
+      CALL GNOTR2(X,Y,Z,1,5,6,2,P(20))
+      CALL GNOTR2(X,Y,Z,3,4,8,7,P(24))
+      CALL GNOTR2(X,Y,Z,1,3,7,5,P(28))
+      CALL GNOTR2(X,Y,Z,2,6,8,4,P(32))
+*
+      END