Last version:
, October 16th 2024 (last respin 2024.1.2) 06-May-2024
|
[ <--- prev -- ] [ HOME ] [ -- next ---> ] [ full index ]
PART-THRes
sets different energy transport cutoffs for hadrons, muons and
neutrinos
See also EMFCUT, LOW-BIAS, THRESHOLd
The meaning of WHAT(1) depends also on the value of WHAT(5).
For WHAT(5) = 0.0 :
WHAT(1) < 0.0 : kinetic energy cutoff (GeV)
> 0.0 : momentum cutoff (GeV/c)
For WHAT(5) >= 1.0 :
WHAT(1) < 0.0 : gamma cutoff (Lorentz factor, = E/mc**2)
> 0.0 : eta cutoff (= beta*gamma = v/c E/mc**2)
Default (WHAT(1) = 0.0): the cutoff is 0 for neutrinos,
and 1.E-14 GeV for neutrons.
For any other hadrons, and for muons:
if option DEFAULTS is missing, or is present with
SDUM = NEW-DEFAults or SHIELDINg, the default cutoff
kinetic energy is 0.01 GeV.
If SDUM = HADROTHErapy, ICARUS or PRECISIOn, the default
cutoff kinetic energy is 0.0001 GeV.
If SDUM = CALORIMEtry, the default cutoff kinetic energy
is = 0.001 * m/m_p GeV (m = particle mass, m_p = proton mass)
In any other case, the default cutoff is 0.050 GeV.
(For e+e- and photons the threshold is set by EMFCUT, see
Note 3 below).
WHAT(2) = lower bound of the particle id-numbers to which the cutoff
applies ("From particle WHAT(2)...").
Default = 1.0
WHAT(3) = upper bound of the particle id-numbers to which the cutoff
applies ("...to particle WHAT(3)...").
Default = WHAT(2)
WHAT(4) = step length in assigning numbers
("...in steps of WHAT(4)")
Default = 1.0.
WHAT(5) : depending on its value, cutoff values indicated by
WHAT(1) are assigned to kinetic energy, momentum, gamma or
eta (see WHAT(1))
WHAT(6) = 1.0 restricts the given cutoff to charged particles only
Default: the cutoff applies to all particles indicated by
WHAT(2-4)
SDUM : not used
Default (option PART-THR not given): thresholds as described above
for WHAT(1) = 0.0.
Notes:
- 1) If low-energy neutron transport is not requested (explicitly via
LOW-NEUT or implicitly via DEFAULTS), the energy of neutrons
below threshold is deposited on the spot.
- 2) The total momentum cutoffs of heavy ions are derived from that
of a 4He ion (4-HELIUM) by scaling the latter with the ratios
of the atomic weights of the heavy ions and the 4He ion.
The total momentum cutoffs for light ions (4-HELIUM, 3-HELIUM,
TRITON and DEUTERON) can be defined by PART-THRes. If this is not
done, they are derived from that one of a proton by scaling the
latter with the ratios of the atomic weights of the light ions and
a proton.
- 3) Option PART-THR acts on all particles excepted e+ e- and
photons, while EMFCUT option is used to set to transport electrons,
positrons and photons.
- 4) When the energy of a heavy charged particle becomes lower than the
cutoff defined by PART-THR, and if such cutoff is lower
than 100 MeV, the particle is not stopped, but is ranged
out to rest in an approximate way. Its kinetic energy is
deposited uniformly over the residual range if the latter
is contained within a single region; otherwise a new residual
range is calculated at each boundary crossing and the residual
kinetic energy is distributed accordingly. If applicable, such
a particle eventually decays at rest or is captured. All other
forms of transport are ignored excepted curved paths in
magnetic fields (multiple scattering, delta ray production,
inelastic or elastic collisions, and INCLUDING DECAY IN
FLIGHT). Magnetic fields are taken into account but only
very roughly, since the continuous slowing down of the
particles is not simulated. Antiprotons and pi-minus are
always ranged out to rest (without allowance for decay) and
made to annihilate on a nucleus.
- 5) If the cutoff is higher than 100 MeV, however, the
particles are stopped in place without any further treatment.
If this happens at a boundary crossing where the material of
the region entered is vacuum, a printed message warns the
user that energy is being deposited in vacuum.
- 6) By default the neutron threshold is set at 1.E-14 GeV (1.E-5 eV, the
lowest boundary of the group structure). So, normally it is not
necessary to issue a PART-THR command at all for neutrons.
A note of caution: if a PART-THR has been issued spanning all
particles, it is generally necessary to override it with another one
resetting the threshold for neutrons to 1.E-14 GeV.
As a general rule however, if a neutron transport threshold is
set < 20 MeV, it is rounded to the closest lower group boundary.
Example:
* A threshold of 2 MeV (kinetic energy) is requested for heavy charged
* particles with id-numbers between 1 and 11 (protons, antiprotons and
* muons). A threshold of Gamma (E/m) = 2 will apply for pions and kaons
* (numbers from 13 to 16). For all other particles, the defaults will apply.
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
PART-THR -0.002 1.0 11.0 0.0 0.0 1.0
PART-THR -2.0 13.0 16.0 1.0 1.0 0.0
|