[u/mrichter/AliRoot.git] / GEANT321 / gdraw / gdsarc.F

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:20:27  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.27  by  S.Giani
*-- Author :
      SUBROUTINE GDSARC(R,PHI1,PHI2,X,NP)
C.
C.    ******************************************************************
C.    *                                                                *
C.    *       Computes and store in X(1..2,1..NP) arc points           *
C.    *       (only X and Y coordinates) for GDRWSC                    *
C.    *       of the circle of radius R centered on Z-axis,            *
C.    *       from PHI1 to PHI2,                                       *
C.    *       on a plane parallel to X-Y plane                         *
C.    *                                                                *
C.    *    ==>Called by : GDRWSC                                       *
C.    *       Author : P.Zanarini   *********                          *
C.    *                                                                *
C.    ******************************************************************
C.
#include "geant321/gconsp.inc"
#include "geant321/gcdraw.inc"
#include "geant321/gcunit.inc"
      DIMENSION X(2,42)
      SAVE MAXNP
      DATA MAXNP/42/
C.
C.    ------------------------------------------------------------------
C.
      PHIMIN=PHI1
      PHIMAX=PHI2
C
      NP=1
      PH=PHIMIN*DEGRAD
      X(1,1)=R*COS(PH)
      X(2,1)=R*SIN(PH)
C
      IF(ABS(PHIMAX-PHIMIN).LE.9.) GO TO 20
C
      I1=PHIMIN/9.0+2
      I2=PHIMAX/9.0+1
      NAS=I2-I1+1
      IF(NAS.LE.0) NAS=NAS+40
      IC=I1
C
      IF (NAS.GT.MAXNP-2) GO TO 30
C
      DO 10 IAS=1,NAS
         IF(IC.GT.40) IC=IC-40
         NP=NP+1
         ICC=IC
         IF (ICC.LE.0) ICC=40+ICC
         X(1,NP)=R*GCOS(ICC)
         X(2,NP)=R*GSIN(ICC)
         IC=IC+1
   10 CONTINUE
C
   20 CONTINUE
C
      PH=PHIMAX*DEGRAD
      NP=NP+1
      X(1,NP)=R*COS(PH)
      X(2,NP)=R*SIN(PH)
      GO TO 999
C
   30 CONTINUE
      NP=NAS+2
      WRITE (CHMAIL,1000) NP
      CALL GMAIL(0,0)
C
 1000 FORMAT (' GDSARC: INTERNAL ERROR !!!   -   NP =',I5)
  999 RETURN
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1995/10/24 10:20:27 cernlib
	6	* Geant
	7	*
	8	*
	9	#include "geant321/pilot.h"
	10	*CMZ : 3.21/02 29/03/94 15.41.27 by S.Giani
	11	*-- Author :
	12	SUBROUTINE GDSARC(R,PHI1,PHI2,X,NP)
	13	C.
	14	C. ******************************************************************
	15	C. * *
	16	C. * Computes and store in X(1..2,1..NP) arc points *
	17	C. * (only X and Y coordinates) for GDRWSC *
	18	C. * of the circle of radius R centered on Z-axis, *
	19	C. * from PHI1 to PHI2, *
	20	C. * on a plane parallel to X-Y plane *
	21	C. * *
	22	C. * ==>Called by : GDRWSC *
	23	C. * Author : P.Zanarini ********* *
	24	C. * *
	25	C. ******************************************************************
	26	C.
	27	#include "geant321/gconsp.inc"
	28	#include "geant321/gcdraw.inc"
	29	#include "geant321/gcunit.inc"
	30	DIMENSION X(2,42)
	31	SAVE MAXNP
	32	DATA MAXNP/42/
	33	C.
	34	C. ------------------------------------------------------------------
	35	C.
	36	PHIMIN=PHI1
	37	PHIMAX=PHI2
	38	C
	39	NP=1
	40	PH=PHIMIN*DEGRAD
	41	X(1,1)=R*COS(PH)
	42	X(2,1)=R*SIN(PH)
	43	C
	44	IF(ABS(PHIMAX-PHIMIN).LE.9.) GO TO 20
	45	C
	46	I1=PHIMIN/9.0+2
	47	I2=PHIMAX/9.0+1
	48	NAS=I2-I1+1
	49	IF(NAS.LE.0) NAS=NAS+40
	50	IC=I1
	51	C
	52	IF (NAS.GT.MAXNP-2) GO TO 30
	53	C
	54	DO 10 IAS=1,NAS
	55	IF(IC.GT.40) IC=IC-40
	56	NP=NP+1
	57	ICC=IC
	58	IF (ICC.LE.0) ICC=40+ICC
	59	X(1,NP)=R*GCOS(ICC)
	60	X(2,NP)=R*GSIN(ICC)
	61	IC=IC+1
	62	10 CONTINUE
	63	C
	64	20 CONTINUE
65	C
66	PH=PHIMAX*DEGRAD
67	NP=NP+1
68	X(1,NP)=R*COS(PH)
69	X(2,NP)=R*SIN(PH)
70	GO TO 999
71	C
72	30 CONTINUE
73	NP=NAS+2
74	WRITE (CHMAIL,1000) NP
75	CALL GMAIL(0,0)
76	C
77	1000 FORMAT (' GDSARC: INTERNAL ERROR !!! - NP =',I5)
78	999 RETURN
79	END