Data Bank > Computer program services
Computer Programs
NEA-0983 MARTHA.
last modified: 14-MAR-1986 | catalog | categories | new | search |

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.