GEANT321/geocad/gcoset.F

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1999/05/18 15:55:17  fca
   6 * AliRoot sources
   7 *
   8 * Revision 1.1.1.1  1995/10/24 10:20:46  cernlib
   9 * Geant
  10 *
  11 *
  12 #include "geant321/pilot.h"
  13 *CMZ :  3.21/02 29/03/94  15.41.32  by  S.Giani
  14 *-- Author :
  15       SUBROUTINE GCOSET(PAR,ISHAPE)
  16 C-
  17 C-   Created  12-JUN-1991   Nils Joar Hoimyr
  18 C-   Modified 21.02.1992    Jouko Vuoskoski
  19 C-
  20 C-   Describes a contour face of a CONE(Segment) from the GEANT
  21 C-   CONE/CONS shape parameters.  This face is rotated around the
  22 C-   Z-axis to generated the solid shape in SET.
  23 C----------------------------------------------------------C
  24 #include "geant321/gcsetf.inc"
  25
  26 C
  27       DIMENSION PAR(100)
  28       REAL P1X,P1Y,P1Z,P2X,P2Y,P2Z,P3X,P3Y,P3Z,P4X,P4Y,P4Z
  29 C
  30 C---------------------------------------------------------
  31       RMIN1=PAR(2)
  32       RMAX1=PAR(3)
  33       RMIN2=PAR(4)
  34       RMAX2=PAR(5)
  35       Z2=PAR(1)
  36       Z1=-Z2
  37       IF (ISHAPE.EQ.7) THEN
  38          PHIMIN=0.0
  39          PHIMAX=360.0
  40       ELSEIF (ISHAPE.EQ.8) THEN
  41          PHIMIN=PAR(6)
  42          PHIMAX=PAR(7)
  43       ENDIF
  44 C       SET CONVERSION
  45       P1X= RMIN1
  46       P1Y= 0
  47       P1Z= Z1
  48       P2X= RMIN2
  49       P2Y= 0
  50       P2Z= Z2
  51       P3X= RMAX2
  52       P3Y= 0
  53       P3Z= Z2
  54       P4X= RMAX1
  55       P4Y= 0
  56       P4Z= Z1
  57 C *WRITE SET@103,F1#3,3,2,P1X,P1Y,P1Z,P2X,P2Y,P2Z,P3X,P3Y,P3Z,P4X,P4Y,P4Z
  58 C *WRITE SET @100,N..#144,!F1,PHIMIN,PHIMAX.
  59 C----------------------------------------------------------------------
  60 C
  61       WRITE(BLKSTR,10000)N1,P1X,P1Y,P1Z,P2X,P2Y,P2Z, P3X,P3Y,P3Z,P4X,
  62      +P4Y,P4Z
  63       CALL GJWRIT
  64       N1=N1+1
  65       WRITE(BLKSTR,10100)N1,N1-1,PHIMIN,PHIMAX
  66       CALL GJWRIT
  67 C
  68 10000   FORMAT('@103,',I10,',:5,2#3,3,2,',G14.7,',',G14.7,',',G14.7,','
  69      +  ,G14.7,',',G14.7,',',G14.7,',',G14.7,',',G14.7,',',G14.7,','
  70      +  ,G14.7,',',G14.7,',',G14.7)
  71 10100   FORMAT('@100,',I10,',:5,2#144,!',I10,',',G14.7,',',G14.7)
  72 C
  73 C
  74       END