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

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