WARNING! Obsolete version!
Multiple updates and restructuring has been incorporated over VIKAR 3.X. Most significantly, the
reaction geometry has been reworked to be fully rigorous, thereby allowing jet-style targets
to be simulated correctly. This also removes some approximations used in propogating particles
through the two layers of a telescope. Functional improvements include:
Realistic simulation of jet and foil targets
User-defined energy thresholds for detectors
Selectable position-dependent/independent energy thresholds for resistive/non-resisitve detectors in cylindrical coordinates
User-defined offset resisitors for resistive-strip detectors (which impact position resolution and thresholds)
Separate energy thresholds for dE and E detectors
Incomplete charge collection included in detector response (optional)
User-defined beam emergence
Internal Q-value calculations (the code Qvalue is still bundled separately for controlled post-processing)
Basic support for massless and neutral beam particles (for simulation of (n,X) and (g,X) reactions)
Automatic xmgrace formatting on output files (optional)
Up-to-date detector files are bundled, including files for ORRUBA, GODDESS, and superORRUBA with JENSA
Details of file changes are documented in below.
1) Subroutine add_phi (v1.0) created, for rotating the azimuthal angle of a vector
in polar coordinates by an angle phi, from specific to absolute coordinates
2) Subroutine add_phi (v1.0) created, for rotating the azimuthal angle of a vector
in polar coordinates by an angle phi from absolute to specific coordinates
3) Subroutine DetSet_cyl_read (v1.1) updated to include separate energy detection thresholds
for dE and E detectors, offset resistors (note that strips resistance is nominally 4 kOhm) and
telescope separation.
4) Subroutine DetSet_plan_read (v1.1) updated to include separate energy detection thresholds
for dE and E detectors, offset resistors (note that strips resistance is nominally 4 kOhm) and
telescope separation.
5) Subroutine DetSet_ann_read (v1.1) updated to allow X,Y translation of detector, separate energy
detection thresholds for dE and E detectors, and telescope separation.
6) Subroutine jet_out (v1.0) created, to determine the thickness that a particle
sees when exiting a jet-like target.
7) Subroutine detect_cyl (v1.0) created to determine whether a cylindrical polar detector
is hit. This is a substantial upgrade, allowing event-by-event variations in reaction
coordinate (imperative for JENSA simulations). Resolution-smeared vectors are passed back
to parent program, plus a hit register.
8) Subroutine DetHit_plan updated to v1.1 to include reaction point offset
(due to beam and target size) and to return x,y positions of hit on detector
9) Subroutine DetHit_ann (v1.0) created to determine whether an annular detector is hit.
This is a substantial upgrade, allowing event-by-event variations in reaction coordinate
(imperative for JENSA simulations). Resolution-smeared vectors are passed back to parent
program, plus a hit register.
10) Subroutine det_thick_cyl updated to v2.0 which uses subtract_phi function to more generally
calculate effective detector thickness
11) Subroutine strag_dE_ann updated to v2.0 to properly propagate an ion that has hit the dE
detector through to the E detector, accounting for the non-central interaction point in the
target, and the non-zero distance between dE and E layers in the telescope.
12) Subroutine strag_dE_plan updated to v2.0 to properly propagate an ion that has hit the dE
detector through to the E detector, accounting for the non-central interaction point in the
target, and the non-zero distance between dE and E layers in the telescope.
13) Subroutine strag_dE_cyl updated to v2.0 to properly propagate an ion that has hit the dE
detector through to the E detector, accounting for the non-central interaction point in the
target, and the non-zero distance between dE and E layers in the telescope.
14) new_resolution_cyl (v1.0) created to determine the resolution effects of hitting a
detector in cylindrical coordinates. Detector thresholds, incomplete charge
collection, and detector resolution effects are applied.
15) new_resolution_plan (v1.0) created to determine the resolution effects of hitting a
detector in planar coordinates. Detector thresholds, incomplete charge
collection, and detector resolution effects are applied.
16) new_resolution_ann (v1.0) created to determine the resolution effects of hitting a
detector in annular coordinates. Detector thresholds, incomplete charge
collection, and detector resolution effects are applied.
17) Qcal (v1.0) created, based on Qvalue code, to generate reconstructed excitation energy spectra
from the output files. Running this spectrum creation is optional. Spectrum dispersion is
selectable.
18) product_proc updated to v1.2 for several major alterations. Both jet and foil targets
are now handled. Separate dE and E positions are determined, allowing for offset in separation
between dE and E detector (for non-central reactions and dE straggling). Routine now uses
the detect_cyl routine, and includes incomplete charge collection in detectors (on/off selectable).