CONS101 Retrieve material cross-sections and stopping power

CALL GFTMAT (IMATE,IPART,CHMECA,KDIM,TKIN,VALUE*,PCUT*,IXST*)

Retrieve and interpolate the $dE/dx$ and cross-sections tabulated in material banks corresponding to a given array of kinetic energy bins.

IMATE

( INTEGER) material number;

IPART

( INTEGER) particle number;

CHMECA

( CHARACTER*4) name of the mechanism bank to be retrieved;

KDIM

( INTEGER) dimension of the arrays TKIN and VALUE;

TKIN

( REAL) array of kinetic energies in GeV;

VALUE

( REAL) array of energy losses in MeV $cm-1$ , or macroscopic cross sections in $cm-1$ ;

PCUT

( REAL) array containing the 5 energy thresholds for the material in GeV ( CUTGAM, CUTELE, CUTNEU, CUTHAD, CUTMUO);

IXST

( INTEGER) return code:

0

the table was not filled:

1

normal return; values are returned into the VALUE array;

The mechanisms have the following conventional name and a code number ( IMECA):

HADF ( 1)

total hadronic interaction cross-section according to FLUKA;

INEF ( 2)

hadronic inelastic cross-section according to FLUKA;

ELAF ( 3)

hadronic elastic cross-section according to FLUKA;

NULL ( 4)

unused;

NULL ( 5)

unused;

HADG ( 6)

total hadronic interaction cross-section according to GHEISHA;

INEG ( 7)

hadronic inelastic cross-section according to GHEISHA;

ELAG ( 8)

hadronic elastic cross-section according to GHEISHA;

FISG ( 9)

nuclear fission cross-section according to GHEISHA;

CAPG (10)

neutron capture cross-section according to GHEISHA;

LOSS (11)

stopping power;

PHOT (12)

photoelectric cross-section;

ANNI (13)

positron annihilation cross-section;

COMP (14)

Compton effect cross-section;

BREM (15)

bremsstrahlung cross-section;

PAIR (16)

photon and muon direct- pair cross-section;

DRAY (17)

$\delta$ -rays cross-section;

PFIS (18)

photo-fission cross-section;

RAYL (19)

Rayleigh scattering cross-section;

MUNU (20)

muon-nuclear interaction cross-section;

RANG (21)

range in cm;

STEP (22)

maximum step in cm.

Note: Common /GCMULO/ contains an array ELOW(200) (see [CONS199]) with NEK1 kinetic energy values ranging from EKMIN to EKMAX. GPHYSI initialises by default the first 91 locations of ELOW with values of energy from $10-5$ GeV (10 keV) to $104$ GeV (10 TeV) equally spaced on logarithmic scale. This can be controlled by the user via the data record or the interactive command ERAN ( [BASE040]). ELOW can be used as the input argument TKIN, e.g.:

CALL GFTMAT (10, 1,'PHOT',NEK1, ELOW ,VALUE, PCUT, IXST)

CALL GPLMAT (IMATE,IPART,CHMECA,KDIM,TKIN,IDM)

Plot and interpolate the $dE/dx$ and cross-sections tabulated in material banks corresponding to a given array of kinetic energy bins. The mechanisms are the same as the previous routine. If the mechanism name is 'ALL', then all the mechanisms will be plotted.

The first five parameters have the same meaning as in the previous routine and:

IDM

( INTEGER) treatment of the created histogram(s):

>0

fill, print, keep histogram(s);

=0

fill, print, delete histogram(s);

<0

fill, noprint, keep histogram(s).

CALL GPRMAT (IMATE,IPART,CHMECA,KDIM,TKIN)

Print and interpolate the $dE/dx$ and cross-sections tabulated in material banks corresponding to a given array of kinetic energy bins. The mechanisms are the same as the previous routine. If the mechanism name is 'ALL', then all the mechanisms will be printed.

The result is printed in the form of a table.

The five parameters have the same meaning as in the previous routine.

M.Maire CONS110