Bugfix in AliPoints2Memory
[u/mrichter/AliRoot.git] / GEANT321 / doc / history / v_315.F
CommitLineData
fe4da5cc 1*
2* $Id$
3*
4* $Log$
5* Revision 1.1.1.1 1995/10/24 10:19:48 cernlib
6* Geant
7*
8*
9#include "sys/CERNLIB_machine.h"
10#include "pilot.h"
11*CMZ : 3.21/02 29/03/94 15.41.18 by S.Giani
12*-- Author :
13*
14*
15* Maintenance and support for the GEANT MonteCarlo system
16* are performed by the Simulation Section of the CERN
17* Computing and Networks division. For any problem please
18* contact:
19*
20* Federico Carminati
21* CERN-CN
22* 1211 Geneva 23
23* Tel: +41 22 767.4959
24* Telefax: +41 22 767.7155
25* E-mail:
26* BITnet/EARN: FCA@CERNVM
27* DECnet : VXCERN::FCA (Node 22.190)
28* Internet : fca@cernvm.cern.ch
29*
30* An electronic GEANT discussion list has been set up on CERNVM.
31* To be included in the mailing list it is enough to send the
32* following message:
33*
34* TELL LISTSERV AT CERNVM SUBSCRIBE LGEANT Firstname Lastname Inst
35*
36* People who do not dispose of a TELL (or equivalent) bitnet
37* command should send a mail to LISTSERV@CERNVM containing the
38* only line:
39*
40* SUBSCRIBE LGEANT Firstname Lastname Inst
41*
42*================================= History
43*
44*===> 20/NOV/1990
45*
46* Current plans set the release of GEANT 3.15 around the third
47* quarter of 1991.
48*
49*===> 27/NOV/1990
50*
51* Correction in GTRACK to set correctly the value of PREC
52* to match machine precision. Thanks to Stephan Egli and
53* M.Maire.
54* Introduced in the 314 correction cradle.
55*
56*===> 28/NOV/1990
57*
58* Correction in GNTUBE to avoid square root of a negative
59* number. Thanks to M.Maire
60* Introduced in the 314 correction cradle.
61*
62* Eliminate TOFG=0 from GINIT, already done in GTRIGI
63* Thanks to F.Carminati
64* Introduced in the 314 correction cradle.
65*
66*===> 29/NOV/1990
67*
68* UNIX and BSLASH flags correctly set inside *GEANT.
69* Thanks to A.Nathaniel
70* GPMATE now prints also the components in case of
71* mixtures. If NUMB is <0 then materials from 1 to
72* -NUMB are printed. Thanks to Stephan Egli.
73*
74* Corrections in GRGET (obsolete), GGET, GRIN to reset
75* correctly the number in GCNUM. Thanks to Stephan Egli
76* Introduced in the 314 correction cradle.
77*
78*===> 30/NOV/1990
79*
80* New version of GNOELT and GNELTU, thanks to Yoshihisa Iga
81* and R.Nierhaus
82* Introduced in the 314 correction cradle.
83*
84* Correction in GTRAK, all mechanisms turned off in the
85* vacuum. Thanks to S.Egli
86* Introduced in the 314 correction cradle.
87*
88*
89*===> 03/DEC/1990
90*
91* New routine GXPICK for the interactive version (Command
92* DRAWING/PICK) to pick a point in a detector and return
93* volume/material characteristics. Thanks to S.Egli.
94*
95*===> 04/DEC/1990
96*
97* MAXNOD from 2000 to 8000 in GDTR0. Suggestion by S.Egli.
98* Introduced in the 314 correction cradle.
99*
100* Bug corrected in GGCLOS, do not touch ISEARC in case
101* of user search (ISEARC.GT.0). Thanks to D.Ward.
102* Introduced in the 314 correction cradle.
103*
104*===> 05/DEC/1990
105*
106* Correction in TWOB, get correctly Tmin and Tmax and
107* their distribution. Thanks to Marko Mikuz.
108* Introduced in the 314 correction cradle.
109*
110*===> 06/DEC/1990
111*
112* Correction in TWOCLU and GENXPT to improve energy
113* conservation. Thanks to S.Egli.
114* Introduced in the 314 correction cradle.
115*
116*===> 11/DEC/1990
117*
118* Generic function LOG introduced in GHFRAK. Generic
119* functions MAX and MIN introduced in several routines
120* in GGEOM (HYPE) and CGPACK. Generic SIN, COS, ATAN and
121* SQRT introduced in CGPACK. Thanks to Federico Carminati.
122* Introduced in the 314 correction cradle.
123*
124*===> 13/DEC/1990
125*
126* GTAU routine to be obsoleted.
127*
128* Supporting link for the LSCAN data structure mispelled
129* in several places. Corrected now. Thanks to Isabel Josa.
130* Introduced in the 314 correction cradle.
131*
132* Unused variable IPOT eliminated from GHFFER.
133*
134* Changes in the interactive part for ATC GKS which does
135* not use FORTRAN unit for metafiles. (routines GXINT and
136* GXGCON). Thanks to W.Koellner
137*
138*===> 14/DEC/1990
139*
140* REAL*8 changed to DOUBLE PRECISION in CGPACK.
141* Life of the D+/- mesons was 100 times too short.
142* Thanks to Cheng He Sheng.
143* Introduced in the 314 correction cradle.
144*
145*===> 19/DEC/1990
146*
147* Modification in GDRAW and GDSPEC to correct a bug.
148* Modifications in GLVOLU not to print a warning message
149* caused by drawing.
150* Introduced in the 314 correction cradle.
151*
152*===> 20/DEC/1990
153*
154* Correction in GRAYLI affecting Rayleigh effect in
155* mixtures.
156* Introduced in the 314 correction cradle.
157*
158* Recode of GRAYLI and GHESIG to be safer with respect
159* to compiler optimization.
160* Thanks to Rafi Yaari.
161*
162*===> 8/JAN/1991
163*
164* Corrections in GRGET, GGET, GRIN to set correctly
165* NPART, NTMED, NMATE, NVOLUM. Thanks to Lee Roberts.
166* Introduced in the 314 correction cradle.
167*
168* New meaning for ISVOL. If <0 the material is not a
169* a detector but the tracking parameters are the same
170* as in the case ISVOL>0 (sensitive medium).
171* Thanks to Michel Maire.
172* Introduced in the 314 correction cradle.
173*
174*===> 9/JAN/1991
175*
176* TMAXFD not recomputed with IGAUTO=1 if the value given
177* by the user >= 0. Default value reset to 20. instead of
178* 60. degrees. Thanks to W.Mueller and M.Maire
179* Introduced in the 314 correction cradle.
180*
181*===> 10/JAN/1991
182*
183* Corrections in GPFIS, GPFISI and GPHYSI for the fotofission
184* and fotoabsorbtion (Giant resonance) of photons on Uranium
185* and plutonium. Thanks to Harm Fesefeld.
186* Introduced in the 314 correction cradle.
187*
188* Correction in GPHYSI in case ILOSS=2 and more than one
189* tracking medium with the same material number then a
190* spurious error message was printed. Thanks to Kevin Sparks
191* Introduced in the 314 correction cradle.
192*
193* Corrections in GENXPT, TWOCLU, TWOB, CINEMA and EXNU.
194* Thanks to Harm Fesefeldt.
195* Introduced in the 314 correction cradle.
196*
197*===> 11/JAN/1991
198*
199* Correction in GPPART in case of user words. Thanks to Walter
200* Mueller.
201* Introduced in the 314 correction cradle.
202*
203* Correction in GSTMED to protect for EPSIL <=0. Thanks to
204* Michel Lefebvre.
205* Introduced in the 314 correction cradle.
206*
207*===> 14/JAN/1991
208*
209* Correction in GRGET to read the banks in the correct
210* division. Thanks to M.Maire.
211* Introduced in the 314 correction cradle.
212*
213* Improvement of GRIN, GROUT. New global keywords KINE and
214* TRIG can now be used. Thanks to Federico Carminati
215* Introduced in the 314 correction cradle.
216*
217*===> 17/JAN/1991
218*
219* Correction in GSVOLU, GSDVN, GSDVN2, GSDVT, GSDVT2 not to
220* use uninitialized LINATT, set by GDINIT. This should take
221* care of the problem of empty drawings. Thanks to Simone Giani.
222* Introduced in the 314 correction cradle.
223*
224* Correction in GPKINE to print correctly user words and vertex
225* information. Thanks to Kevin Sparks.
226* Introduced in the 314 correction cradle.
227*
228* Correction in GPVERT to print correctly user words and vertex
229* information. Thanks to Soren G. Frederiksen.
230* Introduced in the 314 correction cradle.
231*
232*===> 18/JAN/1991
233*
234* Max. GCKING stack size parametrized. The actual value of the
235* stack is not changed (100), but can be redefined in the
236* sequence GCKMAX. Important for LHC studies.
237* Thanks to Harm Fesefeldt.
238*
239* Correction in GPRINT to call correctly GPDIGI and GPSETS.
240* Thanks to Michel Maire.
241* Introduced in the 314 correction cradle.
242*
243*===> 22/JAN/1991
244*
245* New facilities of DZDOC introduced in the interactive menu.
246* Patchy flag (DZDOC) to assure backward compatibility.
247* Thanks to O.Schaile
248*
249*===> 31/JAN/1991
250*
251* New facility using the JUMPT package to call user routines.
252* Patchy flag (USRJMP) to assure bacward compatibility.
253* Thanks to F.Carminati
254*
255*===> 05/FEB/1991
256*
257* Correction in GINIT to print the level of the correction
258* cradle applied. Thanks to F.Carminati.
259* Introduced in the 314 correction cradle.
260*
261* Correction Cradle 3.14/06 stamped and released.
262*
263*===> 06/FEB/1991
264*
265* Introduction of the AIX370 flag in GEANX and GEANT.
266* Thanks to Roger Howard.
267*
268*===> 13/FEB/1991
269*
270* Protection in GTRACK, GFTRAC and GLTRAC, do not update
271* pointers to cross section banks in the vacuum.
272* Thanks to C.Fuglesang, S.Banerjee, M.Maire.
273* Introduced in the 314 correction cradle.
274*
275*???> Backward incompatibility
276*
277* GFDETH returns NAMESH as a character*4 array and
278* GFDETD returns NAMESD as a character*4 array.
279* Thanks to M.Maire.
280*
281*
282*===> 14/FEB/1991
283*
284* Mods in GPMATE, GPPART, GPTMED, GSDK, GSMATE, GSMIXT, GSPART,
285* GSTMED, GSROTM, GPROTM
286* to warn against replacement of existing objects. Thanks to
287* Steve O'Neale.
288* Introduced in the 314 correction cradle.
289*
290* Force load of user routines in GINIT.
291* Introduced in the 314 correction cradle.
292*
293*===> 15/FEB/1991
294*
295* GSDVN to notify correctly when mother volume does not exists.
296* Thanks to Victor Perevoztchikov.
297* Introduced in the 314 correction cradle.
298*
299*===> 20/FEB/1991
300*
301* Correction in GHEINI for the exponent range to be compatible
302* with IEEE machines. Original limits were for IBM/370.
303* Routine LIMDAT eliminated. Thanks to Hans-Jochen Trost.
304* Introduced in the 314 correction cradle.
305*
306*???> Backward incompatibility
307*
308* Parameter NAMATE returned from GFMATE now a CHARACTER*20
309* variable. Thanks to Michel Maire.
310*
311* Monitoring introduced, flag MONITOR. Thanks to F.Carminati.
312* Mods in GXINT, GINIT, GLAST. Routine GEAMON introduced.
313* Introduced in the 314 correction cradle.
314*
315*===> 22/FEB/1991
316*
317* GWORK in GINIT was (DUMMY) called without arguments.
318* GUDTIM was (dummy) called as a routine and it is a
319* function. Thanks to F.Carminati.
320* Introduced in the 314 correction cradle.
321*
322* GUVIEW to specify 3rd argument as a character. Thanks to
323* Federico Carminati.
324* Introduced in the 314 correction cradle.
325*
326* Call NUCRIN with the right number of arguments (the last
327* was not used anyway). Thanks to Federico Carminati.
328* Introduced in the 314 correction cradle.
329*
330*===> 04/MAR/1991
331*
332* Correction in GMUSIG to avoid division by zero. Thanks to
333* M.Verzocchi.
334* Introduced in the 314 correction cradle.
335*
336*
337*===> 20/MAR/1991
338*
339* Corrections in GPRELM, muon e+e- soft radiation. Correction in
340* GDRELM for the ionization energy loss for muons. Thanks to
341* A.Nathaniel. Introduced in the 314 correction cradle.
342*
343*===> 21/MAR/1991
344*
345* Correction in GPHYSI to initialise the JMULOF bank even in
346* the case of vacuum. Electron tracking in vacuum was not taking
347* into account the TMAXFD parameter. Correction in GRANGI to set
348* the range to BIG in case of 0 energy loss tabulated in JLOSS.
349* Correction in GPHYSI to switch off all the mechanisms in the
350* vacuum but decay. The corresponding code in GTRACK has been
351* removed. The previous corrections for IUPD in GLTRAC, GFTRAC
352* and GTRACK have been removed as well. Thanks to D.Ward, R.Brun,
353* F.Carminati
354* Introduced in the 314 correction cradle.
355*
356* Corrections in GRANGI and in GTHADR to cure small steps in
357* tracking. Thanks to M.Maire.
358* Introduced in the 314 correction cradle.
359*
360*===> 25/MAR/1991
361*
362*???> Backward incompatibility
363*
364* GFPART returns NAPART as a CHARACTER*20 array. Thanks to R.Rui
365* Introduced in the 314 correction cradle.
366*
367*===> 26/MAR/1991
368*
369* GFTMAT will return the correct values for the vacuum, except
370* for hadrons. Thanks to P.Gumplinger.
371* Introduced in the 314 correction cradle.
372*
373* Correction in GPCXYZ when the number of mechanisms is greater
374* than 6. Thanks to Y.Iga.
375* Introduced in the 314 correction cradle.
376*
377*
378*===> 27/MAR/1991
379*
380* GDXYZ corrected to use HIGZ generic line types. Thanks to
381* M.Maire.
382* Introduced in the 314 correction cradle.
383*
384* Correction in the CDF for RZ/IN and RZ/OUT. Thanks to
385* F.Carminati.
386* Introduced in the 314 correction cradle.
387*
388*===> 04/APR/1991
389*
390* For reasons of speed GKS-type lines are used now in
391* GDCXYZ, GDXYZ. Thanks to F.Carminati.
392* Introduced in the 314 correction cradle.
393*
394*===> 18/APR/1991
395*
396* Correction of the logic for TMAXFD. Now is:
397*
398* AUTO 1: 0<TMAXFD<20 --> accepted
399* TMAXFD>20 --> set to 20
400* TMAXFD<=0 --> set to 20
401*
402* AUTO 0: TMAXFD>0 --> accepted
403* TMAXFD<=0 --> set to 20
404*
405* Thanks to Gerry Lynch.
406* Introduced in the 314 correction cradle.
407*
408* Correction in GMULOF. STMIN=BIG for the vacuum. Thanks
409* to F.Carminati.
410* Introduced in the 314 correction cradle.
411*
412* Correction in GRFILE, GRIN called with IDVERS=0 and not
413* 999. Thanks to F.Carminati.
414* Introduced in the 314 correction cradle.
415*
416*===> 22/APR/1991
417*
418* Corrections in GPRELA, GMULOF, GPHYSI and GTNEXT. Now
419* STMIN is set to 0. for vacuum independly of the AUTO
420* flag. Thanks to M.Maire.
421* Introduced in the 314 correction cradle.
422*
423* Format correction in GSDVN. Thanks to Federico Carminati.
424* Introduced in the 314 correction cradle.
425*
426*===> 23/APR/1991
427*
428* Correction in GINIT to initialize LIN via IQTTIN if
429* different from 0. Thanks to F.Carminati
430* Introduced in the 314 correction cradle.
431*
432*===> 24/APR/1991
433*
434* New routines GPLMAT, GPRMAT, GPGKIN from Michel Maire.
435*
436*===> 25/APR/1991
437*
438* geant.metafile is not opened any more by default.
439* The METAFILE command now accepts the metafile name
440* for packages like ATC-GKS which do not use FORTRAN.
441* Thanks to W.Koellner.
442*
443*===> 01/MAY/1991
444*
445* Problem corrected in GHEISH for neutron cascade. A neutron
446* undergoing an hadronic reaction could generate two recoil
447* protons. Thanks to Henk den Bok.
448* Introduced in the 314 correction cradle.
449*
450* Correction in GNEXT/GTNEXT. STMIN not used anymore for
451* MANY volumes. Thanks to M.Maire.
452* Introduced in the 314 correction cradle.
453*
454*===> 02/MAY/1991
455*
456* Call to NVETIM routine removed from GHEISHA. Thanks to
457* F.Carminati.
458*
459*===> 07/MAY/1991
460*
461* RNDM interactive command added to set the random number
462* generator seeds. Thanks to F.Carminati.
463*
464*===> 08/MAY/1991
465*
466* New routine GDTRAK introduced to draw a track and delete
467* it from the JXYZ buffer optionally. Thanks to R.Brun,
468* F.Carminati.
469*
470* Important correction in GTGAMA, GTNEUT to avoid problems
471* due to machine precision. When a push is made to cross a
472* volume boundary the condition that the change in coordinate
473* be not negligible is imposed. Thanks to F.Carminati, R.Brun.
474* Introduced in the 314 correction cradle.
475*
476*===> 09/MAY/1991
477*
478* Correction in GMGAUS due to G.Lynch. The sigma of the multiple
479* scattering in the gaussian approssimation was too small.
480* Introduced in the 314 correction cradle.
481*
482*===> 14/MAY/1991
483*
484* Corrections in GTELEC, GTHADR and GTMUON to introduce the
485* same boundary correction than in GTNEUT, GTGAMA. Thanks to
486* P.Gumplinger.
487* Introduced in the 314 correction cradle.
488*
489* New menu in GXINT to handle FZ files. Thanks to F.Carminati
490*
491*===> 17/MAY/1991
492*
493* Correction in GHFHDN to reset to 0. TOFD. Thanks to
494* R.Rui.
495* Introduced in the 314 correction cradle.
496*
497*===> 31/MAY/1991
498*
499* Correction in GTRACK to allow PREC to decrease and avoid
500* problems when changing particle and position but not medium.
501* Thanks to R.Spiwoks.
502* Introduced in the 314 correction cradle.
503*
504*===> 4/JUN/1991
505*
506* New version of the hidden line removal, thanks to S.Giani.
507* New menu CVOL introduced to cut volumes and new SHIFT
508* command to shift volumes for drawing.
509*
510* Correction in GHEISH. The ENP parameters are reset to 0.
511* at the beginning of every event. Thanks to Henk den Bok.
512* Introduced in the 314 correction cradle.
513*
514*===> 5/JUN/1991
515*
516* Correction introduced in GPHYSI. JTM recalculated after
517* bank push. Thanks to K.Sparks.
518* Introduced in the 314 correction cradle.
519*
520*===> 11/JUN/1991
521*
522* Correction introduced in GTELEC, GTMUON, GTHADR. When the
523* energy loss is below machine precision, it is recalculated
524* via the De/Dx table. Thanks to H.Fesefeld and Ralf Spiwoks.
525* Introduced in the 314 correction cradle.
526*
527* Booking of the bank IMAT-4 delayed till GPHYSI. Will save some
528* words in case of materials which are not inserted in a tracking
529* medium. Thanks to F.Carminati.
530*
531*???> Backward incompatibility
532*
533* The argument DMAXMS in the calling sequence of GSTMED, unused
534* in 314, will be called STEMAX and will assume the meaning of
535* maximum step allowed for a particle in the given material.
536* This will only affect people running with AUTO 0. In case of
537* automatic computation of the tracking media parameters, this
538* number will be set to BIG (=1.E10).
539* Thanks to F.Carminati, M.Maire.
540*
541*===> 12/JUN/1991
542*
543* Corrections in GMUNU to allow the usage of a variable number of
544* energy bins. Corrections in PHASP and GMUSIG for very high
545* energies. Thanks to Hans-Jochen Trost.
546* Introduced in the 314 correction cradle.
547*
548* Correction in GFLTHE to avoid division by 0. Thanks to
549* Roy Bossingham.
550* Introduced in the 314 correction cradle.
551*
552*===> 13/JUN/1991
553*
554* Correction in GINPGO and GINPCO to check correctly the PHI
555* limits. Thanks to R.Bossingham.
556* Introduced in the 314 correction cradle.
557*
558*===> 19/JUN/1991
559*
560* Correction in GMUNU to calculate correctly the angle of the
561* incoming muon. Thanks to H-J.Trost.
562* Introduced in the 314 correction cradle.
563*
564*===> 20/JUN/1991
565*
566* Corrections in GLANDZ for high energy to avoid numerical
567* problems. Thanks to H-J.Trost.
568* Introduced in the 314 correction cradle.
569*
570*===> 26/JUN/1991
571*
572* All the routines and the commons of the HADRIN/NUCRIN
573* package have been renamed. All routines begin now with
574* GHF and all commons with GCF. This to avoid interactions
575* with the interface with FLUKA. Thanks to F.Carminati
576*
577*===> 28/JUN/1991
578*
579* New algorithm for pushing a view bank. Now push of the
580* max between MORGS and the 25% of the size of the view
581* bank. Thanks to S.Giani.
582*
583*===> 17/JUL/1991
584*
585* Common HIATT of HIGZ removed from GXCONT and workstation type
586* retrieved in GXDZ from Workstation ID. Thanks to O.Couet.
587*
588* Modification in GMUSIG to protect against possible division
589* by zero. Thanks to M.Sarris.
590*
591* Correction in GDSHOW to recalculate COSPSI and SINPSI. Thanks
592* to M.Verzocchi.
593* Introduced in the 314 correction cradle.
594*
595*===> 18/JUL/1991
596*
597* Correction in GINVOL. When tracking in magnetic
598* field they could return the wrong volume. Thanks to D.Greiner.
599* Introduced in the 314 correction cradle.
600*
601*===> 24/JUL/1991
602*
603* Corrections in GRIN to allow tracking after reading a data
604* structure. Thanks to R.Brun.
605* Introduced in the 314 correction cradle.
606*
607* Modifications in GXINT to allow running with the MOTIF
608* user interface. Thanks to R.Brun.
609*
610*===> 31/JUL/1991
611*
612* INT=0 suppressed in GHEISH after nuclear fission.
613* Thanks to Kati Lassila.
614*
615*===> 01/AUG/1991
616*
617* Protection introduced in GMUSIG. Thanks to Jochen Trost.
618* Introduced in the 314 correction cradle.
619*
620*===> 02/AUG/1991
621*
622* Correction in GHSTOP to fix the calculation of the time of
623* flight for stopping particles. Thanks to L.Roberts.
624* Introduced in the 314 correction cradle.
625*
626* Correction in GNPGON not to calculate safety for the inner
627* radius when this is 0. Thanks to Andrei Nomerotsky.
628* Introduced in the 314 correction cradle.
629*
630* Correction in GNOTRP to return the correct SNXT. Thanks to
631* V.Innocente.
632* Introduced in the 314 correction cradle.
633*
634* Correction in GHSTOP not to discard hadrons at rest, they
635* may decay. Thanks to H.Fesefeld.
636* Introduced in the 314 correction cradle.
637*
638* Modifications in GHSTOP and GHEISH to handle user defined
639* particles. Thanks to P.Gumplinger.
640* Introduced in the 314 correction cradle.
641*
642*===> 12/AUG/1991
643*
644* Call to GUINTI added in GINTRI to define user commands. Thanks
645* to V.Vercesi.
646* Introduced in the 314 correction cradle.
647*
648*===> 21/AUG/1991
649*
650* Modification in GTMED to print a warning when FIELDM .EQ. 0
651* and IFIELD .NE. 0. Thanks to Federico Carminati
652*
653* Modification in GTHADR to set CFLD = BIG in case FIELDM=0.
654* Thanks to V.Vercesi.
655* Introduced in the 314 correction cradle.
656*
657*===> 27/AUG/1991
658*
659* Call to C dummy routine to initialise X11 introduced in
660* GXINT for IBM/VM. Thanks to M.Marquina.
661*
662*===> 19/SEP/1991
663*
664* Rayleigh effect now the default. IRAYL is set to 1 by default
665* in GINIT. Thanks to M.Maire
666*
667* Changes in the multiple scattering routines. GMOLI1 and GMOLI2
668* suppressed and GMOLIO used instead. Corresponding changes in
669* GMULTS and GMULOF. The code of GMOLS has been put in line inside
670* GMOLIE and GMOLS is obsolete. GPOISS and GMCOUL modified for
671* performance reasons.
672* Gaussian scattering is now generated according to the Rossi
673* gaussian formula and there is no logarithmic term in the
674* sigma of the gaussian. This gives a distribution where the
675* tails are underestimated, but which is consistent over many
676* steps.
677* Moliere and single Coulomb scattering are chosen according to
678* the value of Omega in the Moliere formula.
679* The new meaning of IMULS is the following:
680*
681* IMULS = 0 No multiple scattering
682* IMULS = 1,2 Moliere or single Coulomb scattering
683* IMULS = 3 Gaussian scattering with Rossi formula
684*
685* Thanks to G.Lynch.
686*
687*===> 14/OCT/1991
688*
689* New algorithm for calculation of energy loss. The stopping
690* range is now calculated with a parabolical interpolation
691* instead than with a linear one. This gives a dE/dx curve
692* which is a piecewise linear function and not a step function.
693* Thanks to F.Carminati.
694*
695* Modification of GNOPG1 taking care of the case of a particle
696* which is very near to the surface and it may appear inside due
697* to machine precision. This could cause the volume to be skipped.
698* Thanks to R.Nierhaus.
699*
700* Routine GNPGON rewritten for the calculation of SNXT. The
701* previous routine was returning wrong results. Thanks to
702* R.Nierhaus.
703*
704*===> 28/OCT/1991
705*
706* Modifications in the routines GFLPHI and GFCOOR. The ordering
707* for TUBS in phi was wrong. Volumes could become invisible.
708* Thanks to F.Carminati.
709*
710* New routine GFVERT introduced. Retrieves the parameters of a
711* vertex. Thanks to F.Carminati.
712*
713*===> 29/OCT/1991
714*
715* Modification in NUCREC to zero the whole of the PV array
716* to prevent the use of uninitiated variables. Thanks to
717* F.Ranjard. Introduced in the 3.14 correction cradel.
718*
719*===> 1/NOV/1991
720*
721* Useless code commented out in GTRAIN. Thanks to F.Carminati.
722*
723* GNCONE gone to double precision. Thanks to J.Toth.
724*
725*===> 5/NOV/1991
726*
727* Variable USERW undefined in the RESULT common now set to
728* UPWGHT from GCTRAK common. Thanks to F.Ranjard.
729*
730* Type declarations for GCSTAK completed. Thanks to M.Battle.
731*
732*===> 5/DEC/1991
733*
734* New calculation of the range table. Simpson integration rule
735* used. Modifications in GRANGI, GCOEFF.
736*
737*===> 10/DEC/1991
738*
739* Protection introduced against the reading of a pre-315 data
740* structure. The value of STEMAX is set to BIG in this case.
741* Modification in GPHYSI.
742*
743* Better handling of version numbers in I/O operations.
744* Modifications in GPHYSI and GRIN.
745*
746*===> 08/JAN/1992
747*
748* Correction in GLANDZ to avoid gaussian distribution for
749* very thin layers. The version of GEANT 3.13 has been
750* used for this. Thanks to F.Carminati and M.Maire.
751*
752* Corrections in GTNEXT, GINVOL and GTMEDI to protect against
753* wrong values of INGOTO when using MANY volumes. Thanks to
754* R.Brun.
755*
756*===> 15/JAN/1992
757*
758*???> Backward incompatibility
759*
760* The random numbers seeds are stored at the end of every event
761* in the JRUNG bank at locations 19 and 20. If the JRUNG data
762* structure is read in and the data card RNDM or the interactive
763* command RNDM has not been issued (NRNDM(1), NRNDM(2) <> 0),
764* and if the words 19/20 are not 0, then the random number
765* generator is restarted with these seeds. Thanks to F.Carminati
766*
767* The RNDM command now reads the values of the seeds into
768* NRNDM(1) and NRNDM(2) in common GCFLAG. The values 0 0
769* can be used for the random seeds. These values will not
770* alter the current status of the random number generator,
771* but zeroing the variables NRNDM(1) and NRNDM(2) will
772* allow them to be reinitialised with the values stored
773* in a data structure read from disk. Thanks to F.Carminati
774*
775* Routine GREND now needs an integer as input. This is to be
776* consistent with GRFILE.
777*
778*===> 27/JAN/1992
779*
780* The following changes made to the names of FLUKA routines
781* to avoid clashes with ISAJET. Thanks to L.Roberts.
782*
783* Type Original name New name
784* Routine DECAY FKDECA
785* Routine FLAVOR FKFLAV
786* Routine SIGINT FKSIGI
787*
788*===> 31/JAN/1992
789*
790* Corrections in GMEDIA in case of many volumes to avoid
791* program crash. Thanks to R.Jones.
792*
793* Change of logics in GRUN. Now if NEVENT is <= 0, no event
794* is processed. Thanks to B.Lockman.
795*
796* Common GCFLAX put in the GCFLAG sequence with the BATCH and
797* NOLOG variable from GXINT. Thanks to B.Cole.
798*
799* Problem corrected in GDECAY. The mass of the particle was
800* altered. If the next particle was the same, the mass was
801* not reset to its correct value. Thanks to S.Tonse.
802*
803* TIMINT is now really the time left after initiatlization
804* as is specified in the documentation and not the time
805* USED for initialization as it was till now.
806* Thanks to V.Ivanov.
807*
808* Corrections in GNPGON. Improvement of code safety under
809* optimization in GNOGO1, GNPGON, GNPCON. Thanks to Y.Iga.
810*
811*===> 4/FEB/1992
812*
813* Corrections in GMULOF and GTHADR. SFIELD takes precedence
814* on STMIN in case IFIELD=1 to make sure that the Runge-Kutta
815* approssimation still works. Thanks to R.Hawkings.
816*
817*===> 10/FEB/1992
818*
819* Correction in GDRAY to improve the precision in the calculation
820* of the angle. Thanks to F.Carminati & P.Lubrano.
821*
822* Change in GFKINE. The variable TOFG is not updated any more.
823* Thanks to F.Carminati
824*
825* New GMEDIA, GTMEDI, GINVOL to take care of the problems with
826* MANY volumes. Thanks to R.Jones.
827*
828*===> 14/FEB/1992
829*
830* Bug corrected in GNPGON when the particle was exactly on the
831* wall of the last Z section. Thanks to V.Palichik.
832*
833*===> 17/FEB/1992
834*
835* Updated routines GLUND, GLUNDI and new sequences LUDAT1, LUDAT3,
836* LUJETS. The data cards MSTE, KTYP, PMAS, PWID, IDB have been
837* removed because either obsolete or not aplicable. The new code
838* runs with JETSET 7.3 upward. Thanks to F.Carminati, T.Sjostrand.
839*
840* Update of GLUDKY to work with 7.3. Thanks to F.Carminati.
841*
842*===> 23/FEB/1992
843*
844* New subroutine GFIN to handle sequential input. Routine GGET
845* has been maintained for backward compatibility.
846* Thanks to F.Carminati.
847*
848*===> 01/MAR/1992
849*
850* New Runge-Kutta integration routine for the tracking in
851* magnetic field. If the result of the stepping is not accurate,
852* the step is divided in 2 parts and the integration repeated
853* and so on. Thanks to V.Perevotchikov.
854*
855*===> 05/MAR/1992
856*
857* Correction in GLTRAC. When the particle fetched has an entry
858* in JKINE the correct vertex number is calculated. Thanks to
859* Y.Foka.
860*
861*===> 06/MAR/1992
862*
863* Correction in GTELEC. The calculated range may be slightly
864* larger than the maximum allowed range due to precision
865* problems and this was leading to very small negative steps.
866* Thanks to R.Brun.
867*
868*===> 08/MAR/1992
869*
870* New GRKUTA, GHELIX, GHELX3 from V.Perevotchikov.
871*
872*===> 16/MAR/1992
873*
874* ENERGY renamed to FKENER in FLUKA. Thanks to F.Carminati
875*
876*===> 17/MAR/1992
877*
878* VERTEX->FKVERT, ZEROIN->FKZERO, ERROR->FKERRO in FLUKA.
879* Thanks to F.Carminati
880*
881* Corrections in PBANH. Thanks to M.Sasaki.
882*
883*===> 19/MAR/1992
884*
885* Corrections in GFIN, GFOUT and GRIN, version 0 is now allowed
886* for a data structure in I/O. Thanks to B.Cole.
887*
888* Modification in GBREME to improve precision for small angles.
889* Thanks to F.Carminati.
890*
891* Inlining of rotation routines in GFTRAC, GINVOL, GTMEDI,
892* GTNEXT, GMEPOS, GMEDIA. Thanks to D.Kryn, F.Carminati.
893*
894* New routines GDLENS, GDPLST and GDPRTR in the drawing package.
895* New version of GDTREE. Thanks to S.Giani.
896*
897* Streamlining of GDTOM, GINROT, GITRAN, GMTOD, GRMTD, GRMUL,
898* GROT, GTRMUL and GTRNSF. Thanks to F.Carminati.
899*
900*===> 23/MAR/1992
901*
902*???> Backward incompatibility
903* Tracking medium name in GFTMED changed in CHARACTER variable.
904* Thanks to V.Perevotchikov.
905*
906* IMPULS renamed to FKIMPU
907* DRES renamed to FKDRES
908* ERUP renamed to FKERUP in FLUKA. Thanks to F.Carminati.
909*
910*===> 24/MAR/1992
911*
912* COSI entry point renamed to FKCOSI
913* POLI renamed to FKPOLI in FLUKA. Thanks to F.Carminati
914*
915* Changes in GCOMP to increase the precision of the rotation
916* of the scattered photons and electrons in the reference
917* frame of the incoming particle. Thanks to F.Carminati
918*
919* New version of the routines GNOPG1 and GNPGO1. Speed up of
920* a factor two or more achieved. Thanks to F.Carminati,
921* M.Roethlisberger.
922*
923*===> 29/MAR/1992
924*
925* Improvement of the routines GMEDIA, GTMEDI and GINVOL. A better
926* use is made of the variable INGOTO. Thanks to R.Jones, F.Carminati.
927*
928*===> 01/APR/1992
929*
930* Modification in GTRACK. If a particle tries for more than 5 times
931* consecutively to exit a volume, the precision used for tracking is
932* multiplied by 5 and so on every fifth attempt. Thanks to
933* F.Carminati.
934*
935* Modification in FLUFIN to normalise in double precision the
936* direction cosines given to FLUKA. Thanks to A.Ferrari
937*
938* Modification in GLANDZ to avoid peaks in energy loss for very
939* light materials. Thanks to F.Carminati
940*
941* New routine GETVER to crack the title sequence and to return
942* the correct version number. Thanks to F.Carminati and M.Maire.
943*