GEANT321/doc/cmzt/Versions/v3_21

   1 *
   2 * $Id$
   3 *
   4 * $Log$
   5 * Revision 1.1.1.1  1995/10/24 10:19:48  cernlib
   6 * Geant
   7 *
   8 *
   9 *CMZ :  3.21/04 21/03/95  17.58.14  by  S.Giani
  10 *-- Author :
  11 *
  12 *::>          VERSION  3.21/04   21/03/95 17.58.14
  13 *
  14 *
  15 *::>          VERSION  3.21/03   17/10/94 16.00.29
  16 *
  17 *
  18 *::>          VERSION  3.21/02   25/07/94 15.24.16
  19 *
  20 *
  21 *::>          VERSION  3.21/01   28/03/94 16.10.10
  22 *
  23 *
  24 *::>          VERSION  3.21/00   28/02/94 16.32.43
  25 *
  26 **********************************************************************
  27 *                                                                    *
  28 *                       G E A N T   3 . 2 1                          *
  29 *                                                                    *
  30 *                            C E R N                                 *
  31 *                                                                    *
  32 *                       15/3/94  Geneva  CH                          *
  33 *                                                                    *
  34 **********************************************************************
  35 *
  36 * GEANT 3.21 has a new tracking logic to speed up the search over the
  37 *            volumes positioned inside a given node (Virtual Divisions).
  38 *            The new tracking is initialized by GGCLOS in a completely
  39 *            transparent way and in a negligible amount of time.
  40 *            The GUNEAR and GSORD techniques are automatically disabled
  41 *            and no change in the user code is needed. Anyway, it is
  42 *            still possible to get the old Geant tracking recompiling
  43 *            GEANT and GEANG with the patchy flag OLD.
  44 *            To implement the new logic, the following routines have
  45 *            written or deeply modified: GTNEXT, GGCLOS, GTMEDI, GNEXT,
  46 *            GMEDIA, GINVOL, GVDCAR, GVDRAD, GVDPHI, GVDLIM, GSCVOL.
  47 *                                                               S.Giani.
  48 *
  49 * GEANT 3.21 has a new logic for MANY volumes: the tracking is not
  50 *            based anymore on a step-search, but relies on a algorithm
  51 *            similar to the one used for ONLY volumes (giving then
  52 *            basically the same speed obtained for ONLY). The new logic
  53 *            of MANY allows boolean operations between volumes, clipping
  54 *            of protuding objects, divisions along arbitrary axis,
  55 *            overlaps of 1-D structures to build 2-D or 3-D structures
  56 *            (like honeycombs, spaghetti, etc.). To implement the new
  57 *            logic, the following routines have been written or modified:
  58 *            GTNEXT, GTMEDI, GNEXT, GMEDIA, GTELEC, GTCKOV, GTHADR,
  59 *            GTHION, GTMUON.
  60 *                                                             S.Giani.
  61 *
  62 * GEANT 3.21 contains modifications to GTELEC, GTGAMA, GDRAY in order
  63 *            to exploit the e.m. low energy effects introduced in 3.16
  64 *            without any overhead in performance compared with 3.15.
  65 *            A bug in the delta ray productions is now fixed.
  66 *                                                        Thanks to L.Urban.
  67 *
  68 * GEANT 3.21 provides a ray-tracing as a new visualization tool. The
  69 *            ray-tracing is based on the new Geant tracking and allows
  70 *            the visualization of the boolean operations results.
  71 *            The following routines have been written or modified:
  72 *            GDRACK, GDSTEP, GDNINO, GDOPT, GDRAYT, GDRAW, GDINIT,
  73 *            GTRIG, GTRIGI, GDKINE, GDLUMI, GLIGHT, GFIND, GXDRAW,
  74 *            CDF.
  75 *                                                             S.Giani.
  76 *
  77 * A new GGPERP has been written to cope with almost all the shapes.
  78 *            It is now of general utility both for the Cerenkovs
  79 *            and for the ray-tracing.
  80 *                                               Thanks to R.Jones.
  81 *
  82 * A new set of routines to compute distances to surfaces has been
  83 *            introduced (GVSAFE, GVGRAD, GVPSQR, GVSIGN, GVLONG).
  84 *                                               Thanks to V.Perev.
  85 *
  86 * Important fix in GPHYSI: user words for media were overwriting
  87 *            ILABS, ISYNC and ISTRA.
  88 *                                    Thanks to S.Giani, S.Ravndal.
  89 *
  90 * GXOBJ has been improved to handle DETs as objects of a SET directory
  91 *            in the HITS class.
  92 *                                    Thanks to S.Giani, A.Dellacqua.
  93 *
  94 * GTNINO has been corrected in order to print the mechanism NEXT when
  95 *            crossing boundaries.
  96 *                                    Thanks to S.Giani, M.Maire.
  97 *
  98 * GFMFIN and GFMDIS are the new names for the 316 routines FMUFIN and
  99 *        FMDIST to avoid clashes with the program library. GFMFIN has
 100 *        also been corrected for a 316 bug affecting neutrons.
 101 *                             Thanks to C.Zeitniz,K.Lassila,S.Ravndal.
 102 *
 103 * GINTRP and GINPCO have been corrected to handle sharp edges with no
 104 *        risk of float exception.
 105 *                                Thanks to L.Gibbons.
 106 *
 107 * GGPERP has been corrected in reading the TRAP parameters and a systematic
 108 *        consistency of single precision or double precision variables in
 109 *        the MAX,MIN,MOD operands is provided.
 110 *                                             Thanks to S.Ravndal.
 111 *
 112 * FEREVV and FERHAV are now new from FLUKA.
 113 *                                             Thanks to A.Ferrari.
 114 *
 115 * FEREVV has been adapted to the Geant random numbers generators.
 116 *                                             Thanks to S.Ravndal.
 117 *
 118 * In GRNDMQ the option `SH' is provided to write the seeds ISEED1 and
 119 *         ISEED2 to the event header bank.
 120 *                                             Thanks to S.Ravndal.
 121 *
 122 * In GMORIN,FLINIT and GHEINI the interface messages are removed.
 123 *                                             Thanks to S.Ravndal.
 124 *
 125 * In GGDPAR now PAR(3) is correctly initialized to 1.for slices of
 126 *         a PGON.
 127 *
 128 *
 129 * The seq GCNUM has been added to the deck ERTRGO in Geane.
 130 *
 131 *
 132 *
 133 * G321X and G321M are the new names for the KUIP output of the CDF.
 134 *                                             Thanks to S.Giani.
 135 *
 136 * GMORIN has been modified in order to read correctly xsneut.dat
 137 *          for MSDOS and WINNT as well.
 138 *                                         Thanks to V.Fine.
 139 *
 140 * GEAMON has been modified in order to print Real and CPU Time longer
 141 *          than 999 sec.
 142 *                                         Thanks to R.Brun.
 143 *
 144 *
 145 **********************************************************************
 146 *                                                                    *
 147 *                       G E A N T   3 . 2 1                          *
 148 *                                                                    *
 149 *                      Correction cradle 0.01                        *
 150 *                                                                    *
 151 *                            C E R N                                 *
 152 *                                                                    *
 153 *                       29/3/94  Geneva  CH                          *
 154 *                                                                    *
 155 **********************************************************************
 156 *
 157 * GZINIT and GSCKOV have been corrected (compared with 3.16) in order
 158 *        to initialize correctly the link area for the Zebra pointers
 159 *        of the banks relative to straggling, Cerenkovs, Asho, etc.
 160 *        Four pointers were not seen by Mzlint and, moreover, the link
 161 *        area (containing pointers for straggling, etc.) was initialized
 162 *        only if using Cerenkovs.
 163 *                                   Thanks to M.Verzocchi and R.Brun.
 164 *
 165 * FERHAV in Geanf was provoking a float exception due to a negative sqrt.
 166 *        The problem was fixed thanks to a comparison with the original
 167 *        Fluka code.
 168 *                                   Thanks to B.Slater and S.Giani.
 169 *
 170 * PREPRE has been deoptimized on HP.
 171 *                                   Thanks to A.Kulik.
 172 *
 173 * In GSKINE a pointer was not saved (in the previous versions as well).
 174 *                                   Thanks to A.Dellacqua.
 175 *
 176 * In GTNEXT and GNEXT precision problems are handled more effectively.
 177 *                                   Thanks to P.Ward, C.Markou and S.Giani.
 178 *
 179 * In GTNEXT the blocking of the volume from which particles are exiting
 180 *        has been corrected.
 181 *                                   Thanks to P.Ward.
 182 *
 183 * Problems in the interactive materials, mixtures and media definition
 184 *        have been solved. (Routines GXCONT,GXGEOM). Idem for the 3D
 185 *        Dspec (routine GXDRAW).
 186 *
 187 * In order to handle correctly calls to Gfcoor in phi (for x=0 or y=0)
 188 *        GTMEDI, GMEDIA and GINVOL have been modified.
 189 *                                    Thanks to P.Ward.
 190 *
 191 * In the ray-tracing, GDNINO has been corrected for precision problems.
 192 *                                    Thanks to S.Giani.
 193 *
 194 *
 195 **********************************************************************
 196 *                                                                    *
 197 *                       G E A N T   3 . 2 1                          *
 198 *                                                                    *
 199 *                      Correction cradle 0.02                        *
 200 *                                                                    *
 201 *                            C E R N                                 *
 202 *                                                                    *
 203 *                       25/7/94  Geneva  CH                          *
 204 *                                                                    *
 205 **********************************************************************
 206 *
 207 * In GTNEXT precision problems in phi have been solved.
 208 *                                    Thanks to A.Rosowsky and A.Givernaud.
 209 *
 210 * In GDRAY a missing protection for delta ray production was added.
 211 *                                    Thanks to L.Urban.
 212 *
 213 * In GINVOL the variable ingt is corrrectly reset. It was not reset in
 214 *       all the previous versions of GEANT as well.
 215 *                                    Thanks to R.Jones.
 216 *
 217 * GTCKOV, GTHADR, GTHION, GLISUR have been corrected and extended compared
 218 *       with the first implementation in 3.16 for Cerenkov tracking.
 219 *                                    Thanks to R.Jones.
 220 *
 221 * GSKPHO has been fixed for Cerenkov physics compared with the previous
 222 *       version of GEANT (3.16).
 223 *                                    Thanks to S.Bottai.
 224 *
 225 * GFTMAT had a wrong pointer to plot the muon-nucleus cross sections since
 226 *       version 3.16. Now it is fixed.
 227 *                                    Thanks to C.Wiebusch and S.Ravndal.
 228 *
 229 * GFINDS has been restored to the 3.15 functionality.
 230 *                                    Thanks to A.Dellacqua and S.Tonse.
 231 *
 232 * GFANG has been fixed for a special pathological case.
 233 *                                    Thanks to C.Zosi.
 234 *
 235 * GSROTM  has been fixed for a special pathological case.
 236 *                                    Thanks to A.Vlassov.
 237 *
 238 * GDRAYT does not increment ievent anymore.
 239 *                                    Thanks to P.Jacobs.
 240 *
 241 * GLVOLU has been fixed to recompute the pointers of the geometrical tree
 242 *       even if the names and numbers were identical.
 243 *                                    Thanks to R.Jones, T.Noble, S.Ravndal.
 244 *
 245 * GVDRAD has been fixed to compute the scope of tubs.
 246 *                                    Thanks to K.Albiol and S.Giani.
 247 *
 248 * GDSPEC has been fixed to print degrees for the 3-rd parameter of sphe.
 249 *                                    Thanks to S.Ravndal.
 250 *
 251 * GVDPHI has been fixed for a special pathological case (only for volumes
 252 *       with planar surfaces).
 253 *                                    Thanks to M.Kaspar and S.Giani.
 254 *
 255 * GXCONT is changed to handle correctly the variables ievent and nevent.
 256 *                                    Thanks to P.Jacobs.
 257 *
 258 * SVGAM7 in Gheisha has now a protection against the log(0).
 259 *                                    Thanks to D.Heck.
 260 *
 261 * FEREVV and FERHAV have been updated.
 262 *                                    Thanks to A.Ferrari.
 263 *
 264 * In FLUFIN the correct mechanism is given to lmec in case of energy
 265 *       below cut-off.
 266 *                                    Thanks to B.Slater and S.Ravndal.
 267 *
 268 * In NUDISV the array anpowr is saved.
 269 *                                    Thanks to S.Ravndal.
 270 *
 271 * In GMXSEC imed is used instead of nmed.
 272 *                                   Thanks to K.L.Perini.
 273 *
 274 *
 275 *
 276 **********************************************************************
 277 *                                                                    *
 278 *                       G E A N T   3 . 2 1                          *
 279 *                                                                    *
 280 *                  Correction version 0.03 = 94b release             *
 281 *                                                                    *
 282 *                            C E R N                                 *
 283 *                                                                    *
 284 *                       18/10/94  Geneva  CH                         *
 285 *                                                                    *
 286 **********************************************************************
 287 *
 288 * In GDCGVW traslation problems have been fixed (graphics).
 289 *                                    Thanks to J.Stepanek, S.Giani.
 290 *
 291 * In GLTRAC nlevel is now initialized to 0 and GTNEXT has been improved
 292 *    further for handling precision problems.
 293 *                                    Thanks to A.Dellacqua.
 294 *
 295 * GDRELX and GTELEC have been corrected for energy deposition problems.
 296 *                                    Thanks to L.Urban.
 297 *
 298 * GSDVT2 has been corrected for divisions in phi.
 299 *                                    Thanks to V.Popov.
 300 *
 301 * GPVOLU and GFMOTH are now able to handle cases where copy number = 1
 302 *    is not defined.
 303 *                                    Thanks to S.Giani, S.Nowak.
 304 *
 305 * GTREVE has been changed in order to avoid the loading of the graphics
 306 *    package if not desired.
 307 *                                    Thanks to S.Giani, I.Bird.
 308 *
 309 * GUPHAD, GUHADR and FLINIT give warnings and conditions for the use of FLUKA.
 310 *                                    Thanks to A.Ferrari.
 311 *
 312 * GVDPHI and GVDCAR have been further improved to compute the scope of pgon.
 313 *                                    Thanks to P.Kent.
 314 *
 315 * GTELEC now computes cross sections for SYNC only if there is mag. field.
 316 *                                    Thanks to M.Maire.
 317 *
 318 * GPCXYZ printing format can now handle bigger dimensions.
 319 *                                    Thanks to S.Ravndal.
 320 *
 321 * GDRMAT and GPRMAT headers are now consistent with GFTMAT. In GPLMAT `STEP'
 322 *    and `RANGE' do not belong anymore to `ALL' as physical processes.
 323 *                                    Thanks to M.Maire, S.Ravndal.
 324 *
 325 * GPHYSI is now setting dray=0 if loss=0. Labs, Sync and Stra can be switched
 326 *    interactively. Decay and Sync are now also defined in vacuum.
 327 *                                    Thanks to M.Maire, S.Ravndal.
 328 *
 329 * In GRNDMQ the ' ' and the 'G' options work now as described in the doc.
 330 *                                    Thanks to S.ONeale, S.Ravndal.
 331 *
 332 * GPART calls GSDK anyhow to give flexibility at run time (for the interactive
 333 *    version). The values relatives to the particles are now updated to
 334 *    Phys.Rev. D50 1173 1994.
 335 *                                    Thanks to S.Ravndal, F.H.Heinsius.
 336 *
 337 * In GDEXCA (Cadint) a bug from 316 is now fixed.
 338 *                                    Thanks to J.Vuoskosky.
 339 *
 340 * GBREME and GSYNC (see the new doc for the meaning of sync=0,1,2,3) have
 341 *    been improved and the two relative mechanism can happen at the same time.
 342 *                                    Thanks to M.Maire.
 343 *
 344 *
 345 *
 346 **********************************************************************
 347 *                                                                    *
 348 *                       G E A N T   3 . 2 1                          *
 349 *                                                                    *
 350 *                  Correction version 0.04 = 95a release             *
 351 *                                                                    *
 352 *                            C E R N                                 *
 353 *                                                                    *
 354 *                       28/3/95  Geneva  CH                          *
 355 *                                                                    *
 356 **********************************************************************
 357 *
 358 * In GTCKOV the variable RIN1 is not reset to 0 any longer.
 359 *                                           Thanks to A.Delpozo
 360 *
 361 * In GSTRAC heavy ions have now the correct itrtyp=8.
 362 *                                           Thanks to Yfwang
 363 *
 364 * In GTNEXT the variable SNXT and SNXNEW are compared with BIG1 to cure
 365 *  precision problems.
 366 *                                           Thanks to A.Dellacqua, P.Ward.
 367 *
 368 * In GPVOLU and GPVOLX the variable MOTHER is now declared as CHARACTER*4.
 369 *                                           Thanks to V.Tretyak
 370 *
 371 * In GGCLOS a test on IQUEST(11).le.0 has been changed into IQUEST(11).lt.0.
 372 *                                           Thanks to P.Ward
 373 *
 374 * In GLTRAC nlevel is reset to 0 only if GONLY(NLEVEL).EQ.0.
 375 *                                           Thanks to P.Ward, S.Giani
 376 *
 377 * In GDCGOB a bug on the unit for the angles of spheres has been fixed.
 378 *                                           Thanks to J.Allison
 379 *
 380 * GPRMAT is now protected against floating point exceptions.
 381 *                                           Thanks to R.Ruber
 382 *
 383 * GRNDMQ has been improved storing the seeds.
 384 *                                           Thanks to S.Ravndal.
 385 *
 386 * In GTNINO the time of flight is now updated.
 387 *                                           Thanks to S.Ravndal.
 388 *
 389 * In GXPHYS, GXCS, the interactive facilities have been extended.
 390 *                                           Thanks to M.Maire
 391 *
 392 * GPCXYZ has a corrected format to print.
 393 *                                           Thanks to S.Ravndal
 394 *
 395 * A GUXCS dummy routine has been created to extend the interactive use of COMIS.
 396 *                                           Thanks to M.Maire.
 397 *
 398 * GDPRTR is now able to recognize b/w Xterminals before drawing the tree.
 399 *                                           Thanks to S.Giani, B.Lorazo.
 400 *
 401 * GLVOLU has been corrected in the computation of GTRAN for the slices of
 402 *  divided volumes.
 403 *                                           Thanks to S.Giani.
 404 *
 405 * GXGEOM now allows correctly the definition of the basic shapes interactively.
 406 *                                           Thanks to S.Giani.
 407 *
 408 * GETVER and GINIT are improved to extract and print automatically the version.
 409 *                                           Thanks to M.Maire.
 410 *
 411 * GFTMAT and GPLMAT have been corrected in plotting the cross sections.
 412 *                                           Thanks to S.Ravndal, G.Zholobov.
 413 *
 414 * GLANDZ has been protected against floating point exceptions.
 415 *                                           Thanks to L.Urban, B.Fox.
 416 *
 417 * GLUDKY contains now the time of interaction (Lund).
 418 *                                           Thanks to S.Ravndal, C.Wiebusch.
 419 *
 420 * GPAIRG has been protected against undefined values for the variable X.
 421 *                                           Thanks to L.Urban, H.Ihssen.
 422 *
 423 * GINIT, GFFGO and GTELEC are modified in order to include a new data card,
 424 *  ABAN, which allows to follow the electrons down to the cuts value even
 425 *  when they are not doing physics (ABAN=0) or to propagate them up to the
 426 *  end-point with an optimized algorithm (ABAN=2). The default is ABAN=1
 427 *  (normal 3.21 behaviour).
 428 *                                           Thanks to L.Urban, S.Ravndal.
 429 *
 430 * BRTTP has been protected in case the file flukaaf.dat is not found.
 431 *                                           Thanks to S.Ravndal.
 432 *
 433 * GMORIN has been modified and GMISET has been created in order to interface
 434 *  to the new version of Micap (from Gcalor 1.03.08) which has been introduced
 435 *  in GEANM.
 436 *                                           Thanks to S.Ravndal.
 437 *
 438 * GPRSGG and GBRSGE have been modified to assume constant cross sections
 439 *  above 100 GeV both for pair production and for bremsstrahlung.
 440 *                                           Thanks to G.Battistoni, L.Urban.
 441 *
 442 * GVDPHI has been corrected and improved for the scope computation of tubes,
 443 *  cons and polycons which are rotated and translated at the same time.
 444 *                                           Thanks to S.Giani, J.P.Vialle.
 445