allows to override the standard FLUKA defaults
for physics processes.
(see also DPMJET, EMFCUT, EVENTYPE, POLARIZA, RQMD, THRESHOL)
This codeword concerns the following physics processes:
1) SDUM = DECAYS: decays of pions+/-, muons+/-, Kaon+/- (e-nue, mu-numu,
K+/-mu3 and K+/-e3 channels) and Klong (K0mu3 and K0e3 channels)
2) SDUM = CHARMDECay: flag for charmed hadron transport
3) SDUM = COALESCEnce: flag to activate the coalescence mechanism
4) SDUM = DPMTHREShold: lower energy threshold(s) for DPMJET
5) SDUM = EM-DISSOciation: ion electromagnetic-dissociation
6) SDUM = EVAPORAT: evaporation
7) SDUM = IONSPLITting: activates the superposition model, \ie\ ion
splitting into nucleons
For SDUM = CHARMDECay:
WHAT(1) : flag for charmed hadron decays
=< -1.0 : resets to default (decay at production, no transport)
= 0.0 : ignored
= 1.0 : charmed hadrons are transported
WHAT(2)-WHAT(6) : not used
For SDUM = COALESCEnce:
WHAT(1) : coalescence flag
=< 0.0 : false (no coalescence, default)
= 0.0 : ignored
> 0.0 : true, coalescence is activated
WHAT(2)-WHAT(6) : reserved to developers' use
For SDUM = DECAYS:
WHAT(1) : flag for particle decay
=< -1.0 resets to default
= 0.0 ignored
= 1.0 maximum accuracy, polarisation accounted for in pi/K
--> mu(e)-numu(nue) decays and following mu decays (default)
= 2.0 maximum accuracy, polarisation not accounted for
= 3.0 phase space like decays
200 > WHAT(1) >= 100 : leptonic decays only are allowed (implemented for
tau's only). WHAT(1) - 100 has the same meaning as above
300 > WHAT(1) >= 200 : hadronic decays only are allowed (implemented for
tau's only) WHAT(1) - 200 has the same meaning as above
WHAT(2) and WHAT(3): not used
WHAT(4) = lower bound of the particle numbers to which the decay flag
chosen by WHAT(1) applies
("From particle WHAT(4)...").
Default = 1.0
WHAT(5) = upper bound of the particle numbers to which the decay flag
chosen by WHAT(1) applies
("...to particle WHAT(5)...").
Default = WHAT(4)WHAT(6) = step length in assigning numbers
("...in steps of WHAT(6)")
Default = 1.0.
For SDUM = DPMTHREShold:
WHAT(1) : minimum DPMJET kinetic energy for hadrons (GeV)
=< 0.0 : ignored
WHAT(2) : minimum DPMJET kinetic energy for ions (GeV/n)
=< 0.0 : ignored
WHAT(3) : minimum RQMD kinetic energy for ions (GeV/n)
=< 0.0 : ignored
WHAT(4) : smearing (+/- Delta E, GeV/n) for the RQMD-DPMJET switch energy
< 0.0 : resets to default (0.0)
WHAT(5) : smearing (+/- Delta E, GeV/n) for the FLUKA-DPMJET switch energy
for hA interactions
< 0.0 : resets to default (0.0)
WHAT(6) : flag for restricting DPMJET h-A interactions to primary particles
only
=< -1.0 : resets to default (false)
=< 0.0 : ignored
> 0.0 : sets to true
Default (no DPMTHREShold option): DPMJET is not called for h-A interactions,
but it is called for A-A interactions down to 5 GeV/n. RQMD is
called between 5 and 0.1 GeV/n.
WARNING: to activate interactions of primary and secondary heavy ions,
the EVENTYPE card MUST be issued with SDUM = DPMJET
WARNING: the FLUKA executable must be built with the DPMJET and RQMD
libraries to perform A-A interactions (see ldpmqmd script in $FLUPRO/flutil)
For SDUM = EM-DISSOciation:
WHAT(1) : flag for activating ion electromagnetic-dissociation
=< -1.0 : resets to default (no em-dissociation)
= 0.0 : ignored
= 1.0 : (default) no em-dissociation
= 2.0 : projectile and target em-dissociation activated
= 3.0 : projectile only em-dissociation activated
= 4.0 : target only em-dissociation activated
WHAT(2)-WHAT(6): not used
For SDUM = EVAPORATion:
WHAT(1) : flag for FLUKA evaporation model
=< -1.0 : resets to default (new model, no heavy fragment
evaporation)
= 0.0 : ignored
= 1.0 : old evaporation model (OBSOLETE: kept for developers' use
only)
= 2.0 : new evaporation model, no heavy fragment evaporation
= 3.0 : new evaporation model, with heavy fragment evaporation
(CPU expensive, see Note below)
WHAT(2)-WHAT(6): not used
For SDUM = IONSPLITting:
WHAT(1) : flag for activating ion splitting into nucleons
< 0.0 : false, no ion splitting (default)
= 0.0 : ignored
> 0.0 : true: ion splitting is activated
WHAT(2) : minimum energy for ions (GeV/n) above which splitting into
nucleons will be performed (default: 0.1 GeV/n)
=< 0.0 : ignored
WHAT(3) : maximum energy for ions (GeV/n) below which splitting into
nucleons will be performed (default: 5 GeV/n)
=< 0.0 : ignored
WHAT(4)-WHAT(6) : not used
Default (option PHYSICS not given): standard FLUKA treatment of
physics processes
Note: In order to achieve accurate results for residual nuclei production
the evaporation of heavy fragments MUST be activated. This, however,
is not the default since it brings a heavy CPU burden, and is not
needed for most applications.
Example:
* Only hadronic decays are allowed for tau+ and tau- (id-number 41 and 42)
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+...
PHYSICS 201.0 0.0 0.0 41.0 42.0 0. DECAYS
* Maximum accuracy requested for decay of pi+ and pi-(id-number 13 and 14),
* but without accounting for polarisation
* Phase space
PHYSICS 2.0 0.0 0.0 13.0 14.0 0. DECAYS
* New evaporation model requested
PHYSICS 2.0 0.0 0.0 0.0 0.0 0. EVAPORAT
