IAEA0940 JADSPE.

1. NAME OR DESIGNATION OF PROGRAM:  JADSPE.

3. DESCRIPTION OF PROGRAM OR FUNCTION

JADSPE is a package of eight programs to process multichannel gamma-ray spectra. The programs can be used to:
  
- locate automatically spectral peaks and calculate their positions,   areas, and full widths at half maximum (FWHM);
- plot the spectra on a CALCOMP plotter, TEKTRONIX terminal or a line printer;
- add or subtract several spectra with the possibility of adjusting either their start and end channels or the maxima of the chosen corresponding peaks.
  
The JADSPE package comprises the following programs:
  
- SPECTF: automatic location of peaks and calculation of their psitions, areas and FWHMS. The standard deviations of peak parameters are also determined, and each evaluated region is plotted on the line printer.
- SPECT1: The areas amd FWHMs are calculated for peaks whose positions are known beforehand. The standard deviations of calculated parameters are also determined, and each evaluated region is plotted on the line printer.
- PLOCHA: The peak net area is calculated by summing the channel contents in specified regions and by subtracting a linear background.
- GRAPH: Spectrum plotting on the line printer.
- PLTNEW: Spectrum plotting on CALCOMP plotter or on TEKTRONIX terminal.
- SUMDIF: The channel contents of several gamma-ray spectra are added  or subtracted.
- SSPFP: The channel contents of several gamma-ray spectra are added with adjustment of the maxima of specified peaks.
- SOUCET: The channel contents of several gamm-ray spectra are added with the adjustment of start and end channels of the spectra.

4. METHOD OF SOLUTION:  Non-linear least-square fit.

5. RESTRICTIONS ON THE COMPLEXITY OF THE PROBLEM

The full energy peaks are approximated by a symmetrical Gaussian function and the underlying background is approximated by a first-order polynomial.
A fixed spectrum length of 4096 channels is assumed.
Maxima of:
- number of peaks in one multiplet: 9
- number of peaks identified by the automatic search procedure: 300
- number of channels in one evaluated interval: 500
- dimensions of CALCOMP plots: 54 cam x 36 cm
- number of spectra, which may be added by SSPFP: 5
- number of peaks for spectra adjustment in SSPFP: 20
- number of spectra added by SOUCET: 10.

6. TYPICAL RUNNING TIME

SPECTF, SPECT1: running time depends on the spectrum complexity and  the total number of peaks evaluted. The determination of peak positions, areas, and FWHM's requires from 10 to 60 ETU per spectrum (1 ETU = 3.6 seconds).
GRAPH: plotting on line printer requires aboutg 25 ETU per spectrum. PLTNEW: plotting on CALCOMP requires about 2 ETU per spectrum.
SUMDIF, SSPFP, SOUCET: running time depends on the number of spectra involved. To sum two spectra requires about 4 ETU.

It took a few milliseconds to test each program.

7. UNUSUAL FEATURES OF THE PROGRAM:

8. RELATED AND AUXILIARY PROGRAMS:

10. REFERENCES

- J. R. Wolberg:
  "Prediction Analysis"
  Princeton, D. Van Nostrand Comp, Inc. 1967
- J. Pitha, R.N. Jones:
  Internal Report
  Division of Pure Chemistry, National Research Council, Ottawa
  Canada
- J. Rikovska, V. Zvolska, V. Hnatowicz:
  Radioisotopy 13(5), 1007-1023 (1972)
- E. Stejskalova:
  Acta Polytechnica, 14 (IV,1) (1985)
- "Multijob Programming"
  TP4572, ICL
- "User's Guide to the MULTIJOB GRAFOR routine"
  University of Bath, Computer Unit

- J. Rikovska and E. Stejskalova:
JADSPE - The Program Package for Analysis of the Gamma-Ray
Energy Spectra (December 1987)

11. MACHINE REQUIREMENTS

Program execution requires from 41 to 106 kbytes of mainm storage. CALCOMP plotter or TEKTRONIX 4010 terminal are required for spectrum plotting.

IAEA0940/01 was tested at the NEA Data Bank on a PC with 3GHz INTEL Pentium 4.
Hard disk space requirement is 3.42Mbytes, including 1.5Mbytes of documentation.

13. OPERATING SYSTEM UNDER WHICH PROGRAM IS EXECUTED:
MULTIJOB (ICL/4-72).

IAEA0940/01 was tested at the NEA Data Bank on:
- WINDOWS XP Professional with g95 (g95), Compaq Visual Fortran Professional Edition 6.6 a (df), and Lahey/Fujitsu Fortran 95 Compiler Release 5.50d (lf95).
- SuSE linux 9.2 kernel 2.6.8 with g95 (g95), and Lahey/Fujitsu Fortran 95 Compiler Release 6.10a (lf95).

14. OTHER PROGRAMMING OR OPERATING INFORMATION OR RESTRICTIONS:

Some programs need an impact line printer.

15. NAME AND ESTABLISHMENT OF AUTHORS

     J. RIKOVSKA and E. STEJSKALOVA
     Technical University of Prague
     Faculty of Nuclear Science and Physical Engineering
     Department of Mathematics
     Brehova 7
     115 19 PRAGUE 1 - Stare Mesto
     Czech Republic

source program   mag tapeSPECTF Source program                      SRCTP
test-case data   mag tapeSPECTF JCL                                 DATTP
test-case data   mag tapeSPECTF Sample Problem Input Data           DATTP
source program   mag tapeSPECT1 Source program                      SRCTP
test-case data   mag tapeSPECT1 JCL                                 DATTP
test-case data   mag tapeSPECT1 Sample Problem Input Data           DATTP
source program   mag tapePLOCHA Source program                      SRCTP
source program   mag tapePLOCHA Macro (USERCODE)                    SRCTP
test-case data   mag tapePLOCHA JCL                                 DATTP
test-case data   mag tapePLOCHA Sample Problem Input Data           DATTP
test-case data   mag tapePLOCHA Sample Problem Input Data S0019     DATTP
test-case output mag tapePLOCHA Sample Problem Printed Output       OUTTP
source program   mag tapeGRAPH Source program                       SRCTP
test-case data   mag tapeGRAPH JCL                                  DATTP
test-case data   mag tapeGRAPH Sample Problem Input Data            DATTP
source program   mag tapePLTNEW Source program                      SRCTP
test-case data   mag tapePLTNEW JCL                                 DATTP
test-case data   mag tapePLTNEW Sample Problem Input Data           DATTP
source program   mag tapeSUMDIF Source program                      SRCTP
test-case data   mag tapeSUMDIF JCL                                 DATTP
test-case data   mag tapeSUMDIF Sample Problem Input Data           DATTP
test-case output mag tapeSUMDIF Sample Problem Printed Output       OUTTP
source program   mag tapeSSPFP Source program                       SRCTP
test-case data   mag tapeSSPFP JCL                                  DATTP
test-case data   mag tapeSSPFP Sample Problem Input Data            DATTP
test-case output mag tapeSSPFP Sample Problem Printed Output        OUTTP
source program   mag tapeSOUCET Source program                      SRCTP
test-case data   mag tapeSOUCET JCL                                 DATTP
test-case data   mag tapeSOUCET Sample Problem Input Data           DATTP
test-case output mag tapeSOUCET Sample Problem Printed Output       OUTTP
test-case output mag tapeSPECTF Sample Problem Printed Output       OUTTP
test-case output mag tapeSPECT1 Sample Problem Printed Output       OUTTP
test-case output mag tapeGRAPH Sample Problem Printed Output        OUTTP
report                   JADSPE Report  (December 1987)             REPPT