GEANT321/geocad/gtuset.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:47  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 GTUSET(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   Converts GEANT shapes TUBE&TUBS to SET
  21 C-
  22 C-   Describes a contour face of a TUBE(Segment) from the GEANT
  23 C-   TUBE/TUBS shape parameters.  This face is rotated around the
  24 C-   Z-axis to generated the solid shape in SET.
  25 C
  26 C---------------------------------------------------------------------
  27 C
  28 #include "geant321/gcsetf.inc"
  29
  30 C
  31       DIMENSION PAR(100)
  32       REAL DX1,DZ,RMIN,RMAX,PHIMIN,PHIMAX,Z2
  33 C
  34 C----------------------------------------------------------------------
  35 C
  36       RMIN=PAR(1)
  37       RMAX=PAR(2)
  38       Z2=PAR(3)
  39       IF (ISHAPE.EQ.5) THEN
  40          PHIMIN=0.0
  41          PHIMAX=360.0
  42       ELSEIF (ISHAPE.EQ.6) THEN
  43          PHIMIN=PAR(4)
  44          PHIMAX=PAR(5)
  45       ENDIF
  46 C       SET CONVERSION
  47       Z= 2*Z2
  48       DX1=0.0
  49       DY=0.0
  50       DZ=Z2
  51       IF (RMAX .LE. 0.0) THEN
  52          WRITE (*,*) 'IMPOSSIBLE RADIUS VALUE'
  53          N1=N1-1
  54          GOTO 10
  55       ENDIF
  56 C *WRITE SET @50,N1,:5,2#32,RMAX,Z,PHIMIN,PHIMAX,RMIN
  57 C *WRITE SET @302,N2,:5,2#317,0,0,-Z2
  58 C *WRITE SET @100,N3,:5,2,:9,'MATNAM'#101,!N1,!N2
  59 C-------------------------------------------------------------------
  60 C
  61       WRITE(BLKSTR,10000)N1,RMAX,Z,PHIMIN,PHIMAX,RMIN
  62       CALL GJWRIT
  63       N1=N1+1
  64       WRITE(BLKSTR,10100)N1,-DX1,-DY,-DZ
  65       CALL GJWRIT
  66       N1=N1+1
  67       WRITE(BLKSTR,10200)N1,N1-2,N1-1
  68       CALL GJWRIT
  69 C
  70 10000   FORMAT('@50,',I10,',:5,2#32,',G14.7,',',G14.7,','
  71      +  ,G14.7,',',G14.7,',',G14.7)
  72 10100   FORMAT('@302,',I10,'#317,',G14.7,',',G14.7,',',G14.7)
  73 10200   FORMAT('@100,',I10,',:5,2#101,!',I10,',!',I10)
  74 C
  75    10 RETURN
  76       END