New version SGV 2.31
Functional Changes:
-
A new generator mode has been introduced. It is now possible generate a sample single particles, with equidistant steps momentum, in azimuthal angle, or in polar angle. For detector studies this is useful, as one can simply print out the momentum (or the angle as the case might be) together with a number of properties of the track. This way one gets a list the propoerties ordered in the varying variable. This list can then directly be plotted in PAW.
The patchy-steering flag to activate this mode is "+USE,SCAN", and the variable to scan, and the range to scan is given by steering cards, see 4.1.7 in the User's Guide. A sample steering-file is given in Appendix 7
Enhancements :
-
The feauture to restart the simulation from the previously saved state of the random-generators is now implemented also on GNU/Linux (previously, this feautre was only available on VMS). If fort.40 is present, and is a link to an existing file, this file and the one that fort.41 is linked to are assumed to contain the state of the Lund and RANMAR random-number generators, respectively. The contents will be read, and will replace the state in memory. Note that the SGV by default writes out the state of the generators every 100 generated events (to unit 30 and 31). Hence, a crash or infinite-loop can be quickly reproduced even if it occured after many hours of running: Simply rename the two output-files, soft-link them to fort.40 and fort.41 as appropriate, and restart SGV.
This modification implies that the code in ZXLOOP has changed.
NB. Presently, there seems to be a problem with PYTHIA6.1(36) in this respect: If ISR is activated (MSTP(11)=1), then the same event is NOT generated. In fact, in the first event generated after the status is read in, ISR will not be present. Therefore the sequence imiteatly gets de-syncronised. In the second event, ISR is again generated, but the damage is already done. PYTHIA5, PYTHIA6.2(05) and JETSET works correctly.
-
The event-printout produced for the first MAXPRNT events now also prints out the true event as seen in the analysis part of the code (in ZAUCLU ). This differs from the printout in the generation step in that only particles with non-zero life-time are retained (in the LUJETS sense, ie. particles with either K(*,1) == 1 or with V(*,5) # 0., corresponding to stable particles and weakly decaying ones, respectively). The list, however, is in the same order as the original generated event. It is followed by a cross-reference table between this original order, and the SGV order (which list all seen particles before non-detected ones). With is information, it is possible to with certainty correlate the generated particles with the seen ones.
-
The sample title file for the single particle now is set up to vary momentum and polar angle.
-
The sample assignment script now constructs the name of the hbook-file from the arguments: it is placed in the directory given by argument 2 and with the name given by (argument 1)_sgv.hbook. On VMS, the kumac that creates the file has been modified accordingly (on UNIX, no modification is needed).
Bug fixes :
-
The treatment of mesons was wrong in ZAUJID . This has been fixed. The jet identifcation should now work correctly (Thanks to Volker Adler).
-
A protection has been introduced in ZDORD to avoid an array over-write if to many particles are defined with the SAVE_PRTICLES steering card (Thanks to Volker Adler).
-
In ZADCHC , the variable PARCL1 was not correctly set for PARCL less than 30. This has been fixed. (Thanks to Volker Adler).
-
The syntax of the system call setting up the title-file used by the detector plot kumac file showdet.kumac did not work if the user-shell was csh. This has been fixed. (Thanks to Volker Adler).
-
In ZXINI , the STEER array was not initialised to zero. This has been fixed. (Thanks to Volker Adler and Pawel Luzniak).
-
In ZAUCLU , the jet polar angle was not calculated correctly. This has been fixed. (Thanks to Pawel Luzniak).
-
Also in ZAUCLU , the counter of the number of charged tracks in the jet was incorrect. This has been fixed. (Thanks to Pawel Luzniak and Volker Adler).
Documentation :
-
The explanation of the meaning of the track parameters in section 5.2 of the User's Guide has been corrected (Thanks to Pawel Luzniak).
-
The users guide and reference manual has been updated due to the changes above. Some typos have been fixed.