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

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1999/05/18 15:55:03  fca
* AliRoot sources
*
* Revision 1.1.1.1  1995/10/24 10:20:23  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/04 10/01/95  16.43.07  by  S.Giani
*-- Author :
      SUBROUTINE GDPRTR(NAME,YLEV,XLEV,SIZE,FACHV,
     +IDUP,IDRAW,IONL,IDET,ISEEN)
C.
C.    ******************************************************************
C.    *                                                                *
C.    *       Draws a box centered in (XLEV,YLEV) of size SIZE.        *
C.    *       On the top half box draws its name NAME and, if IDUP<>1, *
C.    *       on the bottom half draw IDUP preceded by a symbol        *
C.    *       (IDUP<0 means division, IDUP>1 means multiple position,  *
C.    *       IDUP=1 means single position).                           *
C.    *       If IDRAW=1 then over the box draws a normal-view of      *
C.    *       volume NAME (with same seen attributes set by GDSPEC).   *
C.    *       If IONL<>1 then draw two other boxes, slightly shifted,  *
C.    *       behind the original one.                                 *
C.    *       If IDET=1 then draw another box, slightly larger,        *
C.    *       around the original one.                                 *
C.    *       If ISEEN=0 then draw the box(es) with dashed lines.      *
C.    *                                                                *
C.    *    ==>Called by : GDTREE                                       *
C.    *       Author : P.Zanarini ; S.Giani  *********                 *
C.    *                                                                *
C.    ******************************************************************
C.
#include "geant321/gcbank.inc"
#include "geant321/gcnum.inc"
#include "geant321/gcdraw.inc"
      DIMENSION X(5),Y(5)
      DIMENSION NAMT0(16),NAM0(4)
      DIMENSION PAR(100)
      DIMENSION NNDM(14),INDM(4,14)
      DIMENSION ATT(10)
      CHARACTER*4 NAME,NAMSEE,ISON
      CHARACTER*14 CHNAM0
      SAVE NNDM,INDM
      COMMON/NPILEV/NPILE
C
      DATA NNDM/0,0,0,4,0,2,0,2,4,3,4,3,0,0/
      DATA INDM/0,0,0,0, 0,0,0,0, 0,0,0,0, 2,3,7,11, 0,0,0,0,
     +          4,5,0,0, 0,0,0,0, 6,7,0,0, 3,4,5,6, 4,5,6,0,
     +          1,2,3,4, 1,2,3,0, 0,0,0,0, 0,0,0,0/
C.
C.    ------------------------------------------------------------------
C.
      IF (ISEEN.EQ.0) THEN
         LINSTY=4
      ELSE
         LINSTY=1
      ENDIF
      CALL MVBITS(LINSTY,0,3,LINATT,10)
C
C            Compute box vertex
C
      SIZ2=SIZE/2.
      SIZ4=SIZE/4.
      SIZ8=SIZE/8.
***      X(1)=XLEV-SIZ2
      X(1)=XLEV-(SIZ2*2.*FACHV)
      X(2)=X(1)
***      X(3)=XLEV+SIZ2
      X(3)=XLEV+(SIZ2*2.*FACHV)
      X(4)=X(3)
      X(5)=X(1)
      Y(1)=YLEV+SIZ2
      Y(2)=YLEV-SIZ2
      Y(3)=Y(2)
      Y(4)=Y(1)
      Y(5)=Y(1)
C
C            Draw box
C
      IWTY=IGIWTY(1)
      CALL IGQWK(0,'HWCO',R)
      IF(IWTY.GT.10.OR.IWTY.LT.1.OR.R.EQ.2.)THEN
         LINCOL=1
         CALL MVBITS(LINCOL,0,8,LINATT,16)
         CALL GDRAWV(X,Y,5)
      ELSE
         CALL IGSET('LWID',1.)
         CALL UCTOH(NAME,IIIN,4,4)
         CALL IGPID(1,'Tree',IIIN,' ')
         CALL IGPAVE(X(1),X(3),Y(1),Y(3),SIZ8,1004,1001,'TRS')
      ENDIF
C
C            If volume is also detector draw another box outside
C
      IF (IDET.EQ.1) THEN
         X(1)=X(1)-SIZ8
         X(2)=X(2)-SIZ8
         X(3)=X(3)+SIZ8
         X(4)=X(4)+SIZ8
         X(5)=X(5)-SIZ8
         Y(1)=Y(1)+SIZ8
         Y(2)=Y(2)-SIZ8
         Y(3)=Y(3)-SIZ8
         Y(4)=Y(4)+SIZ8
         Y(5)=Y(5)+SIZ8
         LINCOL=1
         CALL MVBITS(LINCOL,0,8,LINATT,16)
         CALL GDRAWV(X,Y,5)
      ENDIF
C
C            If volume is not 'ONLY' draw two other boxes behind
C
      IF (IONL.NE.1) THEN
         X(1)=XLEV+(SIZ2*2.*FACHV)-SIZ8
         Y(1)=YLEV+SIZ2
         X(2)=X(1)
         Y(2)=Y(1)+SIZ8
         X(3)=X(2)-(SIZE*2.*FACHV)
         Y(3)=Y(2)
         X(4)=X(3)
         Y(4)=Y(3)-SIZE
         X(5)=X(4)+SIZ8
         Y(5)=Y(4)
         LINCOL=1
         CALL MVBITS(LINCOL,0,8,LINATT,16)
         CALL GDRAWV(X,Y,5)
         DO 10 I=1,5
            X(I)=X(I)-SIZ8
            Y(I)=Y(I)+SIZ8
   10    CONTINUE
         CALL GDRAWV(X,Y,5)
      ENDIF
C
      LINSTY=1
      CALL MVBITS(LINSTY,0,3,LINATT,10)
C
C            Draw name
C
      YL=YLEV-SIZ2+(SIZ8/2.)
      IF(FACHV.EQ.1.)YL=YLEV-SIZ8
      XL=XLEV-(SIZ4*FACHV)+SIZ8
      IF(FACHV.EQ.1.)XL=XLEV
      SZ=1.5*SIZ4*2.
      IF(FACHV.EQ.1.)SZ=1.*SIZ4*2.
      ITXCOL=5
      CALL MVBITS(ITXCOL,0,4,ITXATT,3)
      AITXCO=ITXCOL
      CALL IGSET('TXCI',AITXCO)
      CALL IGSET('TXFP',-60.)
      IWTY=IGIWTY(1)
      IF(IWTY.GT.10.OR.IWTY.LT.1)CALL IGSET('TXFP',-61.)
***      CALL IGPID(1,'Tree',NPILE,' ')
      CALL GDRAWT(XL,YL,NAME,SZ,0.,4,0)
      CALL IGSET('TXFP',2.)
      CALL IGSET('TXCI',1.)
C
      IF (IDUP.EQ.1) GO TO 20
C
C            Draw IDUP in bottom half box
C
      CALL UCTOH('<"4#>',NAMT0(1),1,5)
      IF (IDUP.LT.0) CALL UCTOH('F',NAMT0(3),1,1)
      DUP=FLOAT(IABS(IDUP))
      CALL HBCDF(DUP,8,NAMT0(6))
      CALL UCTOH('$',NAMT0(14),1,1)
      CALL UBUNCH(NAMT0,NAM0,14)
      CALL UHTOC(NAM0,4,CHNAM0,14)
      XL=XLEV+(SIZ4*FACHV)+SIZ2
      IF(FACHV.EQ.1.)XL=XLEV-(3.*SIZ8)
      YL=YLEV-SIZ4-SIZ8
      IF(FACHV.EQ.1.)YL=YLEV-SIZ2+(SIZ8/2.)
      SZ=0.7*SIZ4*2.
      IF(FACHV.EQ.1.)SZ=0.7*SIZ4*1.5
      ITXCOL=5
      CALL MVBITS(ITXCOL,0,4,ITXATT,3)
      AITXCO=ITXCOL
      CALL IGSET('TXCI',AITXCO)
      CALL IGSET('TXFP',2.)
***      CALL IGPID(1,'Tree',NPILE,' ')
      CALL GDRAWT(XL,YL,CHNAM0,SZ,0.,10,-1)
      CALL IGSET('TXCI',1.)
C
   20 CONTINUE
C
      IF (IDRAW.NE.1) GO TO 999
C
C             Draw a normal view of the volume
C
      CALL GLOOK(NAME,IQ(JVOLUM+1),NVOLUM,IVO)
      IF(IVO.LE.0) GO TO 999
C
      JVO=LQ(JVOLUM-IVO)
C
C                GET PARAMETERS
C
      CALL GFPARA(NAME,1,0,NPAR,NATT,PAR,ATT)
      IF(NPAR.LE.0) GO TO 999
C
      ISHAPE=Q(JVO+2)
C
C             CHECK parametER SIZES
C
      PARMAX=-1.
      DO 50  I=1,NPAR
         IF(ISHAPE.GT.14)GO TO 50
         IF(NNDM(ISHAPE).LE.0) GO TO 40
C
         NDM=NNDM(ISHAPE)
         DO 30  IDM=1,NDM
C
            IF(I.EQ.INDM(IDM,ISHAPE)) GO TO 50
C
   30    CONTINUE
   40    CONTINUE
C
         ABSPAR=ABS(PAR(I))
         PARMAX=MAX(PARMAX,ABSPAR)
   50 CONTINUE
      XSCALE=SIZE/(3.*PARMAX)
      YSCALE=XSCALE
      XPOS=XLEV-(SIZE*FACHV*2)+SIZ8
      YPOS=YLEV
*      YPOS=YLEV+SIZE
C
      CALL GFATT(NAME,'SEEN',KSEEN)
C
C             Add local value KSEEN 1 to starting node of tree
C
      KSEEN=KSEEN+110
      ISEEN=KSEEN
      CALL GSATT(NAME,'SEEN',ISEEN)
C
      CALL GDNSON(NAME,NSON,IDIV)
      DO 60  N=1,NSON
         CALL GDSON(N,NAME,ISON)
         CALL GFATT(ISON,'SEEN',KSEEN)
C
C             ISON is a volume with multeplicity;
C             first occurrence has already been set
C
         IF (KSEEN.GT.50) GO TO 60
C
C             Add local value SEEN -2 to each one-level-down node
C
         KSEEN=KSEEN+80
         ISEEN=KSEEN
         CALL GSATT(ISON,'SEEN',ISEEN)
   60 CONTINUE
      CALL GDRAW(NAME,45.,135.,0.,XPOS,YPOS,XSCALE,YSCALE)
C
C             Reset global SEEN values
C
      DO 70  IVO=1,NVOLUM
         CALL UHTOC(IQ(JVOLUM+IVO),4,NAMSEE,4)
         CALL GFATT(NAMSEE,'SEEN',KSEEN)
         IF (KSEEN.LT.50) GO TO 70
         ISEENL=KSEEN/10.+0.5
         ISEENG=KSEEN-ISEENL*10.
         CALL GSATT(NAMSEE,'SEEN',ISEENG)
   70 CONTINUE
C
  999 END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
d43b40e2	5	* Revision 1.1.1.1 1999/05/18 15:55:03 fca
	6	* AliRoot sources
	7	*
fe4da5cc	8	* Revision 1.1.1.1 1995/10/24 10:20:23 cernlib
	9	* Geant
	10	*
	11	*
	12	#include "geant321/pilot.h"
	13	*CMZ : 3.21/04 10/01/95 16.43.07 by S.Giani
	14	*-- Author :
	15	SUBROUTINE GDPRTR(NAME,YLEV,XLEV,SIZE,FACHV,
	16	+IDUP,IDRAW,IONL,IDET,ISEEN)
	17	C.
	18	C. ******************************************************************
	19	C. * *
	20	C. * Draws a box centered in (XLEV,YLEV) of size SIZE. *
	21	C. * On the top half box draws its name NAME and, if IDUP<>1, *
	22	C. * on the bottom half draw IDUP preceded by a symbol *
	23	C. * (IDUP<0 means division, IDUP>1 means multiple position, *
	24	C. * IDUP=1 means single position). *
	25	C. * If IDRAW=1 then over the box draws a normal-view of *
	26	C. * volume NAME (with same seen attributes set by GDSPEC). *
	27	C. * If IONL<>1 then draw two other boxes, slightly shifted, *
	28	C. * behind the original one. *
	29	C. * If IDET=1 then draw another box, slightly larger, *
	30	C. * around the original one. *
	31	C. * If ISEEN=0 then draw the box(es) with dashed lines. *
	32	C. * *
	33	C. * ==>Called by : GDTREE *
	34	C. * Author : P.Zanarini ; S.Giani ********* *
	35	C. * *
	36	C. ******************************************************************
	37	C.
	38	#include "geant321/gcbank.inc"
	39	#include "geant321/gcnum.inc"
	40	#include "geant321/gcdraw.inc"
	41	DIMENSION X(5),Y(5)
	42	DIMENSION NAMT0(16),NAM0(4)
d43b40e2	43	DIMENSION PAR(100)
fe4da5cc	44	DIMENSION NNDM(14),INDM(4,14)
	45	DIMENSION ATT(10)
	46	CHARACTER*4 NAME,NAMSEE,ISON
	47	CHARACTER*14 CHNAM0
	48	SAVE NNDM,INDM
	49	COMMON/NPILEV/NPILE
	50	C
	51	DATA NNDM/0,0,0,4,0,2,0,2,4,3,4,3,0,0/
	52	DATA INDM/0,0,0,0, 0,0,0,0, 0,0,0,0, 2,3,7,11, 0,0,0,0,
	53	+ 4,5,0,0, 0,0,0,0, 6,7,0,0, 3,4,5,6, 4,5,6,0,
	54	+ 1,2,3,4, 1,2,3,0, 0,0,0,0, 0,0,0,0/
	55	C.
	56	C. ------------------------------------------------------------------
	57	C.
	58	IF (ISEEN.EQ.0) THEN
	59	LINSTY=4
	60	ELSE
	61	LINSTY=1
	62	ENDIF
	63	CALL MVBITS(LINSTY,0,3,LINATT,10)
	64	C
	65	C Compute box vertex
	66	C
	67	SIZ2=SIZE/2.
	68	SIZ4=SIZE/4.
	69	SIZ8=SIZE/8.
	70	*** X(1)=XLEV-SIZ2
	71	X(1)=XLEV-(SIZ22.FACHV)
	72	X(2)=X(1)
	73	*** X(3)=XLEV+SIZ2
	74	X(3)=XLEV+(SIZ22.FACHV)
	75	X(4)=X(3)
	76	X(5)=X(1)
	77	Y(1)=YLEV+SIZ2
	78	Y(2)=YLEV-SIZ2
	79	Y(3)=Y(2)
	80	Y(4)=Y(1)
	81	Y(5)=Y(1)
	82	C
	83	C Draw box
	84	C
	85	IWTY=IGIWTY(1)
	86	CALL IGQWK(0,'HWCO',R)
	87	IF(IWTY.GT.10.OR.IWTY.LT.1.OR.R.EQ.2.)THEN
	88	LINCOL=1
	89	CALL MVBITS(LINCOL,0,8,LINATT,16)
	90	CALL GDRAWV(X,Y,5)
	91	ELSE
	92	CALL IGSET('LWID',1.)
	93	CALL UCTOH(NAME,IIIN,4,4)
	94	CALL IGPID(1,'Tree',IIIN,' ')
	95	CALL IGPAVE(X(1),X(3),Y(1),Y(3),SIZ8,1004,1001,'TRS')
	96	ENDIF
	97	C
	98	C If volume is also detector draw another box outside
	99	C
	100	IF (IDET.EQ.1) THEN
	101	X(1)=X(1)-SIZ8
	102	X(2)=X(2)-SIZ8
	103	X(3)=X(3)+SIZ8
	104	X(4)=X(4)+SIZ8
	105	X(5)=X(5)-SIZ8
	106	Y(1)=Y(1)+SIZ8
	107	Y(2)=Y(2)-SIZ8
108	Y(3)=Y(3)-SIZ8
109	Y(4)=Y(4)+SIZ8
110	Y(5)=Y(5)+SIZ8
111	LINCOL=1
112	CALL MVBITS(LINCOL,0,8,LINATT,16)
113	CALL GDRAWV(X,Y,5)
114	ENDIF
115	C
116	C If volume is not 'ONLY' draw two other boxes behind
117	C
118	IF (IONL.NE.1) THEN
119	X(1)=XLEV+(SIZ22.FACHV)-SIZ8
120	Y(1)=YLEV+SIZ2
121	X(2)=X(1)
122	Y(2)=Y(1)+SIZ8
123	X(3)=X(2)-(SIZE2.FACHV)
124	Y(3)=Y(2)
125	X(4)=X(3)
126	Y(4)=Y(3)-SIZE
127	X(5)=X(4)+SIZ8
128	Y(5)=Y(4)
129	LINCOL=1
130	CALL MVBITS(LINCOL,0,8,LINATT,16)
131	CALL GDRAWV(X,Y,5)
132	DO 10 I=1,5
133	X(I)=X(I)-SIZ8
134	Y(I)=Y(I)+SIZ8
135	10 CONTINUE
136	CALL GDRAWV(X,Y,5)
137	ENDIF
138	C
139	LINSTY=1
140	CALL MVBITS(LINSTY,0,3,LINATT,10)
141	C
142	C Draw name
143	C
144	YL=YLEV-SIZ2+(SIZ8/2.)
145	IF(FACHV.EQ.1.)YL=YLEV-SIZ8
146	XL=XLEV-(SIZ4*FACHV)+SIZ8
147	IF(FACHV.EQ.1.)XL=XLEV
148	SZ=1.5SIZ42.
149	IF(FACHV.EQ.1.)SZ=1.SIZ42.
150	ITXCOL=5
151	CALL MVBITS(ITXCOL,0,4,ITXATT,3)
152	AITXCO=ITXCOL
153	CALL IGSET('TXCI',AITXCO)
154	CALL IGSET('TXFP',-60.)
155	IWTY=IGIWTY(1)
156	IF(IWTY.GT.10.OR.IWTY.LT.1)CALL IGSET('TXFP',-61.)
157	*** CALL IGPID(1,'Tree',NPILE,' ')
158	CALL GDRAWT(XL,YL,NAME,SZ,0.,4,0)
159	CALL IGSET('TXFP',2.)
160	CALL IGSET('TXCI',1.)
161	C
162	IF (IDUP.EQ.1) GO TO 20
163	C
164	C Draw IDUP in bottom half box
165	C
166	CALL UCTOH('<"4#>',NAMT0(1),1,5)
167	IF (IDUP.LT.0) CALL UCTOH('F',NAMT0(3),1,1)
168	DUP=FLOAT(IABS(IDUP))
169	CALL HBCDF(DUP,8,NAMT0(6))
170	CALL UCTOH('$',NAMT0(14),1,1)
171	CALL UBUNCH(NAMT0,NAM0,14)
172	CALL UHTOC(NAM0,4,CHNAM0,14)
173	XL=XLEV+(SIZ4*FACHV)+SIZ2
174	IF(FACHV.EQ.1.)XL=XLEV-(3.*SIZ8)
175	YL=YLEV-SIZ4-SIZ8
176	IF(FACHV.EQ.1.)YL=YLEV-SIZ2+(SIZ8/2.)
177	SZ=0.7SIZ42.
178	IF(FACHV.EQ.1.)SZ=0.7SIZ41.5
179	ITXCOL=5
180	CALL MVBITS(ITXCOL,0,4,ITXATT,3)
181	AITXCO=ITXCOL
182	CALL IGSET('TXCI',AITXCO)
183	CALL IGSET('TXFP',2.)
184	*** CALL IGPID(1,'Tree',NPILE,' ')
185	CALL GDRAWT(XL,YL,CHNAM0,SZ,0.,10,-1)
186	CALL IGSET('TXCI',1.)
187	C
188	20 CONTINUE
189	C
190	IF (IDRAW.NE.1) GO TO 999
191	C
192	C Draw a normal view of the volume
193	C
194	CALL GLOOK(NAME,IQ(JVOLUM+1),NVOLUM,IVO)
195	IF(IVO.LE.0) GO TO 999
196	C
197	JVO=LQ(JVOLUM-IVO)
198	C
199	C GET PARAMETERS
200	C
201	CALL GFPARA(NAME,1,0,NPAR,NATT,PAR,ATT)
202	IF(NPAR.LE.0) GO TO 999
203	C
204	ISHAPE=Q(JVO+2)
205	C
206	C CHECK parametER SIZES
207	C
208	PARMAX=-1.
209	DO 50 I=1,NPAR
210	IF(ISHAPE.GT.14)GO TO 50
211	IF(NNDM(ISHAPE).LE.0) GO TO 40
212	C
213	NDM=NNDM(ISHAPE)
214	DO 30 IDM=1,NDM
215	C
216	IF(I.EQ.INDM(IDM,ISHAPE)) GO TO 50
217	C
218	30 CONTINUE
219	40 CONTINUE
220	C
221	ABSPAR=ABS(PAR(I))
222	PARMAX=MAX(PARMAX,ABSPAR)
223	50 CONTINUE
224	XSCALE=SIZE/(3.*PARMAX)
225	YSCALE=XSCALE
226	XPOS=XLEV-(SIZEFACHV2)+SIZ8
227	YPOS=YLEV
228	* YPOS=YLEV+SIZE
229	C
230	CALL GFATT(NAME,'SEEN',KSEEN)
231	C
232	C Add local value KSEEN 1 to starting node of tree
233	C
234	KSEEN=KSEEN+110
235	ISEEN=KSEEN
236	CALL GSATT(NAME,'SEEN',ISEEN)
237	C
238	CALL GDNSON(NAME,NSON,IDIV)
239	DO 60 N=1,NSON
240	CALL GDSON(N,NAME,ISON)
241	CALL GFATT(ISON,'SEEN',KSEEN)
242	C
243	C ISON is a volume with multeplicity;
244	C first occurrence has already been set
245	C
246	IF (KSEEN.GT.50) GO TO 60
247	C
248	C Add local value SEEN -2 to each one-level-down node
249	C
250	KSEEN=KSEEN+80
251	ISEEN=KSEEN
252	CALL GSATT(ISON,'SEEN',ISEEN)
253	60 CONTINUE
254	CALL GDRAW(NAME,45.,135.,0.,XPOS,YPOS,XSCALE,YSCALE)
255	C
256	C Reset global SEEN values
257	C
258	DO 70 IVO=1,NVOLUM
259	CALL UHTOC(IQ(JVOLUM+IVO),4,NAMSEE,4)
260	CALL GFATT(NAMSEE,'SEEN',KSEEN)
261	IF (KSEEN.LT.50) GO TO 70
262	ISEENL=KSEEN/10.+0.5
263	ISEENG=KSEEN-ISEENL*10.
264	CALL GSATT(NAMSEE,'SEEN',ISEENG)
265	70 CONTINUE
266	C
267	999 END