Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

Add Grabmayrs gamma-cascade-reader and a custom neutron capture process #123

Merged
merged 5 commits into from
Sep 12, 2024
Merged

Add Grabmayrs gamma-cascade-reader and a custom neutron capture process #123

Show file tree
Hide file tree
Changes from 2 commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
733a862
Add custom neutron capture to apply Grabmayrs gamma cascades
EricMEsch Sep 2, 2024
7a1a7c9
Add example 06 testing the custom nCapture and employing a user-defin…
EricMEsch Sep 2, 2024
e5e983d
Make reader thread-local to ensure thread-safety
EricMEsch Sep 3, 2024
83eafe1
Remove unnecessary std::move and add custom messenger in RMGGrabmayrG…
EricMEsch Sep 12, 2024
dd71c09
Apply minor style changes for consistency
EricMEsch Sep 12, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
45 changes: 45 additions & 0 deletions docs/rmg-commands.rst
View file
Open in desktop

Some generated files are not rendered by default. Learn more about how customized files appear on GitHub.

10 changes: 10 additions & 0 deletions examples/06-NeutronCapture/CMakeLists.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,10 @@
cmake_minimum_required(VERSION 3.8)
project(TestNeutronCapture)

set(CMAKE_BUILD_TYPE RelWithDebInfo)

find_package(remage)

add_executable(TestNeutronCapture main.cc HPGeTestStand.hh HPGeTestStand.cc IsotopeOutputScheme.cc
IsotopeOutputScheme.hh)
target_link_libraries(TestNeutronCapture PUBLIC RMG::remage)
124 changes: 124 additions & 0 deletions examples/06-NeutronCapture/GammaCascades/156Gd-100Entries.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,124 @@
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% listing of gamma cascades by MAURINA/GAMMOC
%%
%% GD155+N KRK-D, 3LO GAMMA-RAY STRENGTH FUNCTIONS
%% MATNUM: 15511
%% E_n= 5 keV
%%
%% date: 20221101 time: 212412.007/pgrabmayr
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% En: neutron energy [keV]
%% Ex: excitation energy [keV]
%% M: multiplicity of gamma cascade; M_max=20
%% Em: missing energy [keV]
%% :=0 ground state after capture
%% :=E_iso isomeric state after capture
%% & M=M_max :=E_left left out due too many gammas
%% Eg1: energy of first photon [keV]
%% Eg2: energy of 2nd photon [keV] ..... EgM
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% variables:
%% En/keV Ex/kev M Em/keV Eg1/keV Eg2/kev Eg3/kev etc.. EgM
%% format :
%% i6 i6 i3 i6 i6 i6 i6 i6 i6 i6 i6 i6
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5 8540 6 0 1895 1065 2233 1980 1278 89
5 8540 5 0 992 5156 1025 1278 89
5 8540 4 0 2765 2985 2701 89
5 8540 3 0 3360 3938 1242
5 8540 4 0 1800 3485 3166 89
5 8540 1 0 8540
5 8540 4 0 5310 1988 1153 89
5 8540 3 0 3262 2361 2917
5 8540 4 0 1700 5711 1040 89
5 8540 4 0 1618 3282 3551 89
5 8540 6 0 1516 3285 1805 780 1065 89
5 8540 4 0 4097 2252 1037 1154
5 8540 4 0 3600 3664 1187 89
5 8540 4 0 6359 939 1153 89
5 8540 5 0 2704 3576 940 1231 89
5 8540 5 0 1589 4711 1111 1040 89
5 8540 4 0 2732 3591 2128 89
5 8540 4 0 3620 2087 2744 89
5 8540 4 0 2979 3217 977 1367
5 8540 6 0 1928 1191 3910 1223 199 89
5 8540 4 0 3135 2227 2024 1154
5 8540 4 0 3611 2976 711 1242
5 8540 7 0 2320 1918 901 1389 736 1187 89
5 8540 3 0 4459 3992 89
5 8540 7 0 1337 3023 2352 580 960 199 89
5 8540 3 0 2500 2641 3399
5 8540 4 0 5705 565 2181 89
5 8540 5 0 2306 1640 3318 1187 89
5 8540 3 0 6239 2212 89
5 8540 4 0 5317 2935 199 89
5 8540 5 0 1792 3511 1939 1209 89
5 8540 7 0 1840 2958 765 1125 604 1159 89
5 8540 4 0 5441 1857 1153 89
5 8540 6 0 1356 2053 842 2922 1278 89
5 8540 7 0 2423 2431 1734 704 960 199 89
5 8540 3 0 3822 4629 89
5 8540 5 0 4878 1345 1075 1153 89
5 8540 5 0 3986 3278 988 199 89
5 8540 6 0 3901 2609 675 1067 199 89
5 8540 4 0 3287 3053 2111 89
5 8540 3 0 4254 4197 89
5 8540 5 0 3808 2629 1815 199 89
5 8540 5 0 2217 1396 3679 1159 89
5 8540 5 0 1073 4355 2824 199 89
5 8540 5 0 2100 2251 1998 2102 89
5 8540 3 0 2912 942 4686
5 8540 4 0 2298 4875 1278 89
5 8540 3 0 6034 2417 89
5 8540 4 0 2833 3855 698 1154
5 8540 5 0 3252 2173 1760 1266 89
5 8540 6 0 2682 1768 2156 567 1278 89
5 8540 4 0 1295 4192 2964 89
5 8540 6 0 5452 1043 534 263 1159 89
5 8540 4 0 1984 4435 2032 89
5 8540 5 0 2259 4053 861 1278 89
5 8540 4 0 2560 3722 2169 89
5 8540 5 0 2910 1971 2383 1187 89
5 8540 5 0 4593 1844 949 1065 89
5 8540 4 0 5404 1894 1153 89
5 8540 6 0 4208 1509 1565 970 199 89
5 8540 5 0 2864 2649 1898 1040 89
5 8540 5 0 4617 1953 694 1187 89
5 8540 4 0 3722 2793 486 1539
5 8540 7 0 2771 1472 2122 877 1010 199 89
5 8540 4 0 4131 3133 1187 89
5 8540 7 0 2752 1428 1495 1242 1335 199 89
5 8540 6 0 3119 1518 1643 893 1278 89
5 8540 7 0 2827 1275 1391 1578 1181 199 89
5 8540 6 0 1430 781 2105 1832 1150 1242
5 8540 5 0 2368 4069 949 1065 89
5 8540 6 0 3002 1668 2331 1251 199 89
5 8540 4 0 2263 1231 2076 2970
5 8540 4 0 3676 3544 1231 89
5 8540 5 0 3082 2850 1288 1231 89
5 8540 3 0 5109 1284 2147
5 8540 3 0 3614 3684 1242
5 8540 5 0 4121 2075 977 1278 89
5 8540 6 0 2928 3760 497 1067 199 89
5 8540 4 0 4300 2920 1231 89
5 8540 5 0 1874 3372 2036 1169 89
5 8540 5 0 3941 2587 692 1231 89
5 8540 5 0 2498 2964 1924 1065 89
5 8540 2 0 6393 2147
5 8540 4 0 1970 3725 1691 1154
5 8540 3 0 6218 2233 89
5 8540 5 0 4378 879 2041 1153 89
5 8540 3 0 4760 1334 2446
5 8540 7 0 4015 2495 407 1038 297 199 89
5 8540 4 0 3966 3254 1231 89
5 8540 5 0 5017 1743 532 1159 89
5 8540 4 0 1502 5718 1231 89
5 8540 2 0 4525 4015
5 8540 3 0 2571 3646 2323
5 8540 6 0 1595 1430 2344 1050 2032 89
5 8540 4 0 3475 3698 1278 89
5 8540 5 0 5540 1107 1605 199 89
5 8540 5 0 3397 2873 1141 1040 89
5 8540 8 0 1504 1512 3349 767 823 297 199 89
5 8540 4 0 4060 1143 3248 89
5 8540 4 0 1739 5481 1231 89
124 changes: 124 additions & 0 deletions examples/06-NeutronCapture/GammaCascades/158Gd-100Entries.txt
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,124 @@
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% listing of gamma cascades by MAURINA/GAMMOC
%%
%% GD157+N KRK-D 3LO GAMMA-RAY STRENGTH (PAPER)
%% MATNUM: 15711
%% E_n= 5 keV
%%
%% date: 20221218 time: 180456.105/pgrabmayr
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%$$%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% En: neutron energy [keV]
%% Ex: excitation energy [keV]
%% M: multiplicity of gamma cascade; M_max=20
%% Em: missing energy [keV]
%% :=0 ground state after capture
%% :=E_iso isomeric state after capture
%% & M=M_max :=E_left left out due too many gammas
%% Eg1: energy of first photon [keV]
%% Eg2: energy of 2nd photon [keV] ..... EgM
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% variables:
%% En/keV Ex/kev M Em/keV Eg1/keV Eg2/kev Eg3/kev etc.. EgM
%% format :
%% i6 i6 i3 i6 i6 i6 i6 i6 i6 i6 i6 i6
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
5 7941 4 0 1830 2292 3739 80
5 7941 3 0 2443 5418 80
5 7941 5 0 3856 2107 936 962 80
5 7941 4 0 4174 3506 181 80
5 7941 5 0 2329 3345 1290 897 80
5 7941 4 0 4112 2787 962 80
5 7941 2 0 5657 2284
5 7941 4 0 1868 3075 2918 80
5 7941 5 0 4050 1771 1078 962 80
5 7941 3 0 6917 944 80
5 7941 5 0 2026 3967 924 944 80
5 7941 4 0 3559 3358 944 80
5 7941 3 0 3700 3264 977
5 7941 5 0 1380 2895 1730 1856 80
5 7941 4 0 2047 4631 1183 80
5 7941 3 0 5299 2562 80
5 7941 4 0 2829 2787 1348 977
5 7941 4 0 2154 1296 1992 2499
5 7941 4 0 3389 2759 816 977
5 7941 4 0 2298 4203 1360 80
5 7941 3 0 3574 3104 1263
5 7941 5 0 1108 5014 777 962 80
5 7941 5 0 2931 3062 924 944 80
5 7941 5 0 3905 2243 751 962 80
5 7941 4 0 2423 4494 944 80
5 7941 5 0 2874 2742 1301 944 80
5 7941 4 0 1830 3901 2130 80
5 7941 5 0 1264 5100 1316 181 80
5 7941 5 0 3627 2911 361 962 80
5 7941 5 0 2307 3309 1301 944 80
5 7941 5 0 2063 3116 1738 944 80
5 7941 5 0 2831 1939 2147 944 80
5 7941 5 0 1774 4809 1097 181 80
5 7941 3 0 2827 4137 977
5 7941 3 0 3166 2880 1895
5 7941 6 0 3618 2504 660 898 181 80
5 7941 4 0 2502 3786 466 1187
5 7941 4 0 3139 3166 1556 80
5 7941 2 0 5601 2340
5 7941 4 0 1493 2501 3867 80
5 7941 4 0 3311 3367 1183 80
5 7941 4 0 2970 3335 1556 80
5 7941 6 0 1471 1680 2772 1757 181 80
5 7941 5 0 2770 1285 2446 1360 80
5 7941 4 0 5062 2618 181 80
5 7941 5 0 1017 2320 2555 1072 977
5 7941 6 0 1884 4103 540 1153 181 80
5 7941 4 0 3438 3526 897 80
5 7941 4 0 3897 3020 944 80
5 7941 4 0 1094 2447 4320 80
5 7941 5 0 2051 2014 2852 944 80
5 7941 6 0 1895 1904 2983 898 181 80
5 7941 4 0 1794 2136 2751 1260
5 7941 5 0 1767 3728 1469 897 80
5 7941 5 0 2124 2025 1851 1861 80
5 7941 5 0 3069 3011 884 897 80
5 7941 5 0 2181 4601 898 181 80
5 7941 3 0 5977 1884 80
5 7941 5 0 3598 2121 819 1323 80
5 7941 4 0 3409 3555 897 80
5 7941 5 0 1517 4304 706 1334 80
5 7941 3 0 5681 2180 80
5 7941 4 0 5726 1954 181 80
5 7941 4 0 4120 2779 962 80
5 7941 4 0 5657 2023 181 80
5 7941 4 0 3934 2983 944 80
5 7941 5 0 3108 2855 1001 897 80
5 7941 5 0 3706 2832 1142 181 80
5 7941 4 0 3894 2784 1183 80
5 7941 5 0 1749 648 3608 1856 80
5 7941 4 0 2950 3130 884 977
5 7941 4 0 4062 1901 1898 80
5 7941 5 0 3072 1441 2451 897 80
5 7941 6 0 2695 3385 819 781 181 80
5 7941 4 0 1794 5170 897 80
5 7941 3 0 2423 5438 80
5 7941 5 0 3227 3555 898 181 80
5 7941 4 0 2218 4320 1323 80
5 7941 4 0 2948 4016 897 80
5 7941 5 0 2708 1356 2463 1334 80
5 7941 5 0 2318 3804 795 944 80
5 7941 5 0 2295 4232 372 962 80
5 7941 5 0 2580 3708 611 962 80
5 7941 4 0 3801 3098 962 80
5 7941 5 0 3808 3091 781 181 80
5 7941 6 0 2316 3806 643 915 181 80
5 7941 6 0 3384 2603 912 781 181 80
5 7941 4 0 2389 4289 1183 80
5 7941 2 0 7861 80
5 7941 5 0 2015 4474 1191 181 80
5 7941 3 0 3954 3010 977
5 7941 4 0 3229 2760 1872 80
5 7941 5 0 2317 1880 3483 181 80
5 7941 6 0 1842 2133 2924 781 181 80
5 7941 5 0 2697 2977 1243 944 80
5 7941 4 0 2326 4212 1323 80
5 7941 3 0 2810 3868 1263
5 7941 4 0 2869 3436 1556 80
5 7941 5 0 3266 2727 906 962 80
5 7941 4 0 3883 2655 1323 80
Loading