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

*
* $Id$
*
* $Log$
* Revision 1.1.1.1  1995/10/24 10:20:21  cernlib
* Geant
*
*
#include "geant321/pilot.h"
*CMZ :  3.21/02 29/03/94  15.41.26  by  S.Giani
*-- Author :
      SUBROUTINE GDCHIT (IUSET,IUDET,ITRS,ISYMB,SIZMAX,KDHIT,
     +                   HITMIN,HITMAX)
C.
C.    ******************************************************************
C.    *                                                                *
C.    *       Draw Calorimetry Hits of track ITRS in detector IUDET    *
C.    *       of set IUSET, in case that the spatial coordinates       *
C.    *       of the HITS had been stored as the first three words     *
C.    *       in the argument vector of GSAHIT.                        *
C.    *                                                                *
C.    *       The following characters may be chosen :                 *
C.    *          ISYMB = -1   hardware character    (fast),            *
C.    *                =  0   crosses X by TVDRAW   (default),         *
C.    *                       (size = SSYMB in cm, default = 0.1),     *
C.    *                = 840  characters by HPLSOF                     *
C.    *                   .   according to the table p. 25             *
C.    *                   .   of the HPLOT manual                      *
C.    *                  853  (size = SSYMB in cm)                     *
C.    *             (844,854  white, black asterisk * in addition).    *
C.    *                                                                *
C.    *       SIZMAX is the maximum character size                     *
C.    *                                                                *
C.    *       KDHIT is the hit number selected : the character size    *
C.    *       is a function of HITS(KDHIT) and HITMIN and HITMAX are   *
C.    *       its limits                                               *
C.    *                                                                *
C.    *                                                                *
C.    *       JH =LQ(JHITS-ISET)                                       *
C.    *       JHD=LQ(JH-IDET)                                          *
C.    *       IQ(JH+IDET)=pointer to the LAST USED word in JHD         *
C.    *         Each hit was packed into JHD in the following format:  *
C.    *        --Track number ITRA not packed;                         *
C.    *        --Volume numbers packed;                                *
C.    *        --Hits transformed and packed.                          *
C.    *                                                                *
C.    *    ==>Called by : <USER>, <GXINT>                              *
C.    *       Authors : R.Brun, P.Zanarini   *********                 *
C.    *                                                                *
C.    ******************************************************************
C.
#include "geant321/gcbank.inc"
#include "geant321/gcdraw.inc"
      DIMENSION HITS(3,10),ASIZE(10)
      CHARACTER*4 IUSET,IUDET
C.
C.    ------------------------------------------------------------------
C.
      CALL UCTOH('OFF ',IOFF,4,4)
C
C            Set IOBJ to HIT
C
      IOBJ=3
C
      IF (JHITS.LE.0)           GOTO 140
      NSET = IQ(JSET-1)
      NS1        = 1
      NS2        = NSET
      IF (IUSET(1:1).EQ.'*')   GOTO  10
      CALL GLOOK(IUSET,IQ(JSET+1),NSET,ISET)
      IF (ISET .LE.0)           GOTO 140
      NS1        = ISET
      NS2        = ISET
   10 CONTINUE
C
C            loop on all selected sets
C
      DO 130   ISET=NS1,NS2
      JS = LQ(JSET-ISET)
      JH = LQ(JHITS-ISET)
      IF(JS.LE.0)GO TO 130
      IF(JH.LE.0)GO TO 130
      NDET = IQ(JS-1)
      ND1        = 1
      ND2        = NDET
      IF (IUDET(1:1).EQ.'*')   GOTO  20
      CALL GLOOK(IUDET,IQ(JS+1),NDET,IDET)
      IF (IDET .EQ.0)   GOTO 130
      ND1        = IDET
      ND2        = IDET
   20 CONTINUE
C
C            loop on selected detectors for this set
C
      DO 120   IDET=ND1,ND2
      JD = LQ(JS-IDET)
      JHD = LQ(JH-IDET)
      IF (JHD  .LE.0)   GOTO 120
      JDH = LQ(JD-1)
      ILAST = IQ(JH+IDET)
      IF (ILAST.EQ.0)   GOTO 120
      NV = IQ(JD+2)
      NH = IQ(JD+4)
      NW = IQ(JD+1) + IQ(JD+3) + 1
C
C            Now loop on all hits to get their space coordinates
C
      MH         = 0
      DO 110   I=1,ILAST,NW
      ITRA = IQ(JHD+I)
C
C            Is this the right track ?
C
      IF (ITRS.NE.0 .AND. ITRS.NE.ITRA) GOTO 110
C
C            Is THRZ option set on (R-Z projection)
C            and the cut lateral (vertical or horizontal) ?
C
      IF (ITHRZ.NE.IOFF) THEN
        IF (ICUT.EQ.1.OR.ICUT.EQ.2) ITR3D=-ITRA
      ELSE
        ITR3D=0
      ENDIF
C
C            Skip words with volume numbers
C
      NK         = 1
      K          = 1
      IF (NV.LE.0)      GOTO  40
      DO 30   IV=1,NV
      NB = IQ(JD+2*IV+10)
      IF (NB.EQ.0)  NB=32
      K          = NB + K
      IF (K.LE.33)      GOTO  30
      NK         = NK + 1
      K          = NB + 1
   30 CONTINUE
      NK         = NK + 1
C
C            Get unpacked (first three !) HITS.
C            Hits origin is shifted. Division by scale factor.
C
   40 CONTINUE
      IF (NH.LE.2)      GOTO 110
      MH         = MH + 1
C
      IF (KDHIT.GT.3.AND.KDHIT.LE.10) THEN
        KDH=KDHIT
      ELSE
        KDH=3
        HITKH=0.
      ENDIF
C
      K=1
      DO 80   IH=1,KDH
      NB = IQ(JDH+4*IH-2)
      IF (NB.GT.0)      GOTO  50
      IF (K.GT.1) THEN
        NK=NK+1
      ENDIF
      KHIT = IQ(JHD+I+NK)
      K          = 1
      NK         = NK + 1
      GOTO  70
   50 IF (K+NB.LE.33)   GOTO  60
      K          = 1
      NK         = NK + 1
   60 KHIT = IBITS(IQ(JHD+I+NK),K-1,NB)
      K          = K + NB
   70 CONTINUE
      IF (IH.LE.3) THEN
        HITS(IH,MH)= FLOAT(KHIT)/Q(JDH+4*IH) - Q(JDH+4*IH-1)
      ELSE
        HITKH=FLOAT(KHIT)/Q(JDH+4*IH) - Q(JDH+4*IH-1)
      ENDIF
   80 CONTINUE
C
C             Compute symbol size
C
      IF (HITMAX.LE.0) THEN
        SSYMB=SIZMAX
      ELSE
        SSYMB=SIZMAX*(HITKH-HITMIN)/HITMAX
        IF (SSYMB.GT.SIZMAX) SSYMB=SIZMAX
C
C             HPLSOF accepts a minimum SIZE of 0.001,
C             otherwise set SIZE=0.25 : protect here this case
C
        IF (SSYMB.LT.0.0011) SSYMB=0.0011
      ENDIF
      ASIZE(MH)=SSYMB
C
C            Draw the hits (usually 10 at a time);
C            flush HITS vector if
C             --loop on this detector volume is to finish,
C             --loop on current track is to finish,
C             --HITS vector is full.
C
      IF ((ILAST-I).LT.NW)   GOTO 90
      ITRN = IQ(JHD+I+NW)
      IF (ITRA.NE.ITRN)      GOTO 90
      IF (MH  .NE.  10)      GOTO 110
C
   90 CONTINUE
C
      DO 100 II=1,MH
        CALL GDAHIT(HITS(1,II),HITS(2,II),HITS(3,II),ISYMB,ASIZE(II))
  100 CONTINUE
      MH=0
C
  110 CONTINUE
  120 CONTINUE
  130 CONTINUE
  140 CONTINUE
C
C             Reset ITR3D to 'standard projection'
C
      ITR3D=0
      IOBJ=0
C
      END
Commit	Line	Data
fe4da5cc	1	*
	2	* $Id$
	3	*
	4	* $Log$
	5	* Revision 1.1.1.1 1995/10/24 10:20:21 cernlib
	6	* Geant
	7	*
	8	*
	9	#include "geant321/pilot.h"
	10	*CMZ : 3.21/02 29/03/94 15.41.26 by S.Giani
	11	*-- Author :
	12	SUBROUTINE GDCHIT (IUSET,IUDET,ITRS,ISYMB,SIZMAX,KDHIT,
	13	+ HITMIN,HITMAX)
	14	C.
	15	C. ******************************************************************
	16	C. * *
	17	C. * Draw Calorimetry Hits of track ITRS in detector IUDET *
	18	C. * of set IUSET, in case that the spatial coordinates *
	19	C. * of the HITS had been stored as the first three words *
	20	C. * in the argument vector of GSAHIT. *
	21	C. * *
	22	C. * The following characters may be chosen : *
	23	C. * ISYMB = -1 hardware character (fast), *
	24	C. * = 0 crosses X by TVDRAW (default), *
	25	C. * (size = SSYMB in cm, default = 0.1), *
	26	C. * = 840 characters by HPLSOF *
	27	C. * . according to the table p. 25 *
	28	C. * . of the HPLOT manual *
	29	C. * 853 (size = SSYMB in cm) *
	30	C. * (844,854 white, black asterisk * in addition). *
	31	C. * *
	32	C. * SIZMAX is the maximum character size *
	33	C. * *
	34	C. * KDHIT is the hit number selected : the character size *
	35	C. * is a function of HITS(KDHIT) and HITMIN and HITMAX are *
	36	C. * its limits *
	37	C. * *
	38	C. * *
	39	C. * JH =LQ(JHITS-ISET) *
	40	C. * JHD=LQ(JH-IDET) *
	41	C. * IQ(JH+IDET)=pointer to the LAST USED word in JHD *
	42	C. * Each hit was packed into JHD in the following format: *
	43	C. * --Track number ITRA not packed; *
	44	C. * --Volume numbers packed; *
	45	C. * --Hits transformed and packed. *
	46	C. * *
	47	C. * ==>Called by : <USER>, <GXINT> *
	48	C. * Authors : R.Brun, P.Zanarini ********* *
	49	C. * *
	50	C. ******************************************************************
	51	C.
	52	#include "geant321/gcbank.inc"
	53	#include "geant321/gcdraw.inc"
	54	DIMENSION HITS(3,10),ASIZE(10)
	55	CHARACTER*4 IUSET,IUDET
	56	C.
	57	C. ------------------------------------------------------------------
	58	C.
	59	CALL UCTOH('OFF ',IOFF,4,4)
	60	C
	61	C Set IOBJ to HIT
	62	C
	63	IOBJ=3
	64	C
65	IF (JHITS.LE.0) GOTO 140
66	NSET = IQ(JSET-1)
67	NS1 = 1
68	NS2 = NSET
69	IF (IUSET(1:1).EQ.'*') GOTO 10
70	CALL GLOOK(IUSET,IQ(JSET+1),NSET,ISET)
71	IF (ISET .LE.0) GOTO 140
72	NS1 = ISET
73	NS2 = ISET
74	10 CONTINUE
75	C
76	C loop on all selected sets
77	C
78	DO 130 ISET=NS1,NS2
79	JS = LQ(JSET-ISET)
80	JH = LQ(JHITS-ISET)
81	IF(JS.LE.0)GO TO 130
82	IF(JH.LE.0)GO TO 130
83	NDET = IQ(JS-1)
84	ND1 = 1
85	ND2 = NDET
86	IF (IUDET(1:1).EQ.'*') GOTO 20
87	CALL GLOOK(IUDET,IQ(JS+1),NDET,IDET)
88	IF (IDET .EQ.0) GOTO 130
89	ND1 = IDET
90	ND2 = IDET
91	20 CONTINUE
92	C
93	C loop on selected detectors for this set
94	C
95	DO 120 IDET=ND1,ND2
96	JD = LQ(JS-IDET)
97	JHD = LQ(JH-IDET)
98	IF (JHD .LE.0) GOTO 120
99	JDH = LQ(JD-1)
100	ILAST = IQ(JH+IDET)
101	IF (ILAST.EQ.0) GOTO 120
102	NV = IQ(JD+2)
103	NH = IQ(JD+4)
104	NW = IQ(JD+1) + IQ(JD+3) + 1
105	C
106	C Now loop on all hits to get their space coordinates
107	C
108	MH = 0
109	DO 110 I=1,ILAST,NW
110	ITRA = IQ(JHD+I)
111	C
112	C Is this the right track ?
113	C
114	IF (ITRS.NE.0 .AND. ITRS.NE.ITRA) GOTO 110
115	C
116	C Is THRZ option set on (R-Z projection)
117	C and the cut lateral (vertical or horizontal) ?
118	C
119	IF (ITHRZ.NE.IOFF) THEN
120	IF (ICUT.EQ.1.OR.ICUT.EQ.2) ITR3D=-ITRA
121	ELSE
122	ITR3D=0
123	ENDIF
124	C
125	C Skip words with volume numbers
126	C
127	NK = 1
128	K = 1
129	IF (NV.LE.0) GOTO 40
130	DO 30 IV=1,NV
131	NB = IQ(JD+2*IV+10)
132	IF (NB.EQ.0) NB=32
133	K = NB + K
134	IF (K.LE.33) GOTO 30
135	NK = NK + 1
136	K = NB + 1
137	30 CONTINUE
138	NK = NK + 1
139	C
140	C Get unpacked (first three !) HITS.
141	C Hits origin is shifted. Division by scale factor.
142	C
143	40 CONTINUE
144	IF (NH.LE.2) GOTO 110
145	MH = MH + 1
146	C
147	IF (KDHIT.GT.3.AND.KDHIT.LE.10) THEN
148	KDH=KDHIT
149	ELSE
150	KDH=3
151	HITKH=0.
152	ENDIF
153	C
154	K=1
155	DO 80 IH=1,KDH
156	NB = IQ(JDH+4*IH-2)
157	IF (NB.GT.0) GOTO 50
158	IF (K.GT.1) THEN
159	NK=NK+1
160	ENDIF
161	KHIT = IQ(JHD+I+NK)
162	K = 1
163	NK = NK + 1
164	GOTO 70
165	50 IF (K+NB.LE.33) GOTO 60
166	K = 1
167	NK = NK + 1
168	60 KHIT = IBITS(IQ(JHD+I+NK),K-1,NB)
169	K = K + NB
170	70 CONTINUE
171	IF (IH.LE.3) THEN
172	HITS(IH,MH)= FLOAT(KHIT)/Q(JDH+4IH) - Q(JDH+4IH-1)
173	ELSE
174	HITKH=FLOAT(KHIT)/Q(JDH+4IH) - Q(JDH+4IH-1)
175	ENDIF
176	80 CONTINUE
177	C
178	C Compute symbol size
179	C
180	IF (HITMAX.LE.0) THEN
181	SSYMB=SIZMAX
182	ELSE
183	SSYMB=SIZMAX*(HITKH-HITMIN)/HITMAX
184	IF (SSYMB.GT.SIZMAX) SSYMB=SIZMAX
185	C
186	C HPLSOF accepts a minimum SIZE of 0.001,
187	C otherwise set SIZE=0.25 : protect here this case
188	C
189	IF (SSYMB.LT.0.0011) SSYMB=0.0011
190	ENDIF
191	ASIZE(MH)=SSYMB
192	C
193	C Draw the hits (usually 10 at a time);
194	C flush HITS vector if
195	C --loop on this detector volume is to finish,
196	C --loop on current track is to finish,
197	C --HITS vector is full.
198	C
199	IF ((ILAST-I).LT.NW) GOTO 90
200	ITRN = IQ(JHD+I+NW)
201	IF (ITRA.NE.ITRN) GOTO 90
202	IF (MH .NE. 10) GOTO 110
203	C
204	90 CONTINUE
205	C
206	DO 100 II=1,MH
207	CALL GDAHIT(HITS(1,II),HITS(2,II),HITS(3,II),ISYMB,ASIZE(II))
208	100 CONTINUE
209	MH=0
210	C
211	110 CONTINUE
212	120 CONTINUE
213	130 CONTINUE
214	140 CONTINUE
215	C
216	C Reset ITR3D to 'standard projection'
217	C
218	ITR3D=0
219	IOBJ=0
220	C
221	END