GCALOR: A hadronic simulation package for GEANT3

Welcome to the GCALOR Simulation Package

GCALOR general remarks

The GCALOR package contains the original CALOR simulation code (including the low energy neutron code MICAP ) developed and maintained by the Oak Ridge National Laboratory ( ORNL ) and an interface to utilize this code in the framework of GEANT3 .

GCALOR consists of

HETC: The High-Energy-Transport-Code is transporting charged hadrons up to an energy of 10 GeV through the materials of the setup
FLUKA: The FLUKA fragmentation model is utilized for interactions above the HETC limit
MICAP: The neutron code from the Monte-carlo-Ionization-Chamber-Analysis-Program is called for neutrons with a kinetic energy below 20 MeV. The simulation of interactions is based on measured/calculated cross sections and secondary particle energy/angular distributions
Interface: The interface between CALOR and GEANT determines which code to call for the current particle. The material description is automatically extracted from GEANT

The transport of electromagnetic interacting particles is performed by GEANT as well as the tracking and the calculation of the energy deposition.

The code has been The code is maintained by C. ZeitnitzUniversity of Wuppertal) and T.A. Gabriel (ORNL).

How to use the GCALOR package in the frame of GEANT 3

The two user subroutines GUPHAD and GUHADR have to be changed in order to retrieve the interaction probability (cross-section) of the current particle with the material it passed through.

      SUBROUTINE GUPHAD
C get interaction probability
C     CALL GPGHEI
C USE GCALOR INSTEAD OF GHEISHA
      CALL CALSIG
      RETURN

      SUBROUTINE GUHADR
C inelastic hadronic interaction
C     CALL GHEISH
C USE GCALOR INSTEAD OF GHEISHA
      CALL GCALOR
      RETURN

The GCALOR system is automatically initialized on first call. The initializations consists of the following items:

Establish data base for materials defined for the current setup
Read in cross-section data (expected files: chetc.dat and xsneut.dat in the CERN_ROOT/lib directory) for all needed materials
Calculate the macroscopic cross-section for low energy neutrons and write relevant informations for the user to file calor.out

GCALOR Manual, source code and data files

News about the current version of GCALOR

List of downloadable files

GCALOR manual: Postscript file
CMZ file: GCALOR CMZ cards file (gcalor.car)
Fortran Source
Example GUHADR: Subroutine to extract recoil nucleus information
TAR file: Tape ARchieve file contains everything except the data files and the manual (gcalor.tar)
DATA Files: HETC data file (chetc.dat); MICAP neutron cross-section file (xsneut.dat)

Additional Cross-Section Data Files

Here are additional files, which are not included in the data file above!

Cesium cross-section (650kBytes compressed) file including (n,G) data from BNL-database. be aware, that this file has been thinned to 10% precision in order to get the size of the file down.

Cross-Section file 1 (474kBytes compressed) file contains ENDF/B-VI data for K-39 (natural composition), Mn-55, Sr-88, Sm-150 and Tb-159

Ge Cross-Section (546kBytes compressed) file contains Germanium data for the following isotopes: Ge-70, Ge-72, Ge-73, Ge-74, Ge-76

How to include the above files:

Copy file to a local directory (for example /home/users/zeitnit/gcalor/xs_file.dat.gz) and gunzip it.
Include a call to MICFIL in your GEANT material setup routine (for the above example): CALL MICFIL('/home/users/zeitnit/gcalor/xs_file.dat')
You can have an infinite number of MICFIL calls in your program

Additional Requirements

In addition to the GCALOR code and the data files, some CERN library routines and the GEANT library are needed to complete the simulation package

Prof. Dr. Christian Zeitnitz