5 * Revision 1.1.1.1 1995/10/24 10:19:40 cernlib
9 *CMZ : 3.21/02 29/03/94 15.41.35 by S.Giani
13 GEANT specific commands.
20 The hidden line removal technique is necessary to visualize properly
21 very complex detectors. At the same time, it can be useful to visualize
22 the inner elements of a detector in detail. For this purpose, the
23 commands menu CVOL has been developed: these commands allow
24 subtractions (via boolean operation) of given shapes from any part of
25 the detector, therefore showing its inner contents. It is possible
26 to clip each different volume by means of a different shape (BOX ,
27 TUBE, CONE, SPHE are available). If '*' is given as the name of the
28 volume to be clipped, all volumes are clipped by the given shape.
29 A volume can be clipped at most twice (even by
30 different shapes); if a volume is explicitely clipped
31 twice, the '*' will not act on it anymore. Giving '.' as the name
32 of the volume to be clipped will reset the clipping.
36 CNNV ' Name of volume to be clipped ' C D='* '
38 XMIN ' Lower limit of the Shape X coordinate ' R D=-10000.
39 XMAX ' Upper limit of the Shape X coordinate ' R D=-9999.
40 YMIN ' Lower limit of the Shape Y coordinate ' R D=-10000.
41 YMAX ' Upper limit of the Shape Y coordinate ' R D=-9999.
42 ZMIN ' Lower limit of the Shape Z coordinate ' R D=-10000.
43 ZMAX ' Upper limit of the Shape Z coordinate ' R D=-9999.
45 This command performs a boolean subtraction between the volume
46 CNVV and a box placed in the MARS according the values of the given
47 coordinates. See also CVOL.
48 The following commands will clip by a box,
49 with a vertex at the origin, the volume specified by NAME (a valid
50 string for the NAME of the volume can be found using the DTREE command).
54 draw NAME 40 40 0 10 10 .01 .01
56 box NAME 0 1000 0 1000 0 1000
57 draw NAME 40 40 0 10 10 .01 .01
64 CNVV ' Name of volume to be clipped ' C D='* '
66 RMAX ' External radius of tube ' R D=0.1
67 ZDEM ' Half length of tube axis ' R D=0.1
68 XMED ' Center X coordinate ' R D=-10000.
69 YMED ' Center Y coordinate ' R D=-10000.
70 ZMED ' Center Z coordinate ' R D=-10000.
72 This command performs a boolean subtraction between the volume
73 CNVV and a tube; the tube has the given parameters and is placed in
74 the MARS according the given coordinates of its center.
76 The following commands will clip, by a tube,
77 positioned according to the given parameters, the volume specified
78 by NAME (a valid string for the NAME of the volume
79 can be found using the DTREE command).
83 draw NAME 40 40 0 10 10 .01 .01
85 tube * 500 1000 500 0 0
86 draw NAME 40 40 0 10 10 .01 .01
93 CNVV ' Name of volume to be clipped ' C D='* '
95 RMAX1 ' Min external radius ' R D=0.1
96 RMAX2 ' Max external radius ' R D=0.1
97 ZDEM ' Half length of cone axis ' R D=0.1
98 XMED ' Center X coordinate ' R D=-10000.
99 YMED ' Center Y coordinate ' R D=-10000.
100 ZMED ' Center Z coordinate ' R D=-10000.
102 This command performs a boolean subtraction between the volume
103 CNVV and a cone; the cone has the given parameters and is placed in
104 the MARS according to the given coordinates of its center.
106 The following commands will clip by a cone,
107 positioned according the given parameters, the volume specified
108 by NAME (a valid string for the NAME of the volume
109 can be found using the DTREE command).
113 draw NAME 40 40 0 10 10 .01 .01
115 cone * 1 750 1000 0 0 1000
116 draw NAME 40 40 0 10 10 .01 .01
123 CNVV ' Name of volume to be clipped ' C D='* '
125 RMAX ' External radius of sphere ' R D=0.1
126 XMED ' Center X coordinate ' R D=-10000.
127 YMED ' Center Y coordinate ' R D=-10000.
128 ZMED ' Center Z coordinate ' R D=-10000.
130 This command performs a boolean subtraction between the volume
131 CNVV and a sphere; the sphere has the given parameters and is placed in
132 the MARS according to the given coordinates of its center.
133 See also CVOL. The following commands clip by a sphere,
134 positioned according to the given parameters, the volume specified
135 by NAME (a valid string for the NAME of the volume
136 can be found using the DTREE command).
140 draw NAME 40 40 0 10 10 .01 .01
143 draw NAME 40 40 0 10 10 .01 .01
150 XCUT 'x coordinate of cutted value' R D=0.
151 YCUT 'y coordinate of cutted value' R D=0.
152 ZCUT 'z coordinate of cutted value' R D=0.
154 It allows the cutting in the ray-tracing. All the volumes are cutted
155 from XCUT to +BIG along the x axis, from YCUT to +BIG along the y axis
156 and from ZCUT to +BIG along the z axis.
162 Drawing commands. These commands allow the visualization in several ways
163 of the volumes defined in the geometrical data structure. It is possible
164 to draw the logical tree of volumes belonging to the detector (DTREE),
165 to show their geometrical specification (DSPEC,DFSPC), to draw them
166 and their cut views (DRAW, DCUT). Moreover, it is possible to execute
167 these commands when the hidden line removal option is activated; in
168 this case, the volumes can be also either translated in the space
169 (SHIFT), or clipped by boolean operation (CVOL). In addition, it is
170 possible to fill the surfaces of the volumes
171 with solid colours when the shading option (SHAD) is activated.
172 Several tools (ZOOM, LENS) have been developed to zoom detailed parts
173 of the detectors or to scan physical events as well.
174 Finally, the command MOVE will allow the rotation, translation and zooming
175 on real time parts of the detectors or tracks and hits of a simulated event.
176 Ray-tracing commands. In case the command (DOPT RAYT ON) is executed,
177 the drawing is performed by the Geant ray-tracing;
178 automatically, the color is assigned according to the tracking medium of each
179 volume and the volumes with a density lower/equal than the air are considered
180 transparent; if the option (USER) is set (ON) (again via the command (DOPT)),
181 the user can set color and visibility for the desired volumes via the command
182 (SATT), as usual, relatively to the attributes (COLO) and (SEEN).
183 The resolution can be set via the command (SATT * FILL VALUE), where (VALUE)
184 is the ratio between the number of pixels drawn and 20 (user coordinates).
185 Parallel view and perspective view are possible (DOPT PROJ PARA/PERS); in the
186 first case, we assume that the first mother volume of the tree is a box with
187 dimensions 10000 X 10000 X 10000 cm and the view point (infinetely far) is
188 5000 cm far from the origin along the Z axis of the user coordinates; in the
189 second case, the distance between the observer and the origin of the world
190 reference system is set in cm by the command (PERSP NAME VALUE); grand-angle
191 or telescopic effects can be achieved changing the scale factors in the command
192 (DRAW). When the final picture does not occupy the full window,
193 mapping the space before tracing can speed up the drawing, but can also
194 produce less precise results; values from 1 to 4 are allowed in the command
195 (DOPT MAPP VALUE), the mapping being more precise for increasing (VALUE); for
196 (VALUE = 0) no mapping is performed (therefore max precision and lowest speed).
197 The command (VALCUT) allows the cutting of the detector by three planes
198 ortogonal to the x,y,z axis. The attribute (LSTY) can be set by the command
199 SATT for any desired volume and can assume values from 0 to 7; it determines
200 the different light processing to be performed for different materials:
201 0 = dark-matt, 1 = bright-matt, 2 = plastic, 3 = ceramic, 4 = rough-metals,
202 5 = shiny-metals, 6 = glass, 7 = mirror. The detector is assumed to be in the
203 dark, the ambient light luminosity is 0.2 for each basic hue (the saturation
204 is 0.9) and the observer is assumed to have a light source (therefore he will
205 produce parallel light in the case of parallel view and point-like-source
206 light in the case of perspective view).
212 THETA 'Viewing angle theta (for 3D projection)' R R=0.:180.
213 PHI 'Viewing angle phi (for 3D projection)' R R=0.:360.
214 PSI 'Viewing angle psi (for 2D rotation)' R R=0.:360.
215 U0 'U-coord. (horizontal) of volume origin' R
216 V0 'V-coord. (vertical) of volume origin' R
217 SU 'Scale factor for U-coord.' R
218 SV 'Scale factor for V-coord.' R
220 CALL GDRAW(name,theta,phi,psi,u0,v0,su,sv)
221 If optional parameters are missing, the corresponding values are
222 taken from the common /GCDRAW/. This command will draw the volumes,
223 selected with their graphical attributes, set by the SATT
224 facility. The drawing may be performed with hidden line removal
225 and with shading effects according to the value of the options HIDE
226 and SHAD; if the option SHAD is ON, the contour's edges can be
227 drawn or not. If the option HIDE is ON, the detector can be
228 exploded (BOMB), clipped with different shapes (CVOL), and some
229 of its parts can be shifted from their original
230 position (SHIFT). When HIDE is ON, if
231 the drawing requires more than the available memory, the program
232 will evaluate and display the number of missing words
233 (so that the user can increase the
234 size of its ZEBRA store). Finally, at the end of each drawing (with HIDE on),
235 the program will print messages about the memory used and
236 statistics on the volumes' visibility.
237 The following commands will produce the drawing of a green
238 volume, specified by NAME, without using the hidden line removal
239 technique, using the hidden line removal technique,
240 with different linewidth and colour (red), with
241 solid colour, with shading of surfaces, and without edges.
242 Finally, some examples are given for the ray-tracing. (A possible
243 string for the NAME of the volume can be found using the command DTREE).
247 draw NAME 40 40 0 10 10 .01 .01
250 draw NAME 40 40 0 10 10 .01 .01
254 draw NAME 40 40 0 10 10 .01 .01
259 draw NAME 40 40 0 10 10 .01 .01
262 draw NAME 40 40 0 10 10 .01 .01
265 draw NAME 40 40 0 10 10 .01 .01
269 draw NAME 40 40 0 10 10 .01 .01
272 draw NAME 40 40 0 10 10 1 1
278 draw NAME 40 40 0 10 10 5 5
284 XLPOS 'x coordinate of light source' R
285 YLPOS 'y coordinate of light source' R
286 ZLPOS 'z coordinate of light source' R
287 INTEN 'intensity of light source' I
289 This point-like light source can be moved in the space and its intensity
290 can be changed (INTEN going from 0 to 10) relatively to the ambience light.
295 TSEQTO 'total sequential time' R
296 NPROC 'number of processors' I
297 NMPTOT 'number of message passing' I
298 TOTMBY 'total megabytes transfert' R
299 TSEQ 'not parallelized code' R
300 TLAT 'latency time' R
301 TNET 'network speed in Mbytes/sec' R
303 It sets the values of the parameters expressed in the formula and
304 specify which variables must be assumed as x,y,z (setting their value
305 to 1001,1002,1003, respectively).
310 X1MIN 'x coordinate min' R
311 X1MAX 'x coordinate max' R
312 Y1MIN 'y coordinate min' R
313 Y1MAX 'y coordinate max' R
314 Z1MIN 'z coordinate min' R
315 Z1MAX 'z coordinate max' R
317 It sets the range for the x,y,z variables.
322 N 'Number of elements in arrays LNAMES and LNUMBS' I D=1
323 NAMNUM 'Volume names and numbers (ex. "NAME1,NR1,NAME2,NR2")' C
324 CHNRS 'Reference system used' C D='MARS' R='MARS,DRS'
326 THETA 'Viewing angle theta (for 3D projection)' R R=0.:360.
327 PHI 'Viewing angle phi (for 3D projection)' R R=0.:360.
328 PSI 'Viewing angle psi (for 2D rotation)' R R=0.:180.
329 U0 'U-coord. (horizontal) of volume origin' R
330 V0 'V-coord. (vertical) of volume origin' R
331 SU 'Scale factor for U-coord.' R
332 SV 'Scale factor for V-coord.' R
334 CALL GDRVOL(n,lnames,lnumbs,nrs,theta,phi,psi,u0,v0,su,sv)
335 N is the number of levels from the top of the geometry structure
336 to the volume lnames(n),lnumbs(n) to be drawn.
337 NAMNUM contain the arrays lnames and lnumbs,
338 identifying the path, in pairs and separated by commas; for
340 'lname(1),lnumbs(1),lname(2),lnumbs(2) '
341 CHNRS is the name of the reference system used: MARS for MAster Reference
342 System or DRS for Daughter Reference System.
343 NRS=0 for MARS or NRS<>0 for DRS
344 If optional parameters are missing, the current values in /GCDRAW/
351 CAXIS 'Axis value' C R='X,Y,Z'
352 CUTVAL 'Cut plane distance from the origin along the axis' R
354 U0 'U-coord. (horizontal) of volume origin' R
355 V0 'V-coord. (vertical) of volume origin' R
356 SU 'Scale factor for U-coord.' R
357 SV 'Scale factor for V-coord.' R
359 CALL GDRAWC(name,iaxis,cutval,u0,v0,su,sv)
360 The cut plane is normal to caxis (X,Y,Z), corresponding to iaxis (1,2,3),
361 and placed at the distance cutval from the origin.
362 The resulting picture is seen from the the same axis.
363 If optional parameters are missing, the current values in /GCDRAW/
365 When HIDE Mode is ON, it is possible to get the same effect with
366 the CVOL/BOX command.
372 CUTTHE 'Theta angle of the line normal to cut plane' R R=0.:360.
373 CUTPHI 'Phi angle of the line normal to cut plane' R R=0.:360.
374 CUTVAL 'Cut plane distance from the origin along the axis' R
376 THETA 'Viewing angle theta (for 3D projection)' R R=0.:360.
377 PHI 'Viewing angle phi (for 3D projection)' R R=0.:360.
378 U0 'U-coord. (horizontal) of volume origin' R
379 V0 'V-coord. (vertical) of volume origin' R
380 SU 'Scale factor for U-coord.' R
381 SV 'Scale factor for V-coord.' R
383 CALL GDRAWX(name,cutthe,cutphi,cutval,theta,phi,u0,v0,su,sv)
384 The cut plane is normal to the line given by the cut angles
385 cutthe and cutphi and placed at the distance cutval from the origin.
386 The resulting picture is seen from the viewing angles theta,phi.
387 If optional parameters are missing, the current values in /GCDRAW/
393 CNVN ' Name of volume to be shifted ' C D='*'
394 XXXX ' Shift along X axis ' R D=0.
395 YYYY ' Shift along Y axis ' R D=0.
396 ZZZZ ' Shift along Z axis ' R D=0.
398 To draw a volume shifted from its initial position when hidden
399 line removal is ON. It can be useful if you want to extract a
400 volume or some volumes from the detector to show them more clearly.
401 The last requested SHIFT for each volume
402 NAME is performed. Moreover, the SHIFT of
403 each volume will be performed starting from where its mother has
404 been shifted, so that it's easier to SHIFT nicely sets
405 of volumes using the mother-daughter relationships.
406 If '.' is given as the name of the volume
407 to be shifted, the shifts for all volumes will be reset.
408 The following commands will produce the translation along
409 the Z-axis of the previously drawn volume:
413 draw NAME 40 40 0 10 10 .01 .01
420 BOOM ' Exploding factor for volumes position ' R D=0. R=-10.:10.
422 To 'explode' the detector. If BOOM is positive (values smaller
423 than 1. are suggested, but any value is possible)
424 all the volumes are shifted by a distance
425 proportional to BOOM along the direction between their centre
426 and the origin of the MARS; the volumes which are symmetric
427 with respect to this origin are simply not shown.
428 BOOM equal to 0 resets the normal mode.
429 A negative (greater than -1.) value of
430 BOOM will cause an 'implosion'; for even lower values of BOOM
431 the volumes' positions will be reflected respect to the origin.
432 This command can be useful to improve the 3D effect for very
433 complex detectors. The following commands will make explode the
438 draw NAME 40 40 0 10 10 .01 .01
441 draw NAME 40 40 0 10 10 .01 .01
448 NAME 'Volume name' C D=' '
449 LEVMAX 'Depth level' I D=3 R=-15:15
450 ISELT 'Options ' I D=111
452 This command allows the drawing of the logical tree,
453 displaying the name, the multiplicity and other information about the volumes,
454 via a call to GDTREE(name,levmax,isel):
455 if the third parameter is not given (default), the command will
456 produce the drawing of the tree displaying, for each volume, the
457 number of the following levels (red arrows) and of the preceeding
458 levels (green arrows); then the control is automatically given to the
459 mouse: clicking on the left button when the cursor is inside a volume's
460 pave will perform a DSPEC for that volume; doing the same when the cursor
461 is on a red arrow, will perform a DTREE for the relative volume (the
462 number of levels displayed depending on the clicked arrow); doing the
463 same for the 'i-th' green arrow of a given volume, will perform a DTREE
464 for its mother-volume staying 'i' levels before.
465 If running with X-windows, the drawing of the specification (DSPEC)
467 in a different window to speed up the scanning of the tree.
468 Iterating this procedure it is possible to analyse very easily and quickly
469 any kind of tree. Clicking the right button of the mouse will return
470 the control to the command mode.
471 If the ISELT parameter is given,
472 then the TREE will work as in the
473 previous version, with ISELT up to 10001.
474 The following command will perform a drawing of the tree and give the
475 control to the user via the mouse:
485 Trough a call to GDSPEC(name), this command allows one to show three
486 views of the volume (two cut-views and a 3D view), together with
487 its geometrical specifications. The 3D drawing will
488 be performed according the current values of the options HIDE and
489 SHAD and according the current CVOL clipping parameters for that
497 TETA3 'Theta angle' R D=40. R=0.:180.
498 PHI3 'Phi angle' R D=40. R=0.:360.
499 PSI3 'Psi angle' R D=0. R=0.:360.
500 U03 'U-coord. (horizontal) of volume origin' R D=10. R=-40.:40.
501 V03 'V-coord. (vertical) of volume origin' R D=10. R=-40.:40.
502 ZM3 'Zoom factor for current size factors' R D=1. R=0.00001:10.
504 Trough a call to GSPE3D, this command allows one to show
505 the volume (3D views in real time), together with
506 its geometrical specifications (if using MOTIF). The 3D drawing will
507 be performed according the current values of the options HIDE and
508 SHAD and according the current CVOL clipping parameters for that
516 CSORT 'Alphabetic sorting flag' C D='N' R='Y,N,0,1'
517 CINTER 'Interactive/Batch version' C D='I' R='I,B,0,1'
519 CALL GDFSPC(name,isort,inter)
520 Same as DSPEC, but it will draw the specifications for all the volumes.
521 If the alphabetic sorting flag is YES, all pictures will be drawn in ascending
522 alphabetic order; isort is set to 1.
523 If INTERACTIVE, (inter=1), the routine will prompt the user at each plot
524 before doing a clear screen, otherwise it will clear automatically
525 the screen before starting a new frame.
530 X0 'X-coord. (horizontal) of text string' R D=10. R=0.:20.
531 Y0 'Y-coord. (vertical) of text string' R D=10. R=0.:20.
532 TEXT 'Text string' C D='GEANT'
533 SIZE 'Character size (cm)' R D=.5
534 ANGLE 'Rotation angle (deg)' R D=0. R=0.:360.
535 LWID 'Line width' I D=4
536 CENT 'Centering option' C D='CENT' R='CENT,LEFT,RIGH'
538 CALL GDRAWT(x0,y0,text,size,angle,lwid,opt)
539 It allows one to draw some text in the current picture.
540 Now more than 160 colours are available. The text colour
541 must be set via the command IGSET. The size of the
542 text will follow the zooming factors in the view banks.
547 XVECT 'Vector containing X-coord. (horizontal)' C
548 YVECT 'Vector containing Y-coord. (vertical)' C
549 NPOINT 'Number of coord.' I
551 Draw a polyline of 'npoint' point via
552 a call to GDRAWV(xvect,yvect,npoint)
553 where xvect and yvect are two KUIP vectors
558 U 'U-coord. (horizontal) of the centre of scale' R
559 V 'V-coord. (vertical) of the centre of scale' R
562 It draws a scale centered in U,V.
567 X0 'X-coord. of axis origin' R
568 Y0 'Y-coord. of axis origin' R
569 Z0 'Z-coord. of axis origin' R
572 CALL GDAXIS(x0,y0,z0,dx)
573 This commmand superimposes the axis of the MARS on the
574 current picture. It is useful for finding immediately the
575 orientation of the current drawing of the detector in the space.
580 U 'U-coord. (horizontal) of the centre of man' R
581 V 'V-coord. (vertical) of the centre of man' R
582 TYPE 'Man, Wm1, Wm2, Wm3' C D='MAN' R='MAN,WM1,WM2,WM3'
584 CALL GDMAN(u,v),CALL GDWMN1(u,v),CALL GDWMN2(u,v),CALL GDWMN2(u,v)
585 It superimposes the picure of a man or of a woman, chosen among
586 three different ones, with the same scale factors as the detector
587 in the current drawing.
593 ISEL 'Option flag' I D=111110
595 CHRSIZ 'Character size (cm) of title NAME' R D=0.6
597 CALL GDHEAD(isel,name,chrsiz)
599 0 to have only the header lines
600 xxxxx1 to add the text name centered on top of header
601 xxxx1x to add global detector name (first volume) on left
602 xxx1xx to add date on right
603 xx1xxx to select thick characters for text on top of header
604 x1xxxx to add the text 'EVENT NR x' on top of header
605 1xxxxx to add the text 'RUN NR x' on top of header
606 NOTE that ISEL=x1xxx1 or ISEL=1xxxx1 are illegal choices,
607 i.e. they generate overwritten text.
609 and CHRSIZ the character size in cm of text name.
614 Position the cursor on the first point (u1,v1) and hit the space bar(GKS).
615 Position the cursor on the second point (u2,v2) and hit the space bar(GKS).
616 Clicking the left button of the mouse (X11) will have the same effect as
617 hiting the space bar (GKS).
618 The command will compute and print the distance in space separating
619 the two points on the projection view. It can be useful to measure
620 distances either between volumes or between tracks or hits.
626 Activates graphic input to identify detector elements
627 in a cut view. Clicking on the left button of the mouse when
628 the cursor is in a given point of the drawing and clicking again
629 (outside the detector) will produce the following effect:
630 a line joininig the two points will be drawn together with
631 the name and the medium number of the volume picked
632 with the first clicking close to the second point.
637 NAME 'Volume name' C D=' '
639 NOPT 'S=sample mode,T=tracks,H=hits' C D=' '
641 Positioning some daughter volumes inside a 'mother', it can be
642 important to check if overlaps between such volumes have occurred.
643 Instead of putting the drawing in a view bank, zooming, and iterating
644 the process for different viewing angles of the same detector, the
645 MOVE facility has been developed (for machines running with X11):
646 it is sufficient to draw a view of the volumes to be analysed (after
647 setting the proper SEEN, COLO, etc. attributes) and then to enter
648 'MOVE' followed by the same 'NAME' used for the last command DRAW.
649 The detector will appear in a panel with five buttons at the
650 bottom: THETA, PHI, TRASL, ZOOM, OFF. Clicking on the left button
651 of the mouse, when the cursor is inside the THETA area, will rotate the
652 detector along the polar angle theta according to the
653 backward-to-forward movement of the mouse
654 (clicking up and down the left button if
655 not in sample mode); clicking on the right button of
656 the mouse will stop the rotation; clicking now on the
657 left button of the mouse when inside the PHI area will activate a
658 rotation along the polar angle phi. In the same way, activating the
659 TRASL button, the detector can be translated in the u,v plane
660 of the screen according to the 2D-movement of the mouse. Finally,
661 activating the ZOOM button, the detector will be zoomed (or unzoomed)
662 according to the backward-to-forward movement of the mouse. Clicking on the
663 OFF button will return the control to the 'command mode'. The MOVE
664 command will work also with hidden line removal and shading options
665 (when SHAD is on the background will be black);
666 moreover, if the volumes are clipped, exploded, shifted, etc., they
667 will be 'MOVED' with these features as well.
668 Tracks and hits of a previously stored physical event can be moved
669 together with the detector, allowing a dynamical 3-D analysis of the
670 simulated events. Clicking the central button of the mouse when a good
671 view of the event is found, will stop any movement and the mouse will
672 allow the normal picking capabilities first for the tracks and then for
673 the hits. After clicking of the right button, the normal
674 movement will restart to find another interesting view of the event
675 and to iterate the process.
676 The MOVE is also available in sample mode.
677 The following commands will produce a drawing of a volume
678 and then will give the control to the MOVE panel; try the following
683 draw NAME 40 40 0 10 10 .01 .01
688 draw NAME 40 40 0 10 10 .01 .01
695 draw NAME 40 40 0 10 10 .01 .01
702 NAME 'Volume name' C D=' '
704 THETA 'Viewing angle theta (for 3D projection)' R D=40. R=0.:180.
705 PHI 'Viewing angle phi (for 3D projection)' R D=40. R=0.:360.
706 PSI 'Viewing angle psi (for 2D rotation)' R D=0. R=0.:180.
707 U0 'U-coord. (horizontal) of volume origin' R D=10. R=0.:20.
708 V0 'V-coord. (vertical) of volume origin' R D=10. R=0.:20.
709 SU 'Scale factor for U-coord.' R D=0.01
710 SV 'Scale factor for V-coord.' R D=0.01
711 SZ 'Scale zoom factor' R D=1. R=0.1:10.
712 NOPT 'T=tracks,H=hits' C D=' ' R='T,H'
714 Same functionality of the command MOVE interfaced with MOTIF.
719 NAME 'Volume name' C D=' '
720 DISTT 'Volume distance from observer' R D=1000.
722 SAMP 'Control to the mouse' C D='OFF '
724 To control the perspective according to the variation of the distance
725 between the observer and the object (if PROJ has the value PERS).
726 If SAMP is ON the control of the distance is given via the mouse.
731 KNUM 'View bank identifier' I D=1
733 KSAM 'Sample mode ' C D='OFF '
735 Interactive zooming for detectors and events when running
736 with X-windows. Using this command, when showing the contents of a
737 view bank, it is possible to click (left button) in two points of the
738 drawing (which will represent the left upper corner and the right
739 bottom corner of the part to be zoomed). After the second click
740 a new 'window' will appear to fit the frame defined
741 by the two clicks and it will show a zoomed view as seen from a
742 lens with those dimensions. Clicking now the central button will
743 translate the lens over the drawing, while clicking the right button
744 will stop it. Moreover, clicking the left button of the
745 mouse, the lens will increase (or decrease) its magnification
746 power according to the backward-to-forward movement of the mouse.
747 A click on the right button will stop this action and it is possible
748 to restart the translation of the lens or, clicking
749 on the right button again, to make the lens disappear. It is then possible
750 to open another 'window-lens' with different dimensions. Thus,
751 this command can be useful to scan detailed parts of a detector or
752 to scan hits and showers for events. Clicking the right
753 button when no lens is displayed will return the control to the
754 'command mode'. The LENS is also available in sample mode when KSAM is
756 The following commands will fill a view bank and will
757 allow to scan the detector and an event previously stored
758 via the use of LENS (when running
763 draw NAME 40 40 0 10 10 .01 .01
775 ZFU 'Zoom factor for U-coord. (horizontal)' R D=2.
776 ZFV 'Zoom factor for V-coord. (vertical)' R D=2.
778 UZ0 'U-coord. of the centre of zoom rectangle' R R=0.:20. D=10.
779 VZ0 'V-coord. of the centre of zoom rectangle' R R=0.:20. D=10.
780 U0 'U-coord. of the centre of resulting zoomed rectangle' R R=0.:20. D=10.
781 V0 'V-coord. of the centre of resulting zoomed rectangle' R R=0.:20. D=10.
783 CALL GDZOOM(zfu,zfv,uz0,vz0,u0,v0)
784 This command sets the zoom parameters that will be used by
785 subsequent calls to the drawing routines. Each zoom operation is always
786 relative to the status of the current zoom parameters.
787 The scale factors in u,v are respectively zfu,zfv.
788 zfu=0 (or zfv=0) will act as a reset (i.e. unzoomed viewing).
789 The zoom is computed around uz0,vz0 (user coordinates),
790 and the resulting picture will be centered at u0,v0.
791 The use of the space bar is replaced by the left button of the mouse
795 1. position the cursor at (uz0,vz0)
796 2. type the space bar (GKS)
797 (u0,v0 are chosen at centre of screen)
800 1. position the cursor at first corner of zoom rectangle
801 2. type the space bar (GKS)
802 3. position the cursor at second corner of zoom rectangle
803 4. type the space bar (GKS)
804 (zfu,zfv are chosen according to the zoom rectangle;
805 uz0,vz0 are chosen at the centre of the zoom rectangle;
806 u0,v0 are chosen at centre of screen)
809 1. position the cursor at (uz0,vz0)
810 2. type the space bar (GKS)
811 3. position the cursor at (u0,v0)
812 4. type the space bar (GKS)
814 If isel=1000+n and running with X-windows:
815 1. n must be the identifier of an active view bank
816 2. clicking on the left button of the mouse will display
817 a zoomed view (computed around the cursor position) of
818 the previous drawing in a new window
819 3. it is now possible to iterate the zooming from the new window
820 4. clicking on the right button will return the control to the
822 5. clicking on the left button it is possible to open new windows
823 zooming in other points of the detector
824 6. clicking on the right button when the main window is active
825 will return the control to the 'command mode'.
831 ITRA 'Track number' I D=0
834 Draw tracks previously stored via GSXYZ.
840 EPSILO 'Delta angle' R D=0.25
843 The picking of track points requires the JXYZ data structure
844 and is repeated until the character typed is 'Q' or 'q' (GKS)
845 or the right button of the mouse is clicked (X11).
846 EPSILO is the delta angle used for picking; if EPSILO=0
847 there is no optimization performed and
848 over all the track points the one nearest to the pick
855 ITRA 'Track number' I D=0
856 ISEL 'Option flag' I D=11
857 SIZE 'Character size (cm) for particle names' R D=0.25
859 CALL GDPART(itra,isel,size)
860 isel=x1 to draw the track number
861 isel=1x to draw the particle name
867 CHUSET 'User set identifier' C D='*'
868 CHUDET 'User detector identifier' C D='*'
869 ITRA 'Number of the selected track' I D=0
870 ISYMB 'Character selection number' I D=0
871 SSYMB 'Size of characters (cm)' R D=0.1
873 CALL GDHITS(chuset,chudet,itra,isymb,ssymb).
874 The character plotted at each hit point may be chosen by isymb :
875 -1 (small) hardware points (fast)
876 0 software crosses (default)
877 840,850 empty/full circles (slow)
878 841,851 empty/full squares (slow)
879 842,852 empty/full triangles (up) (slow)
880 843,853 empty diamond/full triangle (down) (slow)
881 844,854 empty/full stars (slow)
882 Except for isymb=-1, the size of the character on the screen can be
883 chosen by SSYMB cm. The hit colour will follow the value of TXCI (text
884 colour) for isymb>0, the value of PMCI (polymarkers colour) for isymb<0,
885 the value of PLCI (polyline colour) for isymb=0.
891 CHUSET 'User set identifier' C D='*'
892 CHUDET 'User detector identifier' C D='*'
893 EPSILO 'Pick aperture' R D=0.1
895 CALL GKHITS(chuset,chudet,epsilo)
896 The picking of hit points requires the appropriate JSET data structure
898 and is repeated until the character typed is 'Q' or 'q' (GKS) or the
899 right button of the mouse is clicked (X11).
900 If the character typed to pick is 'K' or 'k' then the
901 kinematics of the corresponding track is also printed.
902 The search is made of all the hits of all tracks in
903 detector CHUDET of set CHUSET.
904 EPSILO is the pick aperture; if EPSILO<0 its absolute value is taken
905 and in addition the pick aperture is drawn; if EPSILO=0
906 there is an infinite pick aperture and
907 over all the hits the one nearest to the pick point is taken.
913 CHUSET 'User set identifier' C D='*'
914 CHUDET 'User detector identifier' C D='*'
915 ITRA 'Number of the selected track' I D=0
916 ISYMB 'Character selection number' I D=0
917 SIZMAX 'Maximum character size (cm)' R D=1
918 IHIT 'Index of array HITS' I D=4
919 HITMIN 'Lower boundary of HITS(IHIT)' R D=0
920 HITMAX 'Upper boundary of HITS(IHIT)' R D=0
922 CALL GDCHIT(chuset,chudet,itra,isymb,sizmax,ihit,hitmin,hitmax)
923 The character plotted at each hit point may be chosen via
924 CSYMB; isymb is composed as:
925 -1 (small) hardware points (fast)
926 0 software crosses (default)
927 840,850 empty/full circles (slow)
928 841,851 empty/full squares (slow)
929 842,852 empty/full triangles (up) (slow)
930 843,853 empty diamond/full triangle (down) (slow)
931 844,854 empty/full stars (slow)
932 Except for isymb=-1 the SIZE of the character on the screen
933 is a function of HITS(IHIT), the array containing the calorimeter
934 quantity, with HITMIN and HITMAX defining its range.
935 The maximum character size (used in overflow) is SIZMAX.
936 SIZE = SIZMAX * ( HITS(IHIT) - HITMIN ) / HITMAX
941 NAME 'Detector name' C
943 CPXTYP 'Complexity name' C
945 IVIEW 'View number where picture is stored' I D=0
947 CALL GUVIEW(name,type,cpxtyp,iview)
952 >Menu /GEANT/GRAPHICS_CONTROL
954 Graphics control commands.
958 IVIEW 'View number' I
961 When a drawing is very complex and requires a long time to be
962 executed, it can be useful to store it in a view bank: after a
963 call to DOPEN and the execution of the drawing (nothing will
964 appear on the screen), and after a necessary call to DCLOSE,
965 the contents of the bank can be displayed in a very fast way
966 through a call to DSHOW; therefore, the detector can be easily
967 zoomed many times in different ways. Please note that the pictures
968 with solid colours can now be stored in a view bank or in 'PICTURE FILES'.
974 IVIEW 'View number' I
977 It shows on the screen the contents of a view bank. It
978 can be called after a view bank has been closed.
983 IVIEW 'View number' I
986 It deletes a view bank from memory.
992 It closes the currently open view bank; it must be called after the
993 end of the drawing to be stored.
999 It open a new workstation (if not already opened) and activate it
1000 (deactivating the default one).
1006 It deactivate the previosly activated workstation and reactivate the
1012 IOPT 'Attribute name' C
1014 Set current attribute.
1020 RVAL 'Attribute value' R D=1. R=-10.:10.
1022 Set current attribute value.
1028 NAME 'Volume name' C D='* '
1029 IOPT 'Name of the attribute to be set' C D='DEFA'
1030 IVAL 'Value to which the attribute is to be set' I D=10000
1032 CALL GSATT(name,iopt,ival)
1033 name='*' stands for all the volumes.
1034 iopt can be chosen among the following :
1036 'WORK' 0=volume name is inactive for the tracking
1037 1=volume name is active for the tracking (default)
1039 'SEEN' 0=volume name is invisible
1040 1=volume name is visible (default)
1041 -1=volume invisible with all its descendants in the tree
1042 -2=volume visible but not its descendants in the tree
1044 'LSTY' line style 1,2,3,... (default=1)
1045 LSTY=7 will produce a very precise approximation for
1048 'LWID' line width -7,...,1,2,3,..7 (default=1)
1049 LWID<0 will act as abs(LWID) was set for the volume
1050 and for all the levels below it. When SHAD is 'ON', LWID
1051 represent the linewidth of the scan lines filling the surfaces
1052 (whereas the FILL value represent their number). Therefore
1053 tuning this parameter will help to obtain the desired
1054 quality/performance ratio.
1056 'COLO' colour code -166,...,1,2,..166 (default=1)
1058 n=2=red; n=17+m, m=0,25, increasing luminosity according to 'm';
1059 n=3=green; n=67+m, m=0,25, increasing luminosity according to 'm';
1060 n=4=blue; n=117+m, m=0,25, increasing luminosity according to 'm';
1061 n=5=yellow; n=42+m, m=0,25, increasing luminosity according to 'm';
1062 n=6=violet; n=142+m, m=0,25, increasing luminosity according to 'm';
1063 n=7=lightblue; n=92+m, m=0,25, increasing luminosity according to 'm';
1064 colour=n*10+m, m=1,2,...9, will produce the same colour
1065 as 'n', but with increasing luminosity according to 'm';
1066 COLO<0 will act as if abs(COLO) was set for the volume
1067 and for all the levels below it.
1068 When for a volume the attribute FILL is > 1 (and the
1069 option SHAD is on), the ABS of its colour code must be < 8
1070 because an automatic shading of its faces will be
1073 'FILL' (1992) fill area -7,...,0,1,...7 (default=0)
1074 when option SHAD is 'on' the FILL attribute of any
1075 volume can be set different from 0 (normal drawing);
1076 if it is set to 1, the faces of such volume will be filled
1077 with solid colours; if ABS(FILL) is > 1, then a light
1078 source is placed along the observer line, and the faces of
1079 such volumes will be painted by colours whose luminosity
1080 will depend on the amount of light reflected;
1081 if ABS(FILL) = 1, then it is possible to use all the 166
1082 colours of the colour table, becouse the automatic shading
1084 for increasing values of FILL the drawing will be performed
1085 with higher and higher resolution improving the quality (the
1086 number of scan lines used to fill the faces increases with FILL);
1087 it is possible to set different values of FILL
1088 for different volumes, in order to optimize at the same time
1089 the performance and the quality of the picture;
1090 FILL<0 will act as if abs(FILL) was set for the volume
1091 and for all the levels below it.
1092 This kind of drawing can be saved in 'picture files'
1094 0=drawing without fill area
1095 1=faces filled with solid colours and resolution = 6
1096 2=lowest resolution (very fast)
1097 3=default resolution
1102 Finally, if a coloured background is desired, the FILL
1103 attribute for the first volume of the tree must be set
1104 equal to -abs(colo), colo being >0 and <166.
1106 'SET ' set number associated to volume name
1107 'DET ' detector number associated to volume name
1108 'DTYP' detector type (1,2)
1113 GSCU 'Scale factor for U-coord.' R
1114 GSCV 'Scale factor for V-coord.' R
1116 Change the scale factors GSCU and GSCV in /GCDRAW/.
1121 ICOL 'Colour code' I D=1
1123 CALL GDCOL(-abs(icol))
1128 LWIDTH 'Line width code' I D=1
1130 CALL GDLW(-abs(lwidth))
1135 Clear screen (start a new picture on graphics file, if opened).
1141 IOPT 'Option name' C D='*'
1142 IVAL 'Option value' C D='*'
1144 CALL GDOPT(iopt,ival)
1145 To set/modify the drawing options.
1148 THRZ ON Draw tracks in R vs Z
1149 OFF (D) Draw tracks in X,Y,Z
1152 PROJ PARA (D) Parallel projection
1154 TRAK LINE (D) Trajectory drawn with lines
1155 POIN " " with markers
1156 HIDE ON Hidden line removal using the CG package
1157 OFF (D) No hidden line removal
1158 SHAD ON Fill area and shading of surfaces.
1159 OFF (D) Normal hidden line removal.
1160 RAYT ON Ray-tracing on.
1161 OFF (D) Ray-tracing off.
1162 EDGE OFF Does not draw contours when shad is on.
1163 ON (D) Normal shading.
1164 MAPP 1,2,3,4 Mapping before ray-tracing.
1166 USER ON User graphics options in the raytracing.
1167 OFF (D) Automatic graphics options.
1174 XSIZE 'Size along X' R D=20.
1175 YSIZE 'Size along Y' R D=20.
1177 Set the size of the picture.
1178 On the terminal, the pictures will have the ratio YSIZE/XSIZE, and,
1179 if a metafile is produced, pictures will be YSIZE by XSIZE cm.
1180 This command sets the parameters for the normalisation transformation
1181 number 1 to [0-XSIZE], [0-YSIZE].
1186 DPERS 'Distance from the origin' R
1188 Set the variable dpers in /GCDRAW/, representing
1189 the distance from the origin when using option PERSpective.
1195 ICADD 'Colour table index' I D=0
1196 ICVAL 'Colour table value' I D=0
1198 Sets the color table LOOKTB(ICADD)=ICVAL.
1199 If ICADD=0 then LOOKTB(1:16) is taken.
1200 If ICVAL is omitted the current value of LOOKTB(ICADD) is shown.
1211 LHSTA_1 'user word' C
1212 LHSTA_2 'user word' C
1213 LHSTA_3 'user word' C
1214 LHSTA_4 'user word' C
1215 LHSTA_5 'user word' C
1216 LHSTA_6 'user word' C
1217 LHSTA_7 'user word' C
1218 LHSTA_8 'user word' C
1219 LHSTA_9 'user word' C
1220 LHSTA_10 'user word' C
1221 LHSTA_11 'user word' C
1222 LHSTA_12 'user word' C
1223 LHSTA_13 'user word' C
1224 LHSTA_14 'user word' C
1225 LHSTA_15 'user word' C
1226 LHSTA_16 'user word' C
1227 LHSTA_17 'user word' C
1228 LHSTA_18 'user word' C
1229 LHSTA_19 'user word' C
1230 LHSTA_20 'user word' C
1232 The command HSTA is similar to the HSTA data records. It can accept
1233 up to 20 4-character words. If the first argument is '.', the number
1234 of words is reset to 0 and all the words to four blanks.
1240 LGET_1 'user word' C
1241 LGET_2 'user word' C
1242 LGET_3 'user word' C
1243 LGET_4 'user word' C
1244 LGET_5 'user word' C
1245 LGET_6 'user word' C
1246 LGET_7 'user word' C
1247 LGET_8 'user word' C
1248 LGET_9 'user word' C
1249 LGET_10 'user word' C
1250 LGET_11 'user word' C
1251 LGET_12 'user word' C
1252 LGET_13 'user word' C
1253 LGET_14 'user word' C
1254 LGET_15 'user word' C
1255 LGET_16 'user word' C
1256 LGET_17 'user word' C
1257 LGET_18 'user word' C
1258 LGET_19 'user word' C
1259 LGET_20 'user word' C
1261 The command GET is similar to the GET data records. It can accept
1262 up to 20 4-character words. If the first argument is '.', the number
1263 of words is reset to 0 and all the words to four blanks.
1269 LSAVE_1 'user word' C
1270 LSAVE_2 'user word' C
1271 LSAVE_3 'user word' C
1272 LSAVE_4 'user word' C
1273 LSAVE_5 'user word' C
1274 LSAVE_6 'user word' C
1275 LSAVE_7 'user word' C
1276 LSAVE_8 'user word' C
1277 LSAVE_9 'user word' C
1278 LSAVE_10 'user word' C
1279 LSAVE_11 'user word' C
1280 LSAVE_12 'user word' C
1281 LSAVE_13 'user word' C
1282 LSAVE_14 'user word' C
1283 LSAVE_15 'user word' C
1284 LSAVE_16 'user word' C
1285 LSAVE_17 'user word' C
1286 LSAVE_18 'user word' C
1287 LSAVE_19 'user word' C
1288 LSAVE_20 'user word' C
1290 The command SAVE is similar to the SAVE data records. It can accept
1291 up to 20 4-character words. If the first argument is '.', the number
1292 of words is reset to 0 and all the words to four blanks.
1298 LSETS_1 'user word' C
1299 LSETS_2 'user word' C
1300 LSETS_3 'user word' C
1301 LSETS_4 'user word' C
1302 LSETS_5 'user word' C
1303 LSETS_6 'user word' C
1304 LSETS_7 'user word' C
1305 LSETS_8 'user word' C
1306 LSETS_9 'user word' C
1307 LSETS_10 'user word' C
1308 LSETS_11 'user word' C
1309 LSETS_12 'user word' C
1310 LSETS_13 'user word' C
1311 LSETS_14 'user word' C
1312 LSETS_15 'user word' C
1313 LSETS_16 'user word' C
1314 LSETS_17 'user word' C
1315 LSETS_18 'user word' C
1316 LSETS_19 'user word' C
1317 LSETS_20 'user word' C
1319 The command SETS is similar to the SETS data records. It can accept
1320 up to 20 4-character words. If the first argument is '.', the number
1321 of words is reset to 0 and all the words to four blanks.
1327 LPRIN_1 'user word' C
1328 LPRIN_2 'user word' C
1329 LPRIN_3 'user word' C
1330 LPRIN_4 'user word' C
1331 LPRIN_5 'user word' C
1332 LPRIN_6 'user word' C
1333 LPRIN_7 'user word' C
1334 LPRIN_8 'user word' C
1335 LPRIN_9 'user word' C
1336 LPRIN_10 'user word' C
1337 LPRIN_11 'user word' C
1338 LPRIN_12 'user word' C
1339 LPRIN_13 'user word' C
1340 LPRIN_14 'user word' C
1341 LPRIN_15 'user word' C
1342 LPRIN_16 'user word' C
1343 LPRIN_17 'user word' C
1344 LPRIN_18 'user word' C
1345 LPRIN_19 'user word' C
1346 LPRIN_20 'user word' C
1348 The command PRIN is similar to the PRIN data records. It can accept
1349 up to 20 4-character words. If the first argument is '.', the number
1350 of words is reset to 0 and all the words to four blanks.
1356 LGEOM_1 'user word' C
1357 LGEOM_2 'user word' C
1358 LGEOM_3 'user word' C
1359 LGEOM_4 'user word' C
1360 LGEOM_5 'user word' C
1361 LGEOM_6 'user word' C
1362 LGEOM_7 'user word' C
1363 LGEOM_8 'user word' C
1364 LGEOM_9 'user word' C
1365 LGEOM_10 'user word' C
1366 LGEOM_11 'user word' C
1367 LGEOM_12 'user word' C
1368 LGEOM_13 'user word' C
1369 LGEOM_14 'user word' C
1370 LGEOM_15 'user word' C
1371 LGEOM_16 'user word' C
1372 LGEOM_17 'user word' C
1373 LGEOM_18 'user word' C
1374 LGEOM_19 'user word' C
1375 LGEOM_20 'user word' C
1377 The command GEOM is similar to the GEOM data records. It can accept
1378 up to 20 4-character words. If the first argument is '.', the number
1379 of words is reset to 0 and all the words to four blanks.
1385 LVIEW_1 'user word' C
1386 LVIEW_2 'user word' C
1387 LVIEW_3 'user word' C
1388 LVIEW_4 'user word' C
1389 LVIEW_5 'user word' C
1390 LVIEW_6 'user word' C
1391 LVIEW_7 'user word' C
1392 LVIEW_8 'user word' C
1393 LVIEW_9 'user word' C
1394 LVIEW_10 'user word' C
1395 LVIEW_11 'user word' C
1396 LVIEW_12 'user word' C
1397 LVIEW_13 'user word' C
1398 LVIEW_14 'user word' C
1399 LVIEW_15 'user word' C
1400 LVIEW_16 'user word' C
1401 LVIEW_17 'user word' C
1402 LVIEW_18 'user word' C
1403 LVIEW_19 'user word' C
1404 LVIEW_20 'user word' C
1406 The command VIEW is similar to the VIEW data records. It can accept
1407 up to 20 4-character words. If the first argument is '.', the number
1408 of words is reset to 0 and all the words to four blanks.
1414 LPLOT_1 'user word' C
1415 LPLOT_2 'user word' C
1416 LPLOT_3 'user word' C
1417 LPLOT_4 'user word' C
1418 LPLOT_5 'user word' C
1419 LPLOT_6 'user word' C
1420 LPLOT_7 'user word' C
1421 LPLOT_8 'user word' C
1422 LPLOT_9 'user word' C
1423 LPLOT_10 'user word' C
1424 LPLOT_11 'user word' C
1425 LPLOT_12 'user word' C
1426 LPLOT_13 'user word' C
1427 LPLOT_14 'user word' C
1428 LPLOT_15 'user word' C
1429 LPLOT_16 'user word' C
1430 LPLOT_17 'user word' C
1431 LPLOT_18 'user word' C
1432 LPLOT_19 'user word' C
1433 LPLOT_20 'user word' C
1435 The command PLOT is similar to the PLOT data records. It can accept
1436 up to 20 4-character words. If the first argument is '.', the number
1437 of words is reset to 0 and all the words to four blanks.
1443 LSTAT_1 'user word' C
1444 LSTAT_2 'user word' C
1445 LSTAT_3 'user word' C
1446 LSTAT_4 'user word' C
1447 LSTAT_5 'user word' C
1448 LSTAT_6 'user word' C
1449 LSTAT_7 'user word' C
1450 LSTAT_8 'user word' C
1451 LSTAT_9 'user word' C
1452 LSTAT_10 'user word' C
1453 LSTAT_11 'user word' C
1454 LSTAT_12 'user word' C
1455 LSTAT_13 'user word' C
1456 LSTAT_14 'user word' C
1457 LSTAT_15 'user word' C
1458 LSTAT_16 'user word' C
1459 LSTAT_17 'user word' C
1460 LSTAT_18 'user word' C
1461 LSTAT_19 'user word' C
1462 LSTAT_20 'user word' C
1464 The command STAT is similar to the STAT data records. It can accept
1465 up to 20 4-character words. If the first argument is '.', the number
1466 of words is reset to 0 and all the words to four blanks.
1472 LRGET_1 'user word' C
1473 LRGET_2 'user word' C
1474 LRGET_3 'user word' C
1475 LRGET_4 'user word' C
1476 LRGET_5 'user word' C
1477 LRGET_6 'user word' C
1478 LRGET_7 'user word' C
1479 LRGET_8 'user word' C
1480 LRGET_9 'user word' C
1481 LRGET_10 'user word' C
1482 LRGET_11 'user word' C
1483 LRGET_12 'user word' C
1484 LRGET_13 'user word' C
1485 LRGET_14 'user word' C
1486 LRGET_15 'user word' C
1487 LRGET_16 'user word' C
1488 LRGET_17 'user word' C
1489 LRGET_18 'user word' C
1490 LRGET_19 'user word' C
1491 LRGET_20 'user word' C
1493 The command RGET is similar to the RGET data records. It can accept
1494 up to 20 4-character words. If the first argument is '.', the number
1495 of words is reset to 0 and all the words to four blanks.
1501 LRSAVE_1 'user word' C
1502 LRSAVE_2 'user word' C
1503 LRSAVE_3 'user word' C
1504 LRSAVE_4 'user word' C
1505 LRSAVE_5 'user word' C
1506 LRSAVE_6 'user word' C
1507 LRSAVE_7 'user word' C
1508 LRSAVE_8 'user word' C
1509 LRSAVE_9 'user word' C
1510 LRSAVE_10 'user word' C
1511 LRSAVE_11 'user word' C
1512 LRSAVE_12 'user word' C
1513 LRSAVE_13 'user word' C
1514 LRSAVE_14 'user word' C
1515 LRSAVE_15 'user word' C
1516 LRSAVE_16 'user word' C
1517 LRSAVE_17 'user word' C
1518 LRSAVE_18 'user word' C
1519 LRSAVE_19 'user word' C
1520 LRSAVE_20 'user word' C
1522 The command RSAV is similar to the RSAV data records. It can accept
1523 up to 20 4-character words. If the first argument is '.', the number
1524 of words is reset to 0 and all the words to four blanks.
1528 >Menu /GEANT/GEOMETRY
1534 IOPTI 'GSORD optimisation level' I D=0 R=-1,2
1536 This flag controls the tracking optimisation performed via the
1538 1 no optimisation at all; GSORD calls disabled;
1539 0 no optimisation; only user calls to GSORD kept;
1540 1 all non-GSORDered volumes are ordered along the best axis;
1541 2 all volumes are ordered along the best axis.
1546 NAME 'Volume name' C
1547 SHAPE 'Volume type' C
1548 NUMED 'Tracking medium number' I
1549 NPAR 'Number of shape parameters' I
1550 PAR 'Vector containing shape parameters' C
1552 CALL GSVOLU(name,shape,numed,par,npar,ivolu)
1553 where par is a KUIP vector.
1554 It creates a new volume in the JVOLUM data structure.
1559 NAME 'Volume name' C
1560 NUMBER 'Copy number of the volume' I
1561 MOTHER 'Mother volume name' C
1562 X0 'X coord. of the volume in mother ref. sys.' R
1563 Y0 'Y coord. of the volume in mother ref. sys.' R
1564 Z0 'Z coord. of the volume in mother ref. sys.' R
1565 IROT 'Rotation matrix number w.r.t. mother ref. sys.' I
1566 ONLY 'ONLY/MANY flag' C
1568 CALL GSPOS(name,number,mother,x0,y0,z0,irot,only)
1569 It positions a previously defined volume in the mother.
1574 NAME 'Volume name' C
1575 MOTHER 'Mother volume name' C
1576 NDIV 'Number of divisions' I
1577 CAXIS 'Axis value' C R='X,Y,Z,1,2,3'
1579 CALL GSDVN(name,mother,ndiv,iaxis)
1580 X,Y,Z of CAXIS will be translated to 1,2,3 for IAXIS.
1581 It divides a previously defined volume.
1589 Prints volumes' specifications.
1594 IROT 'Rotation matrix number' I
1595 THETA1 'Polar angle for axis I' R D=0. R=0.:180.
1596 PHI1 'Azimuthal angle for axis I' R D=0. R=0.:360.
1597 THETA2 'Polar angle for axis II' R D=0. R=0.:180.
1598 PHI2 'Azimuthal angle for axis II' R D=0. R=0.:360.
1599 THETA3 'Polar angle for axis III' R D=0. R=0.:180.
1600 PHI3 'Azimuthal angle for axis III' R D=0. R=0.:360.
1602 CALL GSROTM(irot,theta1,phi1,theta2,phi2,theta3,phi3)
1603 It defines the rotation matrix number IROT.
1611 Print matrixes' specifications.
1617 NTMED 'Tracking medium number' I D=1
1618 NAME 'Tracking medium name' C
1619 NMAT 'Material number' I D=1
1620 ISVOL 'Sensitive volume flag' I D=0
1621 IFIELD 'Magnetic field' I D=0
1622 FIELDM 'Max. field value (Kilogauss)' R D=0
1623 TMAXFD 'Max. angle due to field (deg/step)' R D=0.01
1624 STEMAX 'Max. step allowed' R D=1.E+10
1625 DEEMAX 'Max. fraction of energy lost in a step' R D=0.01
1626 EPSIL 'Tracking precision (cm)' R D=0.01
1627 STMIN 'Min. step due to continuos processes (cm)' R D=0.1
1629 CALL GSTMED(ntmed,name,nmat,isvol,ifield,fieldm,tmaxfd,
1630 + stemax,deemax,epsil,stmin,0,0)
1631 IFIELD = 0 if no magnetic field; IFIELD = -1 if user decision in GUSWIM;
1632 IFIELD = 1 if tracking performed with GRKUTA; IFIELD = 2 if tracking
1633 performed with GHELIX; IFIELD = 3 if tracking performed with GHELX3.
1641 Print tracking media's specifications.
1647 ISEL 'Options' I D=0
1648 NAME 'Volume name' C D=' '
1650 CALL GEDITV(isel,name)
1651 When the routine prompts for input parameters that do not need
1652 to be changed, type return.
1653 ISEL is used to select the editing operation to be performed:
1655 ISEL=1, to modify shape parameters PAR given by GSVOLU
1656 ISEL=2, to modify NAME given by GSVOLU
1657 ISEL=3, to delete NAME given by GSVOLU
1658 ISEL=4, to unlink NAME,NR given by GSPOS/GSDVN/GSDV..
1659 ISEL=5, to modify X0,Y0,Z0 of NAME,NR given by GSPOS
1660 ISEL=6, to modify IROT of NAME,NR given by GSPOS
1661 ISEL=7, to modify NDIV given by GSDVN
1662 ISEL=8, to modify IAXIS given by GSDVN
1667 FNAME 'Name of the SET file' C D='example.set'
1668 ANAME 'Name of the volume' C
1669 NBINS 'Number of the instances' I D=1
1670 LUNIT 'Logical unit number for SET file' I D=66
1671 LUNIT2 'Logical unit number for material file' I D=67
1672 INST 'Name of your institute' C D='CERN'
1673 SITE 'Name of site' C D='MEYRIN'
1674 DEPT 'Name of departement' C D='CN'
1675 RESP 'Name of sender' C D='god_knows_who'
1677 CALL GTXSET(fname,aname,nbins,lunit,lunit2,inst,site,dept,resp)
1678 This command produces a SET file describing the given volume with
1679 the contents currently set visible. (Use the visibility attribute,
1680 see SATT SEEN.) The description is given as a flat assembly
1681 related to the global coordinate system.
1682 The ouput can be read into CAD systems (EUCLID-IS) trough a SET interface.
1683 A list of materials of the volumes in the SET file and the GEANT tree
1684 is written into a file with the same filename as the SET file,
1685 but with extension .mat.
1690 LUN 'Logical unit of the file to be read' I R=1:100
1691 FNAME 'Name of the EUCLID file to be read' C
1693 CALL GREUCL(LUN,FNAME)
1694 Calls the routine to read into GEANT a geometry from an ASCII file
1695 written by the EUCLID-GEANT interface.
1700 LUN 'Logical unit of the file to be written' I R=1:100
1701 FNAME 'Name of the EUCLID file to be written' C
1702 TOPVOL 'Volume name of the starting node' C
1704 NUMBER 'Copy number of TOPVOL (relevant for GSPOSP)' I D=1
1705 NLEVEL 'Number of levels in the tree structure' I D=15
1708 CALL GWEUCL(LUN,FNAME)
1709 Calls the routine to write the current GEANT geometry into an ASCII file
1710 in EUCLID compatible format.
1715 It creates volumes of the given shape interactively.
1716 CALL GSVOLU(name,shape,numed,par,npar,ivolu)
1717 where par is a KUIP vector
1721 NAME 'Volume name' C
1722 NUMED 'Tracking medium number' I
1723 HALFX 'Half X length' R
1724 HALFY 'Half Y length' R
1725 HALFZ 'Half Z length' R
1727 YESNO 'GSPOSP option' C D='NO' R='YES,NO'
1733 NAME 'Volume name' C
1734 NUMED 'Tracking medium number' I
1735 HLFDWX 'Half X length in Lower Z Surface' R
1736 HLFUPX 'Half X length in Upper Z Surface' R
1737 HALFY 'Half Y length' R
1738 HALFZ 'Half Z length' R
1740 YESNO 'GSPOSP option' C D='NO' R='YES,NO'
1747 NAME 'Volume name' C
1748 NUMED 'Tracking medium number' I
1749 HLFDWX 'Half X length in Lower Z Surface' R
1750 HLFUPX 'Half X length in Upper Z Surface' R
1751 HLFDWY 'Half Y length in Lower Z Surface' R
1752 HLFUPY 'Half Y length in Upper Z Surface' R
1753 HALFZ 'Half Z length' R
1755 YESNO 'GSPOSP option' C D='NO' R='YES,NO'
1762 NAME 'Volume name' C
1763 NUMED 'Tracking medium number' I
1764 INRAD 'Inside Radius' R
1765 OUTRAD 'Outside Radius' R
1766 HALFZ 'Half Z length' R
1768 YESNO 'GSPOSP option' C D='NO' R='YES,NO'
1775 NAME 'Volume name' C
1776 NUMED 'Tracking medium number' I
1777 INRAD 'Inside Radius' R
1778 OUTRAD 'Outside Radius' R
1779 HALFZ 'Half Z length' R
1780 SPHI 'Start of section PHI' R R=0.:360.
1781 EPHI 'End of section PHI' R R=0.:360.
1783 YESNO 'GSPOSP option' C D='NO' R='YES,NO'
1790 NAME 'Volume name' C
1791 NUMED 'Tracking medium number' I
1792 INRDW 'Inside Radius in Lower Z Surface' R
1793 OUTRDW 'Outside Radius in Lower Z Surface' R
1794 INRUP 'Inside Radius in Upper Z Surface' R
1795 OUTRUP 'Outside Radius in Upper Z Surface' R
1796 HALFZ 'Half Z length' R
1798 YESNO 'GSPOSP option' C D='NO' R='YES,NO'
1805 NAME 'Volume name' C
1806 NUMED 'Tracking medium number' I
1807 INRDW 'Inside Radius in Lower Z Surface' R
1808 OUTRDW 'Outside Radius in Lower Z Surface' R
1809 INRUP 'Inside Radius in Upper Z Surface' R
1810 OUTRUP 'Outside Radius in Upper Z Surface' R
1811 HALFZ 'Half Z length' R
1812 SPHI 'Start of section PHI' R R=0.:360.
1813 EPHI 'End of section PHI' R R=0.:360.
1815 YESNO 'GSPOSP option' C D='NO' R='YES,NO'
1822 NAME 'Volume name' C
1823 NUMED 'Tracking medium number' I
1824 INRAD 'Inside Radius' R
1825 OUTRAD 'Outside Radius' R
1826 SPHI 'Start of section PHI' R R=0.:360.
1827 EPHI 'End of section PHI' R R=0.:360.
1828 STHETA 'Start of section THETA' R
1829 ETHETA 'End of section THETA' R
1831 YESNO 'GSPOSP option' C D='NO' R='YES,NO'
1838 NAME 'Volume name' C
1839 NUMED 'Tracking medium number' I
1840 HALFX 'Half X length' R
1841 HALFY 'Half Y length' R
1842 HALFZ 'Half Z length' R
1843 AXIS 'Angle of Y mid-faces segment to Y axis' R R=0.:360.
1844 PHI 'PHI angle of Low Z mid-face to High Z mid-face segment' R R=0.:360.
1845 THETA 'THETA angle of mid-low-Z-face to mid-high-Z-face segment' R R=0.:360.
1847 YESNO 'GSPOSP option' C D='NO' R='YES,NO'
1854 >Menu /GEANT/CONTROL
1871 PKINE10 'PKINE(10)' R
1873 Set the variables in /GCFLAG/ IKINE, PKINE(10)
1878 IDRUN 'User run number' I
1879 IDEVT 'User starting event number' I
1881 Set the run number and the starting value for the user event number.
1886 ISTORD 'Flag to control user ordering of the stack' I D=1 R=1,0
1888 If ISTORD is set to 1, the particle with the highest value of the
1889 user weight UPWGHT will be selected to be tracked next.
1894 TIMINT 'Total time after initialisation' R
1895 TIMEND 'Time reserved for the termination phase' R
1896 ITIME 'Frequency of control printing' I
1898 These commands have limited use in the interactive version. In
1899 particular the value of TIMINT is disregarded by GEANT.
1904 Restart tracking, clearing the track and hit
1905 banks, but keeping the kinematics.
1911 N 'Number of events' I D=1
1913 Start one or more new events.
1919 ISEED1 'First seed for the random number generator' I
1920 ISEED2 'Second seed for the random number generator' I
1922 Set the seeds for the random number generator. If no numbers are
1923 given, the currents seeds are printed.
1928 ISWI 'Switch number' I
1929 IVAL 'New switch value' I
1931 Change one element of array ISWIT(10) in /GCFLAG/
1937 LEVEL 'MZ log level' I D=0
1939 Set the log level for the MZ package of ZEBRA: CALL MZLOGL(0,level)
1940 LEVEL = -3 no messages at all
1941 -2 error messages only
1945 +2 log calls to MZ routines
1952 NUMBER 'Number' I D=0
1954 CALL GPRINT(name,number)
1959 LOUT 'New output unit' I
1961 To change lout in /GCUNIT/
1962 Note: unit numbers 5,11,12,13,14,15 are reserved and cannot be used.
1968 CHUSET 'User set' C D='*'
1969 CHUDET 'User detector' C D='*'
1970 NUMHI 'Hit number' I D=0
1972 CALL GPHITS(chuset,chudet)
1978 CHUSET 'User set' C D='*'
1979 CHUDET 'User detector' C D='*'
1981 CALL GPDIGI(chuset,chudet)
1986 IMAT 'Material number' I
1987 NAMATE 'Material name' C
1991 RADL 'Radiation lenght' R
1992 ABSL 'Absorption lenght' R
1996 CALL GSMATE(imat,namate,a,z,dens,radl,absl,ubuf,nwbuf)
2001 IMAT 'Material number' I
2002 NAMATE 'Material name' C
2006 NLMAT 'Flag for WMAT' I
2007 WMAT 'Relative weights or n. of atoms in molecule' R
2009 CALL GSMIXT(imat,namate,a,z,dens,nlmat,wmat)
2014 NUMB 'Material number' I
2021 IMATE 'Material number' I
2022 IPART 'Particle number' I
2025 CALL GPRMAT(imate,ipart,mecan,nekbin,elow)
2030 IMATE 'Material number' I
2031 IPART 'Particle number' I
2034 IDM 'ID mode option' I D=0
2036 CALL GPLMAT(imate,ipart,mecan,nekbin,elow,idm)
2037 IDM convention for histogramming mode :
2038 IDM.gt.0 fill, print, keep histogram(s)
2039 IDM.eq.0 fill, print, delete histogram(s)
2040 IDM.lt.0 fill, noprint, keep histogram(s)
2041 If MECAN = 'ALL' all the mechanisms are histogrammed. If the material number
2042 is negative, the cross sections relative to material ABS(IMATE) will
2043 be histogrammed in barns rather than in 1/cm.
2048 IMATE 'Material number' I
2049 IPART 'Particle number' I
2051 MECAN 'List of mechanism' C D='ALL'
2053 CALL GDRMAT(imate,ipart,mecan,nmec)
2054 If MECAN = 'ALL' all the mechanisms are plotted. If the material number
2055 is negative, the cross sections relative to material ABS(IMATE) will
2056 be plotted in barns rather than in 1/cm.
2057 Note that it is not possible to plot anything if GSTMED has not been called
2058 for the material number IMATE.
2063 ITMED 'Medium number' I
2064 CHPAR 'Cut or mechanism' C
2067 CALL GSTPAR(itmed,chpar,parval)
2072 IPART 'Particle number' I
2073 NAPART 'Particle name' C
2080 BRATIO 'Branching ratios' R
2083 CALL GSPART(ipart,napart,itrtyp,amass,charge,tlife,ubuf,nwbuf);
2084 CALL GSDK(ipart,bratio,mode)
2089 NUMB 'Particle number' I
2096 NUMB 'Track number' I
2104 IDEB 'Debug option' C D='ON' R='ON,OFF'
2106 If ideb='ON ' then :
2107 idebug=1, idemin=1, idemax=1000000, itime=1
2109 idebug=0, idemin=0, idemax=0
2119 NUMBER 'Bank number' I D=1
2121 Print a survey of the structure identified by NAME, NUMBER.
2128 NUMBER 'Bank number' I D=1
2129 CHOPT 'Options' C D='BSV'
2131 Display the contents of a bank or a data structure
2132 identified by its NAME and NUMBER.
2133 The output format of the data part is controlled by the internal
2134 or external I/O characteristic.
2135 CHOPT='B' Print the bank.
2136 CHOPT='S' Print the bank contents from left to right Sideways
2137 with up to ten elements per line.
2138 CHOPT='V' Print the vertical (down) structure.
2139 CHOPT='D' Print the bank contents from top to bottom Downwards
2140 with five elements per line.
2141 CHOPT='L' Print the linear structure.
2142 CHOPT='Z' Print the data part of each bank in hexadecimal format
2148 IDIV 'Division number ' I D=2 R=0:24
2149 CHOPT 'Options' C D='M'
2151 Snap of one or more divisions.
2152 Provides a snapshot of one or more divisions in a ZEBRA store.
2153 The kind of information provided is controlled by CHOPT.
2154 CHOPT='M' Print Map entry for each bank
2155 CHOPT='E' Extend map entry to dump all links of each bank
2156 (otherwise only as many links as will fit on a line)
2157 CHOPT='F' Full. Dump all active banks, links and data
2158 CHOPT='K' Kill. Dropped banks to be treated as active
2159 (dropped banks are not normally dumped under D or F option)
2160 CHOPT='L' Dump all Link areas associated with the store
2161 CHOPT='W' Dump the Working space, links and data
2162 CHOPT='Z' Dump the information in hexadecimal.
2168 IDIV 'Division number ' I D=0 R=0:24
2169 CHOPT 'Options' C D='CLSU'
2171 Check the structure of one or more ZEBRA divisions.
2172 The verification detail depends on the settings in CHOPT.
2173 CHOPT='C' Check chaining of banks only
2174 CHOPT='L' Check validity of the structural links (implies 'C')
2175 CHOPT='S' Check the store parameters
2176 CHOPT='U' Check the validity of the up and origin (implies 'C')
2177 CHOPT='F' Errors are considered fatal and generate a call to ZFATAL
2183 IXSTOR 'Store number' I D=0 R=0:24
2185 Display the structure of the ZEBRA store IXSTOR.
2186 Output the parameters characterizing the store, followed by a
2187 list of all divisions and all link areas associated with the store in
2194 IDIV 'Division number' I D=2
2195 PATH 'Name of the doc file' C D=' '
2197 Facility to display the layout of stores and divisions.
2199 CALL DZDDIV(idiv,LDUMMY,path,'IN',1,0,1,IWTYPE)
2205 BANK 'Name of the bank' C
2207 PATH 'Name of the doc file' C D=' '
2208 NUMBER 'Number of the bank' I D=1
2210 Interactive bank display tool.
2212 CALL DZDISP(IXSTOR,LBANK,path,'N',1,0,1,IWTYPE)
2219 PATH 'Name of the RZ directory to analyse' C
2221 Facility to display RZ directory trees.
2223 CALL DZDIRZ(0,LDUMMY,0,path,'N',1,0,1)
2234 LUN 'Fortran unit of the FZ file' I
2235 KEYSU 'Name of the data structure to be retrieved' C
2237 IDENT 'Version of the data structure to be retrieved' I D=0
2239 Equivalent to a call to:
2241 CALL GFIN(LUN,KEYSU,1,IDENT,' ',IER)
2247 LUN 'Fortran unit with which to open the file' I
2248 FILE 'Name of the file to be opened' C
2249 LUNTYP 'Type of FZ file to be opened by GOPEN' C D='XI'
2250 LEN 'Recordlenght of the file' I D=0
2252 CHOPT 'Optional parameter to specify the action' C D=' '
2254 Equivalent to a call to:
2256 CALL GOPEN(LUN,FILE,LUNTYP,LEN,IER)
2258 If CHOPT = I then a call to GFIN or GFOUT will be performed in addition
2259 according to the value of LUNTYP, with the key INIT to save or retrieve
2260 the whole initialization data structure.
2265 LUN 'Fortran unit of the FZ file' I
2266 KEYSU 'Name of the data structure to be saved' C
2268 IDENT 'Version of the data structure to be saved' I D=1
2270 Equivalent to a call to:
2272 CALL GFOUT(LUN,KEYSU,1,IDENT,' ',IER)
2278 LUN 'Fortran unit of the FZ to close' I
2280 Equivalent to a call to:
2282 CALL GCLOSE(LUN,IER)
2294 IQ1 'Lower limit for IQ index' I D=1
2295 IQ2 'Upper limit for IQ index' I D=20
2297 Print the array IQUEST in /QUEST/.
2302 LUN 'Logical unit number' I
2305 CHOPT 'Options' C D=' ' R=' ,A,N,U'
2308 CHOPT=' ' readonly mode
2309 CHOPT='U' update mode
2310 CHOPT='N' create new file
2311 CHOPT='I' Read all structures from existing file
2312 CHOPT='O' Write all structures on file
2317 LUNRZ 'Logical unit number' I
2319 Close an RZ file opened by GRFILE on logical unit LUNRZ.
2326 CHDIR 'Directory name' C
2328 CHOPT 'Options' C D=' '
2330 To create a new RZ directory below the current directory.
2333 RZTAGS(2)='Idvers-NR '
2339 CHPATH 'Path name' C D=' '
2340 CHOPT 'CHOPT' C D=' '
2342 Change or print the current directory.
2343 Ex. CD dir1 ; make DIR1 the new CWD
2344 CD //file1/dir2 ; make //FILE1/DIR2 the new CWD
2345 CD ; print the name of the CWD
2350 OBJECT 'Structure name' C
2352 IDVERS 'Version number' I D=1
2353 CHOPT 'Option' C D=' '
2355 Read data structure identified by OBJECT,IDVERS into memory.
2356 MATE read JMATE structure
2357 TMED read JTMED structure
2358 VOLU read JVOLUM structure
2359 ROTM read JROTM structure
2360 SETS read JSET structure
2361 PART read JPART structure
2362 SCAN read LSCAN structure
2363 INIT read all above data structures
2368 OBJECT 'Structure name' C
2370 IDVERS 'Version number' I D=1
2371 CHOPT 'Option' C D=' '
2373 Write data structure identified by OBJECT,IDVERS to RZ file.
2374 MATE write JMATE structure
2375 TMED write JTMED structure
2376 VOLU write JVOLUM structure
2377 ROTM write JROTM structure
2378 SETS write JSET structure
2379 PART write JPART structure
2380 SCAN write LSCAN structure
2381 INIT write all above data structures
2387 CHPATH 'Path name' C D=' '
2388 CHOPT 'CHOPT' C D=' '
2390 List the contents of a directory (memory or disk).
2391 To list all RZ files currently open, type 'LD //'.
2397 NKEEP 'Number of cycles to keep' I D=1
2399 Purge an RZ directory.
2404 OBJECT 'Structure name' C
2406 IDVERS 'Version number' I D=1
2408 Delete entry identified by OBJECT,IDVERS on RZ file.
2409 OBJECT may be : MATE,TMED,VOLU,ROTM,SETS,PART,SCAN, *
2410 If OBJECT= * delete all entries with IDVERS.
2415 CHDIR 'Lock identifier' C D='RZFILE'
2417 Lock an RZ directory.
2422 CHDIR 'Lock identifier' C D='RZFILE'
2424 Free an RZ directory.
2430 To define parameters for the SCAN geometry. If the routine GUSTEP
2431 and GUKINE are properly instrumented (see examples in GEANX),
2432 when the TRI command is entered NTETA Geantinos will be
2433 tracked through the real detector starting at the vertex position
2434 defined by the command vertex. A simplified version of the geometry
2435 is automatically generated in (ETA,PHI) or (THETA,PHI) following
2436 the option given in the command TETA. The data structure LSCAN
2437 generated may be saved on an RZ file for subsequent processing.
2438 This data structure may be used for fast parametrization techniques.
2442 NPHI 'Number of PHI divisions' I D=90
2444 PHIMIN 'Minimum PHI in degrees' R D=0. R=0.:360.
2445 PHIMAX 'Maximum PHI in degrees' R D=360. R=0.:360.
2447 To specify number of divisions along PHI. If no parameter is
2448 given, the current values of the parameters are displayed.
2453 NTETA 'Number of TETA divisions' I D=90
2455 TETMIN 'Minimum value of TETA' R
2456 TETMAX 'Maximum value of TETA' R
2457 DIVTYP 'Type of TETA division' I R=1:3
2459 To specify number of divisions along TETA.
2460 If DIVTYP=1 divisions in pseudo-rapidity ETA.
2461 If DIVTYP=2 divisions in degrees following the THETA angle.
2462 If DIVTYP=3 divisions in cos(TETA).
2463 If no parameter is given, the current values of the parameters
2469 LIST 'List of master volumes' C
2471 Only boundary crossings of volumes given in LIST will be seen
2472 in the SCAN geometry. If no parameters are given, the current
2473 SCAN volumes will be listed. If a full stop (.) is given, the list
2474 of scan volumes will be erased.
2479 VX 'Scan X-origin' R D=0.
2480 VY 'Scan Y-origin' R D=0.
2481 VZ 'Scan Z-origin' R D=0.
2483 All Geantinos tracked will start from position VX,VY,VZ.
2488 FACTX0 'Scale factor for SX0' R D=100.
2489 FACTL 'Scale factor for SL' R D=1000.
2490 FACTR 'Scale factor for R' R D=100.
2492 Set scale factors for SX0,SL and R. The given scale factors must be
2494 SX0*FACTX0 < 2**15-1 (32767)
2495 SL*FACTL < 2**10-1 (1023)
2496 SR*FACTR < 2**17-1 (131071)
2501 CHOPT 'SCAN mode setting' C R='ON,OFF,INIT'
2503 Switch on/off SCAN mode. If SCAN mode is on, SCAN geantinos
2504 are generated and tracked to fill (or complete) the current
2505 scan data structure. If SCAN mode is off, normal kinematics
2506 generation and tracking will take place. If INIT is given,
2507 the current SCAN data structure (if any) will be dropped
2508 and SCAN mode will be turned on.
2514 IPARAM 'Parametrization Flag' I R=0:1
2515 PCUTGA 'Parametrization Cut for gammas' R
2516 PCUTEL 'Parametrization Cut for electrons' R
2517 PCUTHA 'Parametrization Cut for charged hadrons' R
2518 PCUTNE 'Parametrization Cut for neutral hadrons' R
2519 PCUTMU 'Parametrization Cut for muons' R
2521 Control parametrization at tracking time.
2523 IPARAM=0 No parametrization is performed
2524 IPARAM=1 Parametrization is performed
2526 If parametrization is active and a particle falls below its
2527 parametrization cut, then the particle will be replaced by
2528 a parametrized shower which will be tracked in the SCAN
2534 ID 'Lego plot identifier' I D=2000
2536 VOLUME 'Volume name' C D='XXXX'
2537 CHOPT 'List of options' C D='OPX' R=' ,O,P,I,X,L'
2539 Generates and plot a table of physics quantities such as
2540 the total number of radiation lengths or interaction lengths
2541 in function of the SCAN parameters TETA,PHI.
2542 CHOPT='O' table is generated at Exit of VOLUME.
2543 CHOPT='I' table is generated at Entry of VOLUME.
2544 CHOPT='X' radiation lengths
2545 CHOPT='L' Interaction lengths
2546 CHOPT='P' Plot the table
2547 If VOLUME='XXXX' Mother volume is used.
2552 IDPHI 'Histogram/phi identifier' I D=1000
2554 VOLUME 'Volume name' C D='XXXX'
2555 CHOPT 'List of options' C D='OPX' R=' ,O,P,I,X,L'
2557 Generates and plots an histogram of physics quantities such as
2558 the total number of radiation lengths or interaction lengths
2559 as a function of the SCAN parameter TETA for a given value of PHI.
2560 CHOPT='O' histogram is generated at Exit of VOLUME.
2561 CHOPT='I' histogram is generated at Entry of VOLUME.
2562 CHOPT='X' radiation lengths
2563 CHOPT='L' Interaction lengths
2564 CHOPT='P' Plot the histogram
2565 If VOLUME='XXXX' Mother volume is used.
2566 The histogram identifier IDPHI is used to also identify which
2567 PHI division to plot: IPHI=MOD(IDPHI,1000).
2568 If IPHI=0, then all PHI divisions are generated (not plotted)
2569 with histogram identifiers IDPHI+PHI division number.
2573 >Menu /GEANT/PHYSICS
2575 Commands to set physics parameters.
2580 IANNI 'Flag IANNI' I D=1 R=0,1,2
2582 To control positron annihilation.
2583 IANNI=0 no annihilation
2584 =1 annihilation. Decays processed.
2585 =2 annihilation. No decay products stored.
2591 IAUTO 'Flag IAUTO' I D=1 R=0,1
2593 To control automatic calculation of tracking medium parameters:
2594 IAUTO=0 no automatic calculation;
2595 =1 automati calculation.
2601 IBREM 'Flag IBREM' I D=1 R=0,1,2
2603 To control bremstrahlung.
2604 IBREM=0 no bremstrahlung
2605 =1 bremstrahlung. Photon processed.
2606 =2 bremstrahlung. No photon stored.
2612 ICKOV 'Flag ICKOV' I D=0 R=0,1,2
2614 To control Cerenkov production
2615 ICOMP=0 no Cerenkov;
2617 =2 Cerenkov with primary stopped at each step.
2623 ICOMP 'Flag ICOMP' I D=1 R=0,1,2
2625 To control Compton scattering
2627 =1 Compton. Electron processed.
2628 =2 Compton. No electron stored.
2634 IDCAY 'Flag IDCAY' I D=1 R=0,1,2
2636 To control Decay mechanism.
2638 =1 Decays. secondaries processed.
2639 =2 Decays. No secondaries stored.
2645 IDRAY 'Flag IDRAY' I D=1 R=0,1,2
2647 To control delta rays mechanism.
2648 IDRAY=0 no delta rays.
2649 =1 Delta rays. secondaries processed.
2650 =2 Delta rays. No secondaries stored.
2656 EKMIN 'Minimum energy of the tables' R D=1E-5
2657 EKMAX 'Maximum energy of the tables' R D=1E+4
2658 NEKBIN 'Number of bins in the tables' I D=90 R=1:200
2660 To define the range and binning of internal tables.
2666 IHADR 'Flag IHADR' I D=1
2668 To control hadronic interactions.
2669 IHADR=0 no hadronic interactions.
2670 =1 Hadronic interactions. secondaries processed.
2671 =2 Hadronic interactions. No secondaries stored.
2677 LABS 'Flag LABS' I D=0
2679 To control absorbtion of Cerenkov photons:
2680 LABS=0 no absorbtion of photons;
2681 LABS=1 absorbtion of photons;
2687 ILOSS 'Flag ILOSS' I D=2 R=0,1,2,3,4
2689 To control energy loss.
2690 ILOSS=0 no energy loss;
2691 =1 restricted energy loss fluctuations;
2692 =2 complete energy loss fluctuations;
2694 =4 no energy loss fluctuations.
2695 If the value ILOSS is changed, then cross-sections and energy loss
2696 tables must be recomputed via the command 'PHYSI'.
2702 IMULS 'Flag IMULS' I D=1 R=0,1,2,3
2704 To control multiple scattering.
2705 IMULS=0 no multiple scattering.
2706 =1 Moliere or Coulomb scattering.
2707 =2 Moliere or Coulomb scattering.
2708 =3 Gaussian scattering.
2714 IMUNU 'Flag IMUNU' I D=1 R=0,1,2
2716 To control muon nuclear interactions.
2717 IMUNU=0 no muon-nuclear interactions.
2718 =1 Nuclear interactions. Secondaries processed.
2719 =2 Nuclear interactions. Secondaries not processed.
2725 IPAIR 'Flag IPAIR' I D=1 R=0,1,2
2727 To control pair production mechanism.
2728 IPAIR=0 no pair production.
2729 =1 Pair production. secondaries processed.
2730 =2 Pair production. No secondaries stored.
2736 IPFIS 'Flag IPFIS' I D=1 R=0,1,2
2738 To control photo fission mechanism.
2739 IPFIS=0 no photo fission.
2740 =1 Photo fission. secondaries processed.
2741 =2 Photo fission. No secondaries stored.
2747 IPHOT 'Flag IPHOT' I D=1 R=0,1,2
2749 To control Photo effect.
2750 IPHOT=0 no photo electric effect.
2751 =1 Photo effect. Electron processed.
2752 =2 Photo effect. No electron stored.
2758 IRAYL 'Flag IRAYL' I D=1 R=0,1
2760 To control Rayleigh scattering.
2761 IRAYL=0 no Rayleigh scattering.
2768 ISTRA 'Flag ISTRA' I D=0 R=0,1,2
2770 To control energy loss fluctuation model:
2771 ISTRA=0 Urban model;
2773 =2 PAI+ASHO model (not active at the moment).
2779 ISYNC 'Flag ISYNC' I D=1 R=0,1
2781 To control synchrotron radiation:
2782 ISYNC=0 no synchrotron radiation;
2783 =1 synchrotron radiation.
2789 CUTGAM 'Cut for gammas' R D=0.001
2790 CUTELE 'Cut for electrons' R D=0.001
2791 CUTHAD 'Cut for charged hadrons' R D=0.01
2792 CUTNEU 'Cut for neutral hadrons' R D=0.01
2793 CUTMUO 'Cut for muons' R D=0.01
2794 BCUTE 'Cut for electron brems.' R D=-1.
2795 BCUTM 'Cut for muon brems.' R D=-1.
2796 DCUTE 'Cut for electron delta-rays' R D=-1.
2797 DCUTM 'Cut for muon delta-rays' R D=-1.
2798 PPCUTM 'Cut for e+e- pairs by muons' R D=0.01
2799 TOFMAX 'Time of flight cut' R D=1.E+10
2800 GCUTS '5 user words' R D=0.
2802 To change physics cuts. If no parameter is given, the list
2803 of the current cuts is printed.
2804 If the default values (-1.) for BCUTE ,BCUTM ,DCUTE ,DCUTM
2805 are not modified, they will be set to CUTGAM,CUTGAM,CUTELE,CUTELE
2807 If one of the parameters from CUTGAM to PPCUTM included
2808 is modified, cross-sections and energy loss tables must be
2809 recomputed via the command 'PHYSI'.
2814 IPART 'GEANT particle number' I
2815 IMATE 'GEANT material number' I
2816 STEP 'step length in centimeters' R
2818 NPOINT 'number of logarithmically spaced energy points' I D=10 R=2:100
2820 This routine prints the relevant parameters linked with the energy loss
2826 Call the GEANT initialisation routine GPHYSI to recompute
2827 the tables of cross-sections and energy loss. This command
2828 must be invoked after CUTS, LOSS or ERAN commands.
2838 The routines in the file FNAME will be compiled by COMIS.
2839 If routines with names: UGEOM,GUKINE,GUOUT,UGLAST are found,
2840 then they will be automatically called by GXINT instead of
2841 the routines with the same names compiled with the standard
2842 Fortran compiler and linked with the application.
2843 The user callable routines from the GEANT library as well as
2844 routines from PACKLIB (HBOOK,HPLOT,HIGZ,ZEBRA) may be called
2845 from these user routines. All GEANT common blocks may be
2847 In case where the routine UGEOM is called several times,
2848 it is important to DROP all the initialisation data structures
2849 JVOLUM,JMATE,JTMED,etc already in memory by using the routine GIDROP.
2850 Example of an interactive session where the routine UGEOM is modified:
2852 GEANT > Edit ugeom.for
2853 GEANT > Fortran ugeom.for
2857 GEANT > Edit ugeom.for
2858 GEANT > Fortran ugeom.for
2863 If FNAME='-', calls to user routines is reset and standard
2864 routines called instead.