Computer Programs

NEA-0983 MARTHA.

NEA-0983 MARTHA.

MARTHA, Nai(Tl) Gamma Scintillation Detector Response by Monte-Carlo

NAME OR DESIGNATION OF PROGRAM, COMPUTER, DESCRIPTION OF PROGRAM OR FUNCTION, METHOD OF SOLUTION, RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM, TYPICAL RUNNING TIME, UNUSUAL FEATURES OF THE PROGRAM, RELATED AND AUXILIARY PROGRAMS, STATUS, REFERENCES, MACHINE REQUIREMENTS, LANGUAGE, OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED, OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS, NAME AND ESTABLISHMENT OF AUTHORS, MATERIAL, CATEGORIES

[ top ]

1. NAME OR DESIGNATION OF PROGRAM: MARTHA.

[ top ]

2. COMPUTERS

To submit a request, click below on the link of the version you wish to order. Rules for end-users are available here.

Program name	Package id	Status	Status date
MARTHA	NEA-0983/01	Tested	14-MAR-1986

Machines used:

Package ID	Orig. computer	Test computer
NEA-0983/01	FACOM M-380	IBM 3084Q

[ top ]

3. DESCRIPTION OF PROGRAM OR FUNCTION

The program MARTHA calculates NaI(Tl) scintillation detector response functions to gamma rays. There are four versions of the program according to the detector shape and gamma-ray energy.

[ top ]

4. METHOD OF SOLUTION

Using a Monte Carlo method, every photon and electron produced by one incident photon are traced in a detector assembly, and the total energy deposited in the crystal is recorded as a pulse height. Repeating this procedure several times, a
response function is obtained.

[ top ]

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

Response functions can be calculated for cylindrical or spherical detectors. The program was originally designed to calculate response functions for gamma energies from 0.01 MeV to 10 MeV. But the program should work up to 20 MeV. Available source conditions are a point source on the detector axis and a parallel beam. The maximum channel number of a
spectrum is 300.

[ top ]

6. TYPICAL RUNNING TIME

Running time depends among other factors on the incident gamma energy. For 10,000 gamma quants of 1.0 MeV incident on a 3 in. x 3 in. cylindrical crystal, execution time is approximately 11 seconds on an FACOM 380.

NEA-0983/01

NEA-DB executed the test cases included in this package on an IBM 3084Q computer. The following CPU times were required: 22 seconds (version 1); 15 seconds (version 2); 53 seconds (version 3); 21 seconds (version 4).

[ top ]

7. UNUSUAL FEATURES OF THE PROGRAM

The program takes the following two effects into account: (1) The housing material of the NaI(Tl) crystal; (2) the non-linear response of a NaI(Tl) detector to electron energy.

[ top ]

8. RELATED AND AUXILIARY PROGRAMS

Five programs to prepare the input data for MARTHA. These data are already given in the sample input.

[ top ]

9. STATUS

Package ID	Status date	Status
NEA-0983/01	14-MAR-1986	Tested at NEADB

[ top ]

10. REFERENCES:

NEA-0983/01, included references:

- K. Saito and S. Moriuchi:
  Monte Carlo Calculation of Accurate Response Functions for a
  NaI(Tl) Detector for Gamma Rays.
  Reprinted from: Nuclear Instruments & Methods, 185(1981)299-308
- K. Saito and S. Moriuchi:
  Monte Carlo Calculation of Accurate Response Functions for a
  NaI(Tl) Detector for Gamma Rays and Analysis of Pulse Height
  Spectrum Formation Mechanism.
  JAERI-M 9741  (In Japanese)  (September 1981)

[ top ]

11. MACHINE REQUIREMENTS:

NEA-0983/01

The test cases were run in 550K bytes of IBM 3084Q main storage.

[ top ]

12. PROGRAMMING LANGUAGE(S) USED

Package ID	Computer language
NEA-0983/01	FORTRAN-IV

[ top ]

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:
FACOM OS 4-F4-MSP.

NEA-0983/01

MVS-SP (IBM 3084Q).

[ top ]

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS:

[ top ]

15. NAME AND ESTABLISHMENT OF AUTHORS

         K. Saito
         Department of Environmental Safety Research
         Japan Atomic Energy Research Institute
         Tokai-mura, Naka-gun, Ibaraki-ken
         Japan

[ top ]

16. MATERIAL AVAILABLE

NEA-0983/01

File name	File description	Records
NEA0983_01.001	This information file	119
NEA0983_01.002	Program to convert EBCDIC to binary	26
NEA0983_01.003	Routine RANDOM from SLATEC library	117
NEA0983_01.004	Program MARTHA source. Version 1	2496
NEA0983_01.005	JCL to run MARTHA version 1	101
NEA0983_01.006	Input for version 1	57
NEA0983_01.007	Library for version 1 (file code 31)	216
NEA0983_01.008	Library for version 1 (file code 32)	594
NEA0983_01.009	Library for version 1 (file code 33)	540
NEA0983_01.010	Test case output, version 1	768
NEA0983_01.011	Program MARTHA source, version 2	1972
NEA0983_01.012	JCL to run version 2	74
NEA0983_01.013	Input for version 2	20
NEA0983_01.014	Test case output,version 2	425
NEA0983_01.015	Program MARTHA source, version 3	1968
NEA0983_01.016	JCL to run version 3	101
NEA0983_01.017	Input for version 3	48
NEA0983_01.018	Library for version 3 (file code 31)	216
NEA0983_01.019	Library for version 3 (file code 32)	594
NEA0983_01.020	Library for version 3 (file code 33)	540
NEA0983_01.021	Test case output, version 3	1076
NEA0983_01.022	Program MARTHA source, version 4	1557
NEA0983_01.023	JCL to run version 4	74
NEA0983_01.024	Input for version 4	15
NEA0983_01.025	Test case output, version 4	663

[ top ]

17. CATEGORIES

O. Experimental Data Processing

Keywords: gamma detectors, response functions, scintillation detectors.